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Some Significant Developments in 
Aerodynamics Since 1946 


The First Daniel and Florence Guggenheim Memorial Lecture 


THEODORE 


[ IS A GREAT HONOR and pleasure for me to give the 
first Daniel and Florence Guggenheim Memorial 
Lecture before this distinguished international gather- 
ing. The name Guggenheim will be connected for- 
ever in the history of aviation with the development of 
the aeronautical sciences. In fact, the Daniel Guggen- 
heim Fund for the Promotion of Aeronautics, which 
Daniel Guggenheim established in 1926, gave a great 
impetus to education and research at a number of 
American universities. I am personally indebted to 
the Guggenheim family by the fact that my change of 
continents from Europe to America originated with 
Daniel Guggenheim and Robert A. Millikan, world- 
famous physicist and then head of the California In- 
stitute of Technology. The story, as I learned it from 
R. A. Millikan, was that, in 1926, as the Daniel Gug- 
genheim Fund started to distribute grants to several 
colleges for the purpose of forming graduate schools for 
aeronautics, Millikan undertook a trip to Long Island 
where Mr. Guggenheim lived. During his visit with 
Mr. Guggenheim, Millikan told him that he would make 
the greatest mistake of his life if a sizable chunk of the 
money of the Fund did not go to California. California 
would be, Millikan insisted, in the near future, the most 
important center of American aircraft production, due 
to climatic conditions and availability of space for fly- 
ing establishments. Mr. Guggenheim answered that 
he would provide the funds for a graduate school if 
Millikan would bring to Pasadena from Europe some- 
one familiar with and active in aeronautical research, 
and especially familiar with the theoretical side of 
research. In this way I received a call to come to the 
United States. I was told later that Robert A. Milli- 
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kan said, ‘First I aimed at Prandtl—then I settled on 
Karman.” Prandtl was, without doubt, the leading gen- 
ius in the early development of modern aerodynamics. 

Daniel Guggenheim was a man of great vision and 
lively interest in progress in aeronautics. I remember 
that after I went through the United States visiting all 
the places where Guggenheim funds were invested for 
the promotion of aeronautical science, he asked me 
what he should do in addition to supporting the schools, 
so that aeronautical science would flourish in the United 
States as it had grown and brought ripe fruits in several 
European countries during the first two decades of our 
century. I suggested that it would be extremely de- 
sirable to have a kind of Handbook of Aerodynamic 
Theory that would lay down the present state of funda- 
mental aeronautical knowledge, so that young American 
scholars would know from which point to start. Then 
I said jokingly, “‘In Europe it was very helpful that 
young scientists could meet in quiet cafés for informal 
discussions. Unfortunately there is no equivalent in- 
stitution in the United States.”’ Mr. Guggenheim an- 
swered, ‘All right, I shall provide the money for the 
book."’ (As a matter of fact, this was the origin of 
Durand’s ‘‘Aerodynamic Theory,” a collective work 
based on international authorship which has had a 
great influence on the development of American scien- 
tific literature in the field of Aerodynamics.) ‘‘How- 
ever,’ he continued, “I will not go into the cafeteria 
business—even for science’s sake!”’ 

Daniel Guggenheim, born in Philadelphia in 1856, 
was of the second generation in the leadership of the 
large Guggenheim enterprises in mining and affiliated 
industries. His son, our great friend Harry, once told 
me jokingly that, in his grandfather's time, a really good 
man, in order to be recognized as such, had to create a 
great fortune. When his fether was in business, the 
main duty was to broaden the scope and to maintain 
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the level of the enterprises. “‘Of me,” he said, “they 
only expect that I spend the money usefully and grace- 
fully.” 

I am sure that this formulation of Harry F. Guggen- 
heim’s activities is extremely modest. He has an im- 
pressive record of public service; he was a successful 
ambassador of the United States to Cuba. There is no 
doubt that his interests in aeronautics and his personal 
experience in aviation, especially as a naval aviator dur- 
ing the first World War, were the decisive factor for 
Daniel Guggenheim’s donations being directed toward 
the promotion of the aeronautical sciences. Also, the 
organization of the Safe Aircraft Competition in 1930, 
the active support of Professor Goddard's early efforts 
in rocket research, and the establishment of the Daniel 
and Florence Guggenheim Jet Propulsion Centers at 
C.I.T. and Princeton University, the Institute of Flight 
Structures at Columbia University, and the Center for 
Aviation Health and Safety at Harvard are far beyond 
the qualification ‘‘spend the money usefully and grace- 
fully.”’ They are proofs of an extraordinary vision for 
the kind of education and research needed for the prog- 
ress of aviation. And last but not least, this ‘First In- 
ternational Congress for the Aeronautical Sciences” 
owes its birth not only to the donation of funds by the 
Daniel and Florence Guggenheim Foundation, but in 
the first place to the initiative and desire of Harry F. 
Guggenheim to spend the money for the organization of 
an international congress open to scientists and engi- 
neers of all the nations of the globe which have an asso- 
ciation or group devoted to the aeronautical sciences. 

I shali now proceed to the technical part of my lec- 
ture. 

For my review of progress in aerodynamics I would 
like to choose as a point of departure my Wright Broth- 
ers lecture delivered on December 17, 1946. I wrote at 
that time: 

“T believe we have now arrived at the stage where 
knowledge of supersonic aerodynamics should be con- 
sidered by the aeronautical engineer as a necessary 
prerequisite to his art. This branch of aerodynamics 
should cease to be a collection of mathematical and 
half-digested, isolated, experimental results. The aero- 
nautical engineer should start to get the same feeling 
for the facts of supersonic flight as he acquired in the 
domain of subsonic velocities by a long process of theo- 
retical study, experimental research, and flight experi- 
ence.” 

I have the impression that this goal has generally 
been achieved. There has been rapid progress in the 
theory of supersonic wings, bodies, and wing-body com- 
binations exposed to supersonic flow, and this progress 
has been reported in the scientific literature. Experi- 
mental work has been carried out in newly created fa- 
cilities in governmental and industrial laboratories. 
Finally, advanced students have had improved train- 
ing at our aeronautical schools. Asa result of these de- 
velopments, a number of engineering firms now have at 
their disposal large staffs with the necessary grasp and 
appreciation of the main features of supersonic aerody- 


namics. As a matter of fact, some of our supersonic 
bombers have been designed on the basis of more de- 
tailed aerodynamic calculations than was possible in 
the case of the best subsonic aircraft. 

However, the honeymoon was short. Soon, new 
problems were facing the aeronautical engineer, who 
nowadays is pleased to call himself a missile engineer 
or even a space technologist. 

The problem of ballistic missiles led us to the range 
of hypersonic speeds, to speeds which are not only com- 
parable with sound velocity but are greater by an order 
of magnitude than the velocity of sound. At first 
sight, the hypersonic range introduces certain simplifi- 
cations into the flow problem; as a matter of fact, it 
was pointed out as early as 1931 by P. S. Epstein, in a 
paper devoted to the problem of the drag of artillery 
projectiles, that for very large Mach Numbers the 
classical Newtonian law of air resistance becomes 
valid. At that time, and even at the time of the mem- 
orable Volta Congress for High Speed, held at Rome, 
in 1935, high Mach Numbers belonged to the realm of 
academic speculation. 

The simplification introduced by very high Mach 
Numbers is largely overbalanced by the complications 
due to the high temperatures caused by shock and fric- 
tion. The production of heat, an annoyance in flight at 
moderate Mach Numbers, becomes a major problem at 
hypersonic speeds. Furthermore, and as a new com- 
plication, one has to take into account the chemical 
changes in the air, such as dissociation and recombina- 
tion. No longer are we dealing with pure aerodynam- 
ics, nor aerothermodynamics; fluid mechanics must now 
be combined not only with thermodynamics, but also 
with chemistry. 

I suggested the term aerothermochemistry for the com- 
bination of these three disciplines. I had mainly in 
mind problems related to combustion, like flame the- 
ory, the theory of quenching, and the like. Hyper- 
sonics has now made it necessary to consider chemical 
reactions which occur without having been planned by 
the chemical engineer. 

Recently, another combination of various disciplines 
has attracted considerable attention—fluid mechanics 
and the theory of the electromagnetic field. Some as- 
pects of the motion of conducting liquids acted upon by 
electromagnetic forces were investigated several dec- 
ades ago. I also remember that Albert Einstein gave 
some thought to a thoroughly practical problem in this 
field: he proposed the design of a refrigerator in which 
the coolant, for example a liquid metal, would be kept 
in circulation by an imposed electromagnetic field. He 
wanted to avoid the use of machinery which needs lubri- 
cation. However, the main interest in ‘“‘magnetofluid- 
mechanics” arose from celestial problems, such as the 
structure and the motions of galaxies, wave motions 
and turbulence in cosmic systems and the like. Ques- 
tions of space flight and problems related to the possible 
utilization of thermonuclear reactions led to increasing 
interest in “‘plasmadynamics.’’ I believe that a sys- 


tematic nomenclature and classification of these new 
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branches of combined fluidmechanics, electromagnetic 
theory and thermodynamics is still lacking. Never- 
theless, I believe that in a review of advances of aero- 
dynamics, they should at least be mentioned. 

A detailed list of references is appended hereto. I 
now want to restrict myself to making some comments 
on some accomplishments in aerodynamics which, I 
believe, are most significant. One general remark may 
precede this. It is remarkable how far one can go with 
the so-called linear theory of supersonic flow, and how 
many useful conclusions for aircraft design can be de- 
rived from this simplified theory which, after all, is 
nothing else than ‘‘acoustics’’—-1.e., it is based on the as- 
sumption that the flow is composed of a uniform parallel 
flow and a perturbation flow of small magnitude. 


(1) THEORY 


In a paper published in 1946, A. E. Puckett applied the 
method of singularities, which I had introduced jointly 
with N. B. Moore, into the first theory of supersonic 
flow around axially symmetric bodies. We used su- 
personic sources distributed along the axis of the body. 
H. S. Tsien introduced doublets (dipoles) in order to 
include lift in the theory for the case of a body of 
revolution with angle of attack. These supersonic 
sources and doublets are formulated in such a way that 
their effect is restricted to the interior of a Mach cone, 
the apex of which is the point of singularity. Thus, the 
solution is essentially identical with the solution of a 
time-dependent two-dimensional acoustic problem, 
where the length coordinate in the flow direction re- 
places the time coordinate of the acoustic phenomena. 

Subsonic wings are mostly calculated by means of the 
concept of the lifting line. Only in a few cases did it 
prove possible to solve the problem of the lifting sur- 
face. However, for the supersonic case, the situation is 
more favorable. Especially if both the leading and 
trailing edges are of the supersonic type—i.e., if the 
components of the flight speed normal to the edges are 
greater than sound velocity—the local slope of the wing 
determines directly the required distribution of sources 
over the wing plan form. These singularities, together 
with the condition that in undisturbed flow the pressure 
is equal to the ambient pressure, completely determine 
the disturbance potential. 

Two complications occur if some portions of the per- 
imeter of the wing plan form are edges of the subsonic 
type: 

For the case of a subsonic leading edge, the flows on 
the upper and the lower surfaces near the edge are no 
longer independent. The fluid flows around the edge 
with subsonic velocity, and one obtains a leading-edge 
suction which has to be calculated. Furthermore, in 
the sector enclosed between leading edge and limiting 
Mach line, the pressures resulting from the upper and 
lower half spaces must be balanced. 

In the case of a subsonic trailing edge the Kutta-Jou- 
kowski condition has to be satisfied; in other words, the 
velocity components normal to the edge on the upper 
and the lower surface must be equalized. This makes 


Fic. 1. Plan form of elliptic supersonic wing. 


it necessary to compute and compensate the velocities 
induced by the trailing vortices in certain sectors of the 
plan form. 

Fig. 1 is reproduced from an excellent review of the 
supersonic wing theory published by R. T. Jones and 
Doris Cohen in Volume VII of the Princeton series 
“High Speed Aerodynamics and Jet Propulsion.’ The 
Figure shows, for the case of an elliptic plan form, the 
various domains which have to be treated in different 
ways. Sector I has a supersonic leading edge and is 
not influenced by any other portion of the wing plan 
form. Hence, the Puckett method can be applied di- 
rectly. Sector II has subsonic leading edges; Sectors 
IV, V, VI, and VII are evidently influenced by the trail- 
ing vortices leaving the subsonic trailing edges FE and 
CD. The vortex wake can be built up by superposition of 
horseshoe vortices (Fig. 2). The Figure also shows the 
domain in which the flow is influenced by such a vortex. 

The problem of the partially subsonic leading edge 
was resolved in an ingenious way by a method appar- 
ently independently proposed by J. C. Evvard in the 
United States and E. A. Krasilshchikova in the USSR. 

Before that, H. J. Stewart treated the special case of 
the delta wing with subsonic straight leading edges by 
means of the method of conical flows. This method, 
originally suggested by A. Busemann, became one of 
the most powerful methods of supersonic aerodynamics. 
It reduces the problem of three-dimensional flow to a 
two-dimensional problem which can be solved by means 
of Laplace’s equation. 

For example, R. Legendre of France and M. C. 
Adams, C. E. Brown, W. H. Michael, and R. H. Ed- 
wards of the United States treated the case of delta 
wings with flow separation at the apex as an example of 
conical flow. 

The linearized wing theory made it possible for R. T. 
Jones to arrive at a series of important results concern- 
ing the minimum drag of supersonic wings. He made 
very ingenious use of the concepts of the reverse flow 
and the combined flow. In earlier studies I had found 
that the drag of a thin wing due to thickness remains 
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Fic. 2. Horseshoe vortex in supersonic flow. 


unchanged if we reverse the flight direction. If we 
superpose both perturbations, which correspond to the 
two opposite flight directions, on one parallel uniform 
flow we arrive at a so-called combined flow. Using 
this concept and the fact that the drag is independent 
of the flight direction, R. T. Jones has shown that the 
drag of a symmetric flat wing of given plan form reaches 
a minimum if the thickness is distributed over the plan 
form in such a way that the drag per unit volume is 
equal for every wing element. 

Jones also investigated the condition for the minimum 
of the drag caused by lift. He found that the optimum 
lift distribution over a given plan form is such that the 
downwash is constant over the planform. It is evident 
that this result is a kind of generalization of Munk’s 
rule for the optimum distribution of lift along a lifting 
line, which leads to the elliptical distribution in classical 
subsonic wing theory. 


"DEVELOPMENT OF CONICAL CAMBER 
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Many special cases have been investigated from the 
viewpoint of minimizing drag. Thus, Doris Cohen has 
shown that, for triangular wings with subsonic edges, 
the minimum drag due to lift is obtained by using nega- 
tive spanwise camber over most of the plan form. It is 
interesting that the idea of conical chamber was used in 
the design of the wing of the Convair B-58 bomber (also 
called the Hustler). This development is attributed to 
C. F. Hall (Ames Laboratory). It is seen from the ex- 
perimental results represented on Fig. 3 that the coni- 
cal camber produces leading-edge suction and thus es- 
sentially reduces drag. 


(II) SLENDER WINGS AND AIRCRAFT 


With the increase of the speeds of flight, aircraft in 
general became more and more slender; large aspect 
ratios disappeared, and present-day aircraft mostly con- 
sist of long bodies and wings of small aspect ratios. 
For such cases, the application of the linearized theory 
had to be revised. Considerable confusion had been 
caused by formal application of the Prandtl-Glauert 
rules for compressibility effects to bodies of revolu- 
tion and slender wings and bodies in general. This 
question was clarified by Goethert and others in a satis- 
factory way. It was shown that the similarity theory 
of compressible flow which was correct for wings of large 
aspect ratio gives false results for slender bodies. It 
was also shown by M. J. Lighthill and G. N. Ward that 
the computation of the pressure based on the axial com- 
ponent of the disturbance velocity is inexact, since, for 
example, in the case of axisymmetric bodies, the influ- 
ence of the radial velocity is often larger. Hence, it 
became necessary to describe the flow by the appropri- 
ate singularities and to compute the pressure including 
the quadratic terms of the disturbance velocities. In 
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Fic. 3. Development of conical camber. 
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order to simulate the flow around bodies which are 
neither thin wings nor axisymmetric bodies, multipoles 
have to be introduced. 

A great simplification was achieved by extension of 
Munk’s classical airship theory to slender bodies in 
general both for subsonic or supersonic speeds by R. T. 
Jones. This theory assumes that the flow around a 
slender body (wing or aircraft) can be considered as a 
superposition of a parallel uniform flow and a sequence 
of two-dimensional flows in planes perpendicular to the 
flight direction as shown in Fig. 4. To some extent 
this is the opposite of Prandtl’s classical assumption for 
his wing theory. Asa matter of fact, for wings of large 
aspect ratios it was assumed that the flow in planes per- 
pendicular to the wing axis can be identified with a 
two-dimensional flow with circulation around the wing 
section. 

The slender-body theory in general gives a better 
approximation to the pressure, forces, and moments 
acting on the system than the pure linearized theory. 


(III) INTERFERENCE EFFECTS 


We shall now consider interference effects between 
the wing and the body structures of the airplane which 
are necessary to carry payload, passengers, or fuel. The 
most significant example of favorable interference is 
the arrangement popularly known as the area rule. W. 
D. Hayes has shown in his doctoral thesis at the Cali- 
fornia Institute of Technology that the resulting flow 
around a system consisting of bodies and thin wings can, 
at large distances, be represented as originating from 
singularities distributed along the axis. At Mach 
Numbers near one the drag is that caused by a single 
equivalent body of revolution. Now we know—for 
example, according to the theories of W. Haack and W. 
R. Sears—the shape of the most favorable bodies of 
revolution as far as minimum drag is concerned. Thus, 
the components of the aircraft can be arranged so that, 
at least for a given Mach number, the equivalent body 
of revolution approaches the shape corresponding to 
minimum drag. This principle of equivalence was 
clearly stated by Hayes and was extended by kK. Oswa- 
titsch to the nonlinear case; the extension is important 
because of the limitations of linearized theory near sonic 
speed. The principle leads, for example, to the con- 
clusion that in the case of nacelles arranged near the 
fuselage the additional drag of the nacelles can be com- 
pensated by a reduction of the diameter of the fuselage 
in the appropriate section. We obtain in this way the 
shape of the fuselage known popularly as the ‘‘Marilyn 
Monroe body.” The name “‘area rule”’ originates from 
the fact that for a flow near Mach one the section 
of the equivalent body of revolution is simply equal to 
the sum of the areas cut out by a plane laid through the 
section considered. The agreement of these theo- 


retical ideas with drag measurements was conclusively 
shown by Richard Whitcomb. 

Another class of interference phenomena makes it 
possible to create favorable interference by reflection of 
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Fic. 4. Wing of low aspect ratio. 
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Fic. 5. Busemann biplane. 


Wedge-wing lifting configuration. 


Fic. 6. 


compression and expansion waves on components of a 
lifting system. 

The oldest example for such a procedure is the so- 
called Busemann biplane (see Fig. 5). In this case we 
may transfer pressure from one wing to the rear portion 
of the other wing (for example, from AB to EF and from 
DE to BC), so that the pressures acting on the forward 
and rear portions of the inner surface are balanced. 
Thus, the drag due to thickness effect can be eliminated, 
and we gain volume without additional drag. 

An interesting example of similar favorable inter- 
ference effects is represented by a combination of a hori- 
zontal delta wing and a vertical surface of wedge shape 
arranged below the wing parallel to the flight direction 
(Fig. 6). One can transfer the pressure created on the 
surface of the vertical wedge to the lower surface of the 
wing. Thus the total lift is increased, and one can 
show that the resulting lift/drag ratio is more favorable 
than for the case of a simple horizontal wing that would 
furnish the same lift. 


(IV) AEROELASTIC THEORIES 


Concerning aeroelastic problems like flutter at high 
flying speeds, we have to mention the so-called “piston 
theory” initiated by M. J. Lighthill and worked out by 
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Holt Ashley and others. This theory assumes that the 
magnitude of the reaction of a high-speed flow on a 
wing surface can be approximated at every instant by 
considering the one-dimensional motion of an air col- 
umn under the action of a moving piston. It seems 
that for Mach Numbers superior to = 1.5 this ap- 
proximate theory furnishes usable results and greatly 
simplifies the computations. 


(V) NONLINEAR THEORIES—HIGHER-ORDER 
APPROXIMATIONS 


It is not possible to enter into exact discussion of the 
highly intricate mathematical problems connected with 
the development of higher-order approximations for a 
solution of the equations of supersonic flow even for the 
cases of plane and axisymmetric motions. However, 
the main trends of the pertinent approaches will be in- 
dicated. 

For the purpose of a two-dimensional airfoil theory, 
Ackeret’s fundamental and simple results correspond to 
a linearized theory. Ata relatively early date (1935), 
Busemann showed the possibility of obtaining a correc- 
tion of second order in the perturbance velocities, and 
in 1948 K. Friedrichs provided a means for obtaining 
shock shapes to the same degree of approximation. 
Therefore, a rather complete and accurate theoretical 
treatment of the two-dimensional airfoil problem is 
available. This applies as long as the effects of entropy 
changes are relatively unimportant. For practical air- 
foil shapes the condition seems to be satisfied up to 
Mach Numbers of 5 or 6. Above this limit, tempera- 
ture and real-gas effects make the application of ideal- 
gas laws illusory anyway. 

The attempts to develop higher approximations for 
three-dimensional flows, at least for specific cases, are 
reviewed by M. J. Lighthill in his article ‘““‘Higher Ap- 
proximations” published in Section E, Vol. VI of the 
Princeton series. 

In the linearized theory, simplifications are intro- 
duced in three aspects of the problem: first the equa- 
tion of motion is only approximately correct; second, 
the boundary conditions at the body are simplified; 
and third, the pressure is computed on the basis of a sim- 
plified relation between pressure and velocity compo- 
nents. The simplest idea is to keep the linearized 
character of the equation of motion, but to correct the 
boundary conditions and the pressure-velocity rela- 
tions. One is tempted to call such theories ‘‘hybrid”’ 
theories. In general, the procedure leads to results 
which are in better agreement with exact solutions and 
or the experimental evidence, but it can work also in 
the opposite sense. 

J. B. Broderick and M. D. Van Dyke worked out 
second-order solutions on the basis of this general idea; 
Van Dyke uses the first-order linearized solutions as a 
starting point, whereas Broderick starts from the 
slender-body solution. 

One of the main weaknesses of the linearized theory 
is the fact that it describes the flow quite inexactly at a 
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certain distance from the body or wing In fact, ac- 
cording to the linearized theory the characteristic lines 
are always straight lines, and their inclination corre- 
sponds to the Mach Number of the undisturbed flow. 
G. B. Whitham has provided a relatively simple pro- 
cedure for improving the flow-field representation. He 
assumes that the flow variables predicted by the linear- 
ized theory are satisfactorily exact as far as their mag- 
nitude is concerned, but must be relocated spatially ac- 
cording to a corrected Mach wave shape. He uses the 
first-order corrections obtained from the linearized the- 
ory for the local sonic velocity and the values for the to- 
tal velocity obtained from the same theory, and uses 
these data to construct corresponding Mach waves, in- 
cluding eventual shock surfaces. 

The Whitham theory, in which the positions of the 
Mach waves are changed without changing the values 
of the flow variables, can be considered as a special case 
of the general method or “‘technique’”’ introduced by 
Lighthill, which he calls a “technique for rendering ap- 
proximate solutions to physical problems uniformly 
valid.” In this technique—which is based on _ the 
methods used by H. Poincaré in his investigations on 
celestial mechanics—in addition to developing the de- 
pendent unknowns in series of successive terms as 
functions of some parameters, the independent vari- 
ables or a number of independent variables are also de- 
veloped as functions of the same parameters. By this 
trick, in domains where the dependent variables have a 
singular behavior as functions of the independent varia- 
bles, they become easy to handle as functions of the new 
parameter which is used for a series development of the 
independent variable itself. This technique has proved 
especially useful in some boundary-layer problems which 
we shall mention later. As far as applications to non- 
viscous problems are concerned, the treatment of the 
mixed flow (hyperbolic and elliptic) about a sharp- 
edged conical wing by R. Vaglio-Laurin is worth men- 
tioning. 

In connection with higher-order approximations, we 
have to mention also the so-called linearized method of 
characteristics used by A. Ferri to build up three-di- 
mensional flow fields from nonlinear two-dimensional 
flow fields which are known by application of classical 
characteristic methods as used in two-dimensional su- 
personic aerodynamics since the fundamental works of 
Prandtl. Ferri finds that the characteristic surfaces of 
the three-dimensional flow field can be approximated by 
the envelopes of the characteristic surfaces of the sys- 
tem of two-dimensional flow fields. Hence the analy- 
sis of the three-dimensional flow is reduced to that of a 
flow field with given characteristic surfaces. 


(VI) TRANSONIC FLow 


The transonic case is characterized by the fact that 
the difference between flight velocity LU’ and sound veloc- 
ity c; i.e., the quantity (U — c); is small in compari- 
son with either U’ orc. It is known that, at least in the 
case of a two-dimensional flow, for the disturbance po- 
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tential we arrive at a relatively simple nonlinear equa- 
tion called Tricomi’s equation: 


(1 — + (0%¢/dy*) = 
(y + 1) (M./c.) (O¢/dx) (0%¢/ 0x?) 


where ./., denotes the Mach Number of the main flow, 
C., the sound velocity at infinity, and y the ratio of the 
specific heats. The nonlinear term of this equation 
contains the product (0¢/0x) (0°¢/0x?); Oswatitsch 
and Keune proposed to replace the derivative of the u 
component of the horizontal velocity, 0°¢/Ox*, by a 
value which is taken as independent of the coordinate 
y, but may vary stepwise as we proceed along the flow.* 
They arrived in this way at very useful results for com- 
putation of transonic flows around wings and bodies. 

To a certain extent, useful results have been achieved 
for transonic flow calculations by the so-called Karman- 
Tsien method, which introduced a “‘hypothetic fluid”’ 
with a simplified pressure-density relation, valid in a 
limited Mach number range. This idea was further de- 
veloped and improved by S. Tomotika and K. Tamada. 

On the other hand, theoretical studies concerning 
exact solutions of the Tricomi equation—such as the 
problem of uniqueness of the solution and the questions: 
in which cases do we obtain continuous solutions and 
in which cases do shocks appear?—have not progressed 
sufliciently in the last ten years. 

In an attempt to obtain approximate but useful solu- 
tions of transonic flow problems, one replaces the dif- 
ferential equation for the disturbance potential by a 
system of difference equations. K. Friedrichs and 
his collaborators used this method; they put the prob- 
lem into a computing machine. If the problems are 
properly formulated for the computing device, one can 
obtain fair approximations even for flows containing 
discontinuities (shocks) without an elaborate discus- 
sion concerning the existence or nonexistence of con- 
tinuous solutions. Such formulations can also be 
used in the case of detached shocks in hypersonic flow 
around blunt bodies. 


(VII) Hypersonic FLow 


The hypersonic speed range is characterized by the 
fact that the sonic velocity ¢ is small compared with the 
flight velocity U’. The interest in this speed range was 
recently very much enhanced by the problems con- 
nected with missile design, especially with the design of 
re-entry vehicles and nose cones. 

If we first consider the flow against an inclined sur- 
face, for example, the case of a wedge with small apex 
angle, the theory of inviscid fluids predicts attached 
and straight shock waves which slightly deflect the 
flow so that the streamlines practically become parallel 
to the wedge surfaces. Hence the flow picture is very 
similar to that assumed in the model treated by Isaac 
Newton. The essential difference is that according to 
Newton’s model the fluid particles move in parallel 
straight lines until they hit the surface and then are 


* This idea has been developed further by J. R. Spreiter. 


deflected into a motion along the body surface, whereas 
according to the ideal-flow theory they are deflected at 
the shock surface which, however, lies very near to the 
body surface. Thus, for a small angle of inclination, 
both for the case of the two-dimensional wedge and that 
of the axisymmetrical cone, we obtain the result that 
the pressure acting on the surface element is given by 


pb = (p/2)U* - 2 sin? B 


where p is the density in the undisturbed fluid, U the 
velocity of the same, and 8 the angle of inclination of 
the surface. 

We may call a flow of this nature a Newtonian flow. 
The actual evidence is, of course, strongly modified by 
two factors. First, the presence of the solid surface 
produces viscous effects—i.e., a boundary layer; and 
second, no mathematically exact sharp edge exists, and, 
therefore, one always obtains a kind of detached shock 
with very large curvature near the front portion. The 
real picture of the flow therefore looks more like the one 
for the case of a body with a blunt edge, where the 
“equivalent body”’ corresponds in a broad sense to the 
domain occupied by the boundary layer. This de- 
scription even applies to the flow observed in the case 
where the undisturbed flow is parallel to one surface of 
the unsymmetric wedge (see Fig. 7). Furthermore, the 
large curvature of the detached shock wave introduces 
vorticity in the inviscid domain of the flow between the 
shock and the boundary layer. 

Nevertheless, the rules of Newtonian flow can be ac- 
cepted as fair approximations for two- and _ three-di- 
mensional flows. If the body surface is curved, the 
flow in the inviscid region is necessarily curved and a 
correction for the Newtonian law can be obtained by 
computing the effect of ‘‘centrifugal forces’’ between 
the shock and the body surface. Such a correction was 
suggested by A. Busemann at an early date (//andbuch 
der Naturwissenschaften, 1933) and recently further 
developed by H. R. Ivey and R. R. Morisette, and fur- 
ther by H. R. Ivey, E. B. Klunker, and E. N. Bowen. 

If we consider a blunt body—for example, a cylindri- 
cal or spherical body exposed to a flow hitting the body 
normal to its surface-—two considerations have to be 
introduced. 

First, the value of the pressure at the stagnation 
point is, strictly speaking, unknown, since it is influ- 
enced by the thermodynamic process. This was un- 
known to Newton, as is evident if we consider that he 
computed the sound velocity, for example, assuming an 
isothermal compression process. 

In fact, there are in addition changes in the distribu- 
tion of the thermal energy between the several kinds of 
degrees of freedom in the gas molecules; at the high 
temperatures produced by the shock we certainly will 
encounter dissociation and recombination of atoms. 
Heat-transfer effects will modify the temperature and 
pressure distribution between the shock and the sur- 
face. 

The second consideration, which further complicates 
the problem, is the fact that we have a mixed-flow prob- 
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lem. Whereas, in the case of the sharp wedge or ogive, 
the flow in the inviscid region was supersonic on both 
sides of the shock wave, we now have downstream of 
the shock a subsonic region followed by a supersonic do- 
main. 

Nevertheless, if we consider the general type of flow 
shown schematically in Fig. 8, we can see that, at least 
in the neighborhood of the stagnation region, the fate of 
the individual gas particles is not very different from 
that assumed in Newton's model. The shock surface 
approximately follows the body surface and the individ- 
ual streamlines show that the particles, broadly speak- 
ing, are deflected from their initial direction in the di- 
rection along the body surface. Correspondingly, we 
obtain a quasi-Newtonian pressure distribution in the 
sense that the local pressure coefficient C, is given ap- 
proximately by the formula 


= Gis mar sin? B 


where C, maz is the local pressure coefficient at the stag- 
nation point and 8 is the local inclination of the surface 
relative to the flow direction. 

Beyond this quasi-Newtonian pressure distribution, 
attempts to determine the exact flow conditions in the 
region between the shock wave and the blunt body were 
not too successful until recently. There are many un- 
certainties: first, the location of the shock wave is un- 
certain; second, the sonic line can be determined only 
by a kind of iteration process, so that the subsonic and 
supersonic flow regions really match. Even if we 
stick to ideal-gas laws, we cannot consider the flow as a 
potential flow because of the vorticity produced in the 
shock region. The best approximation is that of con- 
stant vorticity along streamlines in the two-dimensional 
case and constant-vorticity flow in the axially sym- 
metric case. The computation usually requires the help 
of electronic computers. Among the iteration meth- 
ods, one used by the Soviet mathematicians Dorodnit- 
syn and Belotserkovski is reported to be quite success- 
ful. 

The flow beyond the sonic line can be calculated by 
the method of characteristics. Recent publications by 
H. M. Lieberstein and P. R. Garabedian, M. D. Van 
Dyke, and R. Vaglio-Laurin and A. Ferri deal with this 
problem. They consider the inverse case—to assume a 
shock shape and find the corresponding body shape. 

Fig. 9 shows experimental pressure measurements for 
the case of a body consisting of a long circular cylinder 
with a hemispherical frontpiece. The measurements 
were made at M/ = 7.7. 

Near the stagnation point the pressure distribution 
corresponds to the quasi-Newtonian rule; on the shoul- 
der, where the curvature of the meridian section 
changes discontinuously, apparently the rate of the 


Fic. 7 (left). Interferograms of flat plate model at various 
values of the Reynolds Number Re, based on the leading-edge 
thickness. M = 12.7. Top to bottom: (a) Re: = 177; (b) 
Re, = 578; (c) Re: = 1,510; (d) Re: = 2,480; (e) Re, = 4,050; 
(f) Re: = 7,480; (g) Re, = 15,010; (h) Rez = 35,200; (i) Rex — 
62,200 (Princeton Univ. Dept. of Aero. Engineering, Report No. 


326). 
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pressure drop can be described by the assumption of a 
Prandtl-Meyer expansion process. Then, however, one 
would expect that the pressure would rather quickly ap- 
proach the final value (which because of the expansion 
in the wake is slightly below the value corresponding to 
the undisturbed flow). However, observation shows a 
very slow decrease; evidently this corresponds to an ex- 
panding shock-wave surface over the cylindrical portion 
of the body. 

This process was made understandable by application 
of the analogy between an unsteady motion, in which 
the flow picture remains similar to itself, and a steady 
flow, in which the length coordinate in the flight direc- 
tion replaces the time coordinate of the unsteady proc- 
ess A characteristic example of such unsteady solu- 
tions is the well known solution which G. I. Taylor ob- 
tained during World War II for the problem of a violent 
spherical explosion. The propagation of a weak explo- 
sion with sound velocity can be considered as one limit- 
ing case, namely the case where the energy introduced 
by the explosion is small relative to the enthalpy of the 
gas involved. Taylor’s problem is the other extreme: 
he assumed that the energy introduced is very large in 
comparison with the enthalpy of the gas in which the ex- 
plosion cecurs. However, Taylor’s unsteady solution 
is three-dimensional in space; thus it cannot be used 
immediately for the description of a steady flow (except 
perhaps in a four-dimensional space). Recently S. C. 
Lin solved the problem of a violent cylindrical explosion 
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Fic. 8. 


Blunt body in hypersonic flow. 


and this computation led to the ‘‘blast-wave theory”’ of 
hypersonic motion. Especially, the variable distance be- 
tween the shock wave and the surface of a long body can 
be computed by this method, in that one builds up the 
flow picture in consecutive perpendicular planes of a 
sequence of solutions for the unsteady cylindrical ex- 
plosion at various states of the propagation process. The 
location of the shock wave of the steady flow prob- 
lem is identical with the location of the expanding wave 
at the corresponding time element. One assumes that 
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the explosion occurs at the cylindrical body surface and 
the energy introduced by the explosion corresponds to 
the energy introduced into the fluid by the resistance of 
the body; it is, in general, time-dependent. It is seen 
from Fig. 9 that the results of the blast-wave theory are 
in fair agreement with the observations. From the 
rather extensive literature on the subject, the contribu- 
tions of A. Sakurai, S. C. Lin, L. Lees and T. Kubota, 
and also H. K. Cheng and A. J. Pallone can be men- 
tioned. Independently, the theory was also developed in 
Russia, where the self-similar process was named an 
“auto-model.”’ In particular, L. Sedov and his collabo- 
rators worked on the theory of strong explosions. Grod- 
zovski, Chernyi, and Stanynkovich found and employed 
the analogy between the unsteady and the steady mo- 
tion. 

The true problem of the blunt body exposed to hy- 
personic flow involves aerothermochemistry, since the 
temperatures reached behind the shock are so high that 
the air dissociates. However, before discussing this 
problem, we want to consider some general questions 
related to boundary-layer theory. 


(VIII) BounpaRy-LAYER THEORY 


In 1954 we celebrated the fiftieth anniversary of the 
concept of the boundary layer, in view of the fact that 
Ludwig Prandtl presented the fundamentals of the 
boundary-layer theory in 1904 to the International 
Congress for Mathematics assembled in Heidelberg, 
Germany. For half a century the new concept proved 
to be one of the most fruitful ideas in fluid mechanics. 
The anniversary volume, entitled “Fifty Years of 
Boundary Layer Research,”’ can give an approximate 
picture of the main aspect of this development. We also 
want to mention that, in 1947, Loitsianski published a 
review of the contributions of Soviet scientists to bound- 
ary-layer theory. 

Most investigations of boundary-layer problems re- 
fer to plane or axisymmetric flows. We will mention 
later more general three-dimensional cases. At this 
point we only want to mention that Mangler succeeded, 
by means of a transformation, in reducing the axisym- 
metric case to a corresponding two-dimensional prob- 
lem. 

There are two different approaches which were used 
in the practical solution of boundary-layer problems, es- 
pecially in the case of incompressible fluids. One is 
the so-called integral method suggested by me and 
first used by K. Pohlhausen. This method reduces the 
problem of finding a solution of a partial differential 
equation to that of a solution of an ordinary differen- 
tial equation. Recently, an essential improvement of 
the method was achieved by I. Tani. The second ap- 
proach consists also in a reduction of the partial differ- 
ential equation to an ordinary differential equation by 
looking for special pressure distributions along the wall 
which allow similar solutions through the cross sections 
of the boundary layer. V. M. Falkner, S. W. Skan, D. 
R. Hartree, and B. Thwaites in the early thirties ex- 


celled in the development of such ‘similarity solu- 
tions.”’ 

As far as laminar compressible boundary-layer prob- 
lems are concerned, some special cases were solved be- 
fore 1946 by A. Busemann, myself, H. S. Tsien, W. 
Hantsche and H. Wendt, L. Crocco, and by H. W. Em- 
mons and J. G. Brainerd. Special attention was given 
to the heat-transfer aspect of the problems. 

The treatment of the compressible boundary layer in 
more general cases was greatly facilitated by a clever 
transformation of the independent variable across the 
boundary layer, which takes into account the variable 
density. Such a transformation apparently was first 
suggested by A. A. Dorodnitsyn and independently dis- 
covered by L. Howarth, K. Stewartson, and C. R. II- 
lingworth. 

Since this transformation essentially reduces the prob- 
lem of the compressible boundary layer to the incom- 
pressible case, the methods mentioned above could be 
applied to a broad field of problems. Thus, the integral 
method could be used for the flow about airfoils and 
bodies of revolution in supersonic flight. We can men- 
tion the contributions of L. E. Kalikhman, H. Weil, P. 
A. Libby and M. Morduchow, I. E. Beckwith, D. N. 
Morris and J. W. Smith, and others. 

Also, the method of similar solution was combined 
with the Stewartson-Illingworth transformation by L. 
Crocco and C. B. Cohen; they also extended the method 
of Thwaites by putting together sequence of similarity 
solutions. This procedure was further developed by 
C. B. Cohen and E. Reshotko. It appears that the 
work of S. Levy allows a great degree of generality in 
the formulation of the problem, viz., large tempera- 
ture changes, including viscous heating and arbitrary 
values of the Prandtl Number. This nondimensional 
quantity was taken in many previous investigations to 
be equal to unity. 

For a long time the theory of the laminar boundary 
layer was considered as having more academic than 
practical value. Recently, it was found that there are 
two reasons why the study of the laminar boundary 
layer is interesting also from the viewpoint of practical 
applications. First, there is the problem of re-entry 
of blunt bodies into the atmosphere. The highest 
value of heat transfer occurs at the stagnation point, and 
at the beginning the boundary layer probably has lami- 
nar structure. Second, the flight of missiles and ve- 
hicles at high altitude—i.e., in a medium of extremely 
low density—enhances the interest for conditions at 
small Reynolds Numbers. 

As we mentioned before, the problem of the boundary 
layer behind a detached shock produced by hypersonic 
motion of a blunt body is complicated by changes in the 
physical and chemical nature of the gas at the high 
temperatures produced by the shock. The main 
changes to be expected are changes in the distribution of 
the thermal energy over the degrees of freedom of 
the gas molecule, dissociation, and, finally, ionization. 
The first phenomenon leads especially to a variation in 
the value of the specific heat and particularly of the ra- 


tio 

fun 
by 


trar 
duc 
tica 
by 
Ina 
wav 
tect 
of th 
acte! 
ple 
assu 
hind 
cons: 
woul 
the t 
tion 
ble t 
none 
mate 
critic 
shock 
Th 
tion 
comb 
can b 
tion | 
natiot 
and } 
the er 
by L. 
rium 
condu 
make 
possib 
for th 
betwe 
by a si 
fer at « 
In g 
cernin; 
transfe 
Inve 
transfe 
cooling 
heat by 
called ‘ 
A re 
ments 
flat pla 
flow dit 
pears tl 
the usu 
layer is 
quite w 


apparer 


| 
| 
| 
Me 
| 


SIGNIFICANT DEVELOPMENTS IN AERODYNAMICS SINCE 


tio between the specific heats; the dissociation may 
fundamentally change the mechanism of heat transfer 
by introducing the possibility of diffusion. It is found 
that ionization does not essentially influence the heat 
transfer; however, the fact that the gas becomes con- 
ductive raises the question of whether there is a prac- 
tical way to change the structure of the boundary layer 
by artificially imposed electromagnetic-field effects. 
In addition, the behavior of an ionized gas against radar 
waves is of practical interest in view of problems of de- 
tection and communication in general. 

Concerning the methods for theoretical investigation 
of the influence of variable physical and chemical char- 
acteristics, two limiting cases appear as relatively sim- 
ple—-equilibrium state and frozen state. For example, 
assume a high degree of dissociation immediately be- 
hind the shock and consider this degree of dissociation 
constant over the flow field; recombination in this case 
would be restricted to the immediate neighborhood of 
the body surface. Then, of course, an important ques- 
tion of great practical value arises: whether it is possi- 
ble to prevent recombination at the surface by so-called 
noncatalytic surface coating. The application of such 
materials would essentially diminish heat transfer in the 
critical region and distribute the heat produced by the 
shock over a larger region. 

The more exact investigation of variable concentra- 
tion according to some rate law of dissociation and re- 
combination is more complicated, but the calculation 
can be carried out numerically at least for the stagna- 
tion point. The most extensive investigation of stag- 
nation-point conditions may be attributed to J. A. Fay 
and F. R. Riddell. The most complete discussion of 
the entire problem of the hypersonic boundary layer is 
by L. Lees. He found that two assumptions—equilib- 
rium state and Lewis Number (ratio between thermal 
conductivity and diffusion coefficient) equal to unity— 
make the computation relatively simple. Thus, the 
possibility arises of carrying out an exact investigation 
for the stagnation-point region and of using the ratio 
between the heat transfer at various points furnished 
by a simplified theory for an estimate of the heat trans- 
fer at other points. 

In general it seems that the various assumptions con- 
cerning the dissociation process do not change the heat- 
transfer values more than about 40 to 50 per cent. 

Investigations of boundary-layer flow with mass 
transfer are important in view of the possibilities of 
cooling by injection of material and also absorption of 
heat by melting or evaporation of the surface, generally 
called ‘‘ablation.”’ 

A remarkable phenomenon was found by experi- 
ments on boundary layers in hypersonic flow over a 
flat plate parallel to the flow direction or inclined to the 
flow direction at a small angle (sharp wedge). It ap- 


pears that in such cases the fundamental assumption of 
the usual theory—that the pressure in the boundary 
layer is equal to the pressure in the external flow—is 
quite wrong. One finds considerable pressure increase, 
apparently induced by the boundary layer itself. We 
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Illustration to Gértler’s theory of transition between 
laminar and turbulent boundary-layer flow. 
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mentioned that, probably due to the finite dimension of 
the leading edge, the flow picture corresponds to a de- 
tached shock caused by a kind of equivalent body corre- 
sponding to the domain occupied by the boundary layer. 
Then the motion inside this equivalent body has the 
general character of a flow through a pipe with increas- 
ing cross section, and this kind of motion produces the 
pressure increase. This phenomenon was called hy- 
personic boundary-layer - shock-wave interaction. 

There is no satisfactory criterion that would permit 
the prediction of the transition from laminar flow in the 
boundary layer to turbulent flow—for example, in the 
case of a blunt body. It is found both theoretically and 
experimentally that cooling of the surface helps to keep 
the flow laminar. This point has great practical im- 
portance because the heat transfer through a turbulent 
boundary layer is, in any case, much greater than the 
transfer through a laminar layer. 

Concerning the general problem of stability of the 
laminar boundary layer, the mathematical theory shows 
excellent agreement with experiment as far as the damp- 
ing or the increasing of oscillations is concerned, espe- 
cially after C. C. Lin corrected Schlichting’s original 
numerical calculations. However, this does not mean 
that we really understand the mechanism of transition. 
Recently, Emmons suggested the following mechanism: 
turbulence is created at isolated spots, and these spots 
cause contamination of the downstream flow. Ac- 
cording to Theodorsen the turbulence originates with 
the formation of horseshoe vortices which are formed 
near the walls and penetrate into the fluid. Pfenninger 
and Lachmann observed similar phenomena in their 
experiments to keep the boundary-layer flow laminar 
by means of suction. Géortler proposed the Schlichting- 
Tollmien waves, as a starting point for the understand- 
ing of the transition phenomena, but pointed out that 
the circulation around the curved streamlines plays a 
role in the transition, which according to him is an es- 
sentially three-dimensional phenomenon (Fig. 10). I. 
Tani and F. R. Hama have studied the influence of 
isolated roughnesses on the transition. 


Most investigations mentioned above refer to two- 
dimensional planes of axisymmetrical flow. In recent 
years, there has been considerable effort devoted to the 
interesting problems of the boundary-layer behavior 
when three velocity components are involved. Such 
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flows arise in a variety of ways; of particular interest 
for aeronautical applications are the flows over swept 
wings, over axisymmetric bodies at angle of attack, and 
in rotating blade systems. Recently, F. K. Moore pro- 
vided an excellent survey of developments in three-di- 
mensional boundary-layer theory. 

The boundary layer on swept wings influences the 
stalling and lateral control of the wing and has been the 
subject of considerable research. This problem is ideal- 
ized by considering the wing to be infinite, so that 
changes in the spanwise direction are neglected. One 
can distinguish between a boundary layer normal to the 
leading edge and a spanwise boundary layer. This 
consideration leads for incompressible flows to the im- 
portant “independence principle.” According to this, 
the two boundary layers are independent; one calcu- 
lates first the chordwise boundary layer in the usual 
way and then the spanwise boundary layer. The com- 
bination of the two boundary layers permits the stream- 
lines within the boundary layer to be constructed. 
In general, these streamlines differ from those in the 
external flow so that secondary flows arise. Research 
on this problem was carried out by L. Prandtl, V. 
Struminsky, R. T. Jones, and W. R. Sears. 

For the compressible case the density depends on 
both the chordwise and spanwise velocity components, 
so that the momentum equations in the two directions 
are “coupled” through the energy equation. The well- 
known Crocco integral applies in this case under cer- 
tain restrictions. However, the most general case of 
chordwise pressure gradient and heat transfer requires 
simultaneous solution of the chordwise and spanwise 
momentum equations. Recently E. Reshotko and 
I. E. Beckwith provided a solution for the boundary- 
layer and heat-transfer characteristics in the neighbor- 
hood of the stagnation line, such as arise at the leading 
edge of a blunt-nosed, swept wing. This analysis shows 
that significant reductions in heat transfer can be 
achieved by sweeping the leading edge of a wing; thus, 
the concept of sweep is seen to be useful also for hyper- 
sonic flight. 


(IX) THEORY OF TURBULENCE 


The theory of turbulence has two separate aspects. 
Evidently there is some analogy between the random 
motion of molecules in laminar flow and the random 
motion of eddylike formations in turbulent flow. The 
random motion of molecules leads to definite laws 
of molecular viscosity, heat conduction, and diffusion. 
In a similar way, the random motion we observe, on a 
much larger scale, in the turbulent flow of rivers, canals, 
pipes, and boundary layers apparently also results in 
definite laws for momentum transfer, energy transfer, 
and the transfer of matter, which we call turbulent 
friction, heat transfer, and diffusion, respectively. 
This analogy between the molecular and the turbulent 
processes was recognized and treated at a rather early 
date by Osborne Reynolds. The great difference be- 
tween the two concepts is the fact that in the case of 
the molecular random motion the elements are well de- 


fined, as molecules, whereas in the case of the turbulent 
motion they are not given a priori. Hence, whereas in 
the first case the Boltzmann theory shows the right way 
to a definite theoretical solution of the problem, in 
the case of turbulent motion new principles must be 
found. 

Because of this situation, the turbulence theory de- 
veloped in two essentially different directions. One 
school of thought endeavored to find half-empirical re- 
lations which would lead to definite rules for the predic- 
tion of the turbulent friction, heat transfer, and diffu- 
sion phenomena. On the other hand, very interesting 
ideas were proposed for building up a systematic sta- 
tistical theory of turbulence. 

As far as the first line of development is concerned, for 
the case of an incompressible fluid a satisfactory state 
of affairs was reached by the introduction of a concept of 
the mixing length—as a kind of generalization of the 
mean free path of the kinetic gas theory—due to L. 
Prandtl, and then, I believe, by the logarithmic law for 
the velocity distribution, which I found in 1930. As- 
suming the validity of the logarithmic law for the por- 
tion of the flow field which is mainly influenced by the 
wall, and the so-called ‘‘velocity decrement” law for the 
rest of the field, the turbulent flow in the boundary 
layer, in a circular pipe, and in a_ two-dimensional 
channel can be predicted for given Reynolds Numbers. 
The application of the integral method used by myself 
and Pohlhausen, also mentioned in connection with the 
laminar boundary layer, makes it also possible to calcu- 
late general boundary-layer problems. In this direc- 
tion, F. Clauser recently proposed interesting new con- 
cepts and computing methods. 

We have to mention that the concept of local similar- 
ity of the turbulent flow picture is an important factor 
as a lead for further development of the theory. When 
I first presented the idea in 1930, its practical formula- 
tion may have been somewhat over-simplified, but I 
believe—especially after the concept was clarified in 
some of my own publications together with C. C. Lin 
and in those of C. C. Lin and his collaborators—that it 
will have considerable influence on further development 
of the theory of turbulence. 

The main need today is for a reliable prediction of the 
practically important quantities, especially skin fric- 
tion and heat transfer for compressible fluids, particu- 
larly for supersonic boundary layers at higher Mach 
Numbers. Fig. 11 convincingly shows the discrepan- 
cies between the predictions derived from various pro- 
posed ‘‘theories.”’ 

The diagram, taken from NACA TN 3097 (1954), 
presents the ratio between the predicted value of the 
skin-friction coefficient in a compressible boundary layer 
and the well-known value of the same coefficient for in- 
compressible flow. The experimental value for IJ = 
4 is about 0.5. I notice that the value which goes with 
my name is closer to the experimental value than many 
others deduced by later elaborate theories. I con- 


sidered my so-called theory which I proposed more than 
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twenty years ago as a “‘guestimate.’’ I simply sub- 
stituted for the density, in the formula for incompres- 
sible flow, the value corresponding to the stagnation 
temperature which develops at the wall. This proce- 
cedure naturally underestimates the skin friction, since 
to be correct some appropriately determined average 
temperature should be used. Unfortunately, the re- 
sults of the theories of the Soviet scientists A. A. Do- 
rodnitsyn and L. A. Kalikhman are not included in 
the comparison. 

The systematic development of a statistical theory of 
turbulence was initiated by G. I. Taylor in 1935. He 
introduced the concept of isotropic turbulence, which 
was later modified and somewhat generalized by G. K. 
Batchelor to homogeneous turbulence. Batchelor’s 
book, published by the Cambridge University Press in 
1953, contains extremely valuable information on the 
developments of the statistical theory. Taylor pro- 
posed the average values of velocity correlations as a 
means of description of the turbulent flow. In 1938, L. 
Howarth and I gave a general analysis of the correla- 
tions and initiated the branch of the statistical theory 
which can be designated as the “dynamics of turbu- 
lence.’’ We arrived at a differential equation for the 
prediction of the variations of the correlation values 
with time—i.e., the propagation of turbulence. G. I. 
Taylor himself introduced the representation of the 
turbulent fluctuations by spectral functions as an al- 
ternative for the correlation functions. J. Kampé de 
Feriet investigated the concept of the spectral tensor 
especially from the mathematical point of view. The 
spectral viewpoint was adopted by Obukoff and Heisen- 
berg in their general investigations. Perhaps the most 
important result of these developments is a universal 
law for the spectral function—.e., turbulent energy vs. 
wave number—which is valid for high frequencies and 
was found independently by A. N. Kolmogoroff and W. 
Heisenberg. S. Chandrasekhar made significant con- 
tributions in this direction of the development. 

The theory of homogeneous turbulence—which 
covers the case of a turbulent fluctuation field without 
mean flow or the turbulent field superposed upon a 
uniform stream—is in relatively good shape. How- 
ever, the more interesting and practically more impor- 
tant case, the theory on the mechanism of turbulent 
shear, is still in its initial phase. From a physical point 
of view, the main difficulty seems to be that for the 
momentum transfer the large eddies are mostly respon- 
sible, whereas the energy dissipation is performed by 
the small eddies. This causes a diffusion of energy 
from low to high frequencies which does not comply 
with Onsager’s general energy-transfer scheme which 
provides the real basis for Kolmogoroff’s and Heisen- 
berg’s theories. 

From the mathematical point of view, one could ex- 
pect that the dynamic equation first written up by me, 
jointly with L. Howarth, would apply also to the general 
case of nonisotropic turbulence. However, this equa- 


tion essentially predicts the fate of the second-order and 
third-order correlations only if the higher correlations 
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are known. I believe real progress was made by the 
hypothesis that a so-called ‘‘quasi-normal joint-prob- 
ability’’ (Gaussian probability) can be assumed for the 
velocity field. This assumption makes it possible to 
express the fourth-order correlations essentially by 
products of the second-order correlations. It was first 
suggested by M. Millionschikoff. Attempts were re- 
cently made to apply the hypothesis for the solution of 
the shear problem by J. Proudman, W. H. Reid, T. 
Talsumi, and A. Craya. Other approaches for the 
solution of the shear problem were made by S. Chan- 
drasekhar, W. V. R. Malkus, J. M. Burgers, and M. 
Mitchner. All these theories are in an initial state. 
We are still far away from a theoretical deduction of 
the basic laws for shear and heat transfer which would 
give results in agreement with experimental evidence. 
I believe that the refinement of the similarity concept 
has to play a role in the establishment of a final theory. 

The experimental evidence on turbulent boundary 
layer and turbulent shear was very much enriched by 
the excellent work of many experimenters such as 
L. Kovasznay, J. Laufer, A. A. Townsend, G. B. Schu- 
bauer, and others. 

I mentioned before that ‘heavenly turbulence’’—.e., 
the turbulent phemonena observed in gaseous clouds of 
cosmic dimensions—fascinated the imagination of quite 
a number of mathematicians, astronomers, and aero- 
dynamicists. It was one of the problems which con- 
tributed to the interest in the most recent branch of 
fluid mechanics: ‘‘magnetofluidmechanics.”’ 

I may close with the remark that it appears to me 
that even in this so-called nuclear or space age, aerody- 
namics is not a science to be shelved as an obsolescent 
branch of the physical sciences. Of course, we may 
have to study the methods and results of many sister 
sciences more than in the past; viz., in addition to 
mathematics, mechanics, and thermodynamics, also 
chemistry and electromagnetism. However, I hope 
that the future training of young aerodynamicists will 
enable them to cope with the problems of the future. 

It is my agreeable duty to express my thanks and ap- 
preciation to my friends P. A. Libby and C. C. Lin for 
their assistance in collecting the material treated in 
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this lecture. Professor Libby collected the papers 
dealing with supersonic wing theory, slender aircraft, 
interference effects, higher approximations, hypersonic 
flow, and boundary-layer theory, and Professor Lin 
those dealing with the theory of turbulence. The 
choice of authors mentioned in the paper is somewhat 
arbitrary. I am also indebted to Prof. Wallace D. 
Hayes for letting me have the galley proofs of his book 
on hypersonic flow which he has written jointly with R. 
F. Probstein, and which will appear in the near future. 
Prof. I. Tani has prepared abstracts for me of impor- 
tant Japanese papers published in the last decade. The 
problems of flow in rarefied gases, sometimes referred to 
as problems of superaerodynamics are not treated here 
at all; the reader may find good summaries in the Pre- 
ceedings of an International Conference held in Nice, 
France, in July, 1958. Also, the aeroelastic problems— 
e.g., flutter at high speed—are only briefly mentioned. 
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Direct-Analog Method of Analysis of the 
Influence of Aerodynamic Heating on the 
Static Characteristics of Thin Wings 


Micnaet A. Basin,* RicHarp H. MacNeat,** AND H. Suiecps*** 


Computer Engineering Associates 


SUMMARY 


Application of the direct-analog computer method to the 
analysis of thin wings subjected to the influence of aerodynamic 
heating is described. A general method of analysis of temper- 
ature effects on structures is presented. The method is then 
applied to the example of a square, cantilever, constant-depth 
wing box, which is symmetrically heated along the leading and 
trailing edges. The electric analog computer is used to calculate 
the thermal in-plane stresses and the lateral load influence coeffi- 
Results are compared for the heated and unheated wing. 


cients 
SYMBOLS 

x,y, W = rectangular coordinates defined in Fig. 1 

U,, Uy, Uw = linear deflections, defined positive in the posi- 
tive coordinate direction 

0,, Oy = angular rotations defined in Fig. 1 

nN = dimensionless deflection parameter, Eq. (29) 

Ex, €y = axial strain 

Yzy = shear strain 

Ox, Sy = axial stress 

Try = shear stress 

OrT, SyT = equivalent axial stress produced by temper- 
ature distribution, Eq. (11) 

O27, QyT = thermally induced load density, Eqs. (12) and 
(13) 

a = linear coefficient of temperature expansion 

= temperature 

E = Young’s modulus, voltage 

G = shear modulus, current generator 

m = Poisson’s ratio 

D = flexural rigidity of a plate 

Nz, Ny = normal force per unit length 

Ney = shearing force per unit length 

q = externally applied lateral load density 

Pp = equivalent lateral load density, Eq. (15) 

P,, Ps = components of Pp, Eqs. (17) and (18) 

F.,, Fe, = total junction lateral loads corresponding to 
P, and P2, Eqs. (19) and (20) 

Ax, Ay = finite difference cell dimensions 

l = structure span and chord 

t = skin thickness 

h = height, the distance from top skin centerline to 
bottom skin centerline 

6 = dimensionless temperature parameter, Eq. (30) 

-L = symmetric or antisymmetric unit load 

F., F, = forces in axial load members 

R = resistance 

Ps, Fy = equivalent thermal loads in axial load members 

Sus Oy = horizontal panel shear loads 
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Vv = vertical shear force in vertical shear web 
P = primary of a transformer 
S = secondary of a transformer 


INTRODUCTION 


| geese THE PAST DECADE, direct analogy computer 
techniques have been applied to a wide spectrum 
of physical problems. In particular, such techniques 
have been extensively used in the areas of aircraft 
static, dynamic, and temperature distribution anal- 
ysis.'~* This paper describes a recent application of 
the direct-analog method to the static analysis of 
structures in which thermal stresses are significant. 

The direct-analog computer consists of a large num- 
ber of discrete lumped-parameter, passive, electrical 
elements (resistors, capacitors, inductors, transformers), 
a number of electronic devices such as operational 
amplifiers and arbitrary function generators, and con- 
trol circuitry. The basic underlying principle of oper- 
ation of this computing system is the direct behavioral 
correspondence between electrical circuit elements and 
discrete mechanical ‘‘building blocks’ such as levers, 
masses, and springs. The computer can be thought of 
as a large ‘tinker toy’’ which can be used to simulate 
numerous types of physical systems. 

There are essentially two basic methods for deriving 
electrical analogies for physical systems. The first 
method makes use of the correspondence between a 
set of simultaneous differential equations (or their 
finite difference equivalents) representing a physical 
system and a set of equations describing an electrical 
circuit. Individual coefficients in the system equa- 
tions are then identified with the coefficients of the 
electrical circuit equations. The second method, which 
is used for the analysis presented in this paper, makes 
use of the topological similarity of electrical and me- 
chanical systems. In this method, the continuous me- 
chanical system is first idealized into a prototype model 
consisting of discrete ‘building blocks.’’ Analogies 
are then derived for each component and are inter- 
connected on the computer in such a way as to satisfy 
the interconnection equations at all junctions and 
boundaries. 

This paper considers the application of the direct 
analogy method to the general thermal stress problem 
of thin wings. In particular, a square, cantilever, con- 
stant-depth wing box subjected to aerodynamic heat- 
ing is analyzed by means of the direct-analog com- 


‘ 
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puter, and results are presented and compared with 
those obtained for an unheated wing. 


EFFECT OF TEMPERATURE DISTRIBUTION ON 
STRUCTURES 


It is well known that the application of a temperature 
distribution to a structure produces several significant 
effects. First, stresses are generated in the structure 
by the thermally induced strains. Second, the re- 
sponse of the structure to externally applied loads is 
altered due to the presence of these induced stresses.® 
Numerous authors have analyzed the effect of thermal 
stresses on the torsional properties of beams. Vosteen 
and Fuller® present experimental results substantiating 
this phenomenon for a rectangular plate subjected to 
rapid heating conditions. The approach presented in 
this paper is to treat the temperature effects in light 
of the general equations of elasticity, and to arrive at a 
lumped prototype of the mechanical structure. The 
advantage of this method is its generality and its ap- 
plicability to the analysis of highly redundant and 
geometrically complex structures. 


The In-Plane Problem for Plates 


The equations representing the plane stress problem 
which include the effects of an externally applied tem- 
perature distribution are presented below. 


Strain Relationships: 

e, = OU,/Ox (1) 
= OU,/oy (2) 
= (0U,/Oy) + (OU,/Ox) (3) 

Stress-Strain Relationships: 
er = (1/E) (0, — woy) + aT (4) 
¢ = (1/E) (0, — woz) + aT (5) 
= (1/G)try (6) 

Equilibrium Relationships: 
(O0,/0x) + (O7,,/d¥) = 0 (7) 
(00,/Oy) + (O77,/0x) = 0 (8) 


The temperature effects can be interpreted in either 
of two ways. One way, as shown by Eqs. (4) and (5), 
is to add a strain, a7, to both e, and «,. The alternate 
approach is to rewrite Eqs. (4) and (5) as 

= (1/E) [(o, + + oyr) | (9) 
ey = (1/E) + oyr) — + (10) 
where Orr = oyr = alTE/(1 — pw) (11) 


The load densities per unit volume that must be 
applied at interior points of the ‘‘cold”’ structure to 
produce the thermally induced stresses are 


— (O¢z7/0x) (12) 
— (13) 


QrT 
QyT 


Normal stresses equal to a7 E/(1 — «) must be applied 
at the free edges of the structure. 

The suggested procedure for analyzing a plane stress 
problem including the effects of temperature is 

(1) Construct a lumped model or prototype which 
satisfies the above equations with temperature set 
equal to zero. Such a prototype will be composed of 
axial load carrying members (including Poisson’s ratio 
effect) and of idealized shear panels. 

(2) Subject the prototype to a loading distribution 
calculated to satisfy Eqs. (11), (12), and (13). 

(3) Measure or calculate the apparent internal 
stresses 0, + o;7, a, + oyr, and Tzy. 

(4) Subtract o,7 and o,7 from the apparent internal 
direct stresses to obtain the true direct stresses. 


The Lateral Load Problem 


It has already been noted that the in-plane thermally 
induced stresses alter the response of a plate-like struc- 
ture when it is subjected to lateral loads. This effect 
can best be studied by an examination of the plate 
equation,’ 


(01U,,/Ox*) + + (01U,,/dy*) = 
(1/D) [q + N,(0?U,,/0x?) + N,(0?U,,/Oy?) + 
(14) 


where the positive coordinate system for displacements 
and internal loads is illustrated in Fig. 1. N,, N,, and 
N,, are in-plane loads while q is the externally applied 
lateral load distribution. 

Examination of the right-hand side of Eq. (14) 
shows that the planar stresses in the skin effectively 
add ‘‘kick”’ load contributions to the applied distribu- 
tion g. The effective lateral load distribution can then 
be defined as 


P = q + + + 
2N,,(0?U,/Oxdy) (15) 


The load distribution of Eq. (15) can be divided into 
two parts. Making use of the definitions of 6, and 6, 
(see Fig. 1), 


P=P,+P, (16) 

where 
P, = q + N,(06,/dx) + N,(06,/dy) (17) 
= + (06,/dy) (18) 


The lateral problem can thus be solved by setting up 
a lumped prototype for the ‘“‘cold’’ mechanical system 
and then subjecting it to the load density prescribed 
by Eqs. (17) and (18). 

Many methods of analysis such as the one presented 
by Benscoter and MacNeal* make use of finite differ- 
ence techniques. Fig. 2 presents the points of defi- 
nition of the lateral deflection U,, and of the two 
orthogonal slopes 6, and 6, for a finite difference de- 
rived prototype. Since such a model can only be 
loaded at defined points, all lateral loads must be intro- 
duced at rib-spar intersections. The total junction 
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DIRECT-ANALOG METHOD OF ANALYSIS 


RIB NUMBER 


SPAR NUMBE 


P AND BOTTOM SKIN, EACH OF UNIFORM THICKNES +,” |S REMOVED IN ORDER TO SH 
2 S THE TOTAL HEIGHT FROM TOP SKIN CENTERLINE * BOTTOM SKIN CENTERLINE 
3. THE STRUCTURE 'S SQUARE WITH CHORE PAN = “Z." 


Fic. 1. Mechanical system. 


[RIB NUMBER 


OW WEB DETA 


SYMBOLS: Fw, FORCE AT RIB-SPAR INTERSECTION 
Fw, FORCE AT PANEL CENTER 


j+1+—SPAR NUMBER 


© DENOTES POINT OF DEFINITION OF ©, OR Oy 


IN ANALOG MODEL 


4 DENOTES POINT OF DEFINITION OF Uw 


IN ANALOG MODEL 


Fic. 2. Typical rib web-spar web intersection. 
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Fic. 3. Applied temperature distribution. 
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Fic. 4. In-plane problem applied thermal loads. 
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Fic. 5. Lumped model of a plane sheet. 
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(c) U, CIRCUIT (d) Uy CIRCUIT 


Fic. 6. Electrical analogy for axial load carrying elements, 
including Poisson’s ratio effect. 


load component corresponding to Eq. (17) can thus be 
written as 

Fy, = Pi Ax Ay = qi,; Ax Ay + 

j + (1/2) 

Nay. 


Jj 


— Oy, _ (1/2), (19) 


Yi+ (1/2), 9 


In order to calculate the load component described by 
Eq. (18), it should be noted that N,, is defined at panel 
centers and not at rib spar intersections. Consider 
the finite difference expansion of Eq. (18) about panel 
center [7 + (1/2), 7 + (1/2)]: 


Fy, = Ay = 4 (1/2), (ala) X 
[Ax (0, O25 54 am) + 
4 (1/2), (20) 


The equivalent load F,, can now be calculated at 
each panel center and then be prorated equally to the 
four adjacent rib-spar intersections. These prorated 
contributions can be added directly to the F,,, contribu- 
tions. 

Since the total equivalent lateral load distribution 
as given by Eqs. (19) and (20) is a function of measured 
slopes, the procedure for calculating or measuring the 
solution to the problem of a laterally loaded heated 
structure is as follows: 

(1) Apply load distribution, g, to the cold structure 
and measure appropriate slope differences for use in 
Eqs. (19) and (20). 

(2) Calculate the new value of the equivalent lateral 
load distribution by means of Eqs. (19) and (20). 


i+1, 7+ (1/2) 


AY (Oy, 541 


(3) Apply the new lateral load distribution to the 
structure. 

(4) Repeat the above iteration procedure until the 
solution converges. For problems where the effects of 
thermal stresses are reasonably small, one or two iter- 
ations should be satisfactory. The iterative process 
will fail to converge if the structure buckles under the 
action of thermal stresses. 


AN ILLUSTRATIVE EXAMPLE 


In order to demonstrate the above techniques, a 
simple problem will be formulated and analyzed by 
means of the direct analogy analog computer. The 
system to be studied (Fig. 1) is a square, cantilever, 
constant-depth wing box covered with a constant 
thickness skin on both top and bottom sides. The 
wing box contains four equally spaced rib webs inter- 
sected orthogonally by five equally spaced spar webs. 
The temperature distribution for the structure, which 
is assumed to be the same at all spanwise stations, is 
shown in Fig. 3. The equivalent external loads to be 
applied to the in-plane structure, which are produced 
by the average temperature distribution of Fig. 3, are 
calculated with the aid of Eqs. (11), (12), and (13) and 
are presented in Fig. 4. The thermal stresses induced 
by this temperature distribution will be calculated 
first. The defiection influence coefficients for lateral 
loads will then be calculated, first in the absence of 
thermal stresses (“‘cold’’ wing), and then for two cases 
corresponding to different temperature levels (‘‘hot”’ 
wing cases). 


The Construction of the Lumped Prototype 


In order to analyze a distributed mechanical system, 
a prototype model composed of lumped ‘building 
blocks’’ must first be constructed. In the construction 
of this finite difference model, numerous approxima- 
tions and assumptions are made. For the purposes 
of the sample problem presented in this paper, simpli- 
fying assumptions are used which may not be generally 
valid, but their use should not be considered as a re- 
striction on the general method. The following as- 
sumptions and representations are embodied in the 
lumped model: 

(1) The structure is assumed to possess symmetry 
about the horizontal wing center plane and about the 
vertical plane passing through spar (00) (see Fig. 1). 

(2) The rib and spar webs are assumed to remain 
unstressed when the temperature distribution is applied 
to the wing. This is a valid assumption if the webs 
are assumed to be constructed in such a manner as 
to relieve thermal stresses completely. 

(3) Although the rib and spar webs are allowed to 
carry vertical shear, the strain energy in these webs is 
neglected. 

(4) All axial loads are carried by the skin—i.e., the 
spar and rib cap areas are assumed negligibly small in 
comparison to the equivalent cover skin axial load 
carrying areas. 
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(5) The state of stress in the skin is considered 
separable. The cover skin is thus represented by ideal- 
ized axial load carrying members and by idealized skin 
shear panels which are in equilibrium under shear 
loads. 

(6) The Poisson’s ratio effect between normal stresses 
in the cover skins is included in the analysis. 

(7) In order to emphasize the effects of thermally in- 
duced stresses on the static behavior of the wing sub- 
jected to lateral loads, the values of the elastic con- 
stants & and G are assumed to be invariant with respect 
to temperature. 


Symmetry Considerations 


In constructing prototypes of mechanical systems 
such as the one under consideration, full use may be 
made of structural and load symmetry. In addition 
to a horizontal plane of symmetry which is used to 
separate the in-plane and lateral problems, the struc- 
ture of Fig. 1 has a vertical plane of symmetry about 
spar (00). Making use of this fact, only half of the 
prototype need be constructed. Then, by proper 
application of boundary conditions at this plane of sym- 
metry, the response of the system to symmetrical and 
antisymmetrical loading conditions can be obtained. 
All that is required is to decompose the applied load 
distribution into its symmetric and antisymmetric 
components. 

It should be noted that, for the in-plane problem, 
the assumed thermal loading distribution is symmet- 
rical with respect to the center spar; hence, only 
one set of boundary conditions along spar (00) is 
required for the in-plane problem. The loading in the 
lateral load problem is, however, unsymmetrical. 
The lateral problem is analyzed by applying unit sym- 
metric and antisymmetric lateral loads to the struc- 
ture. Results are presented in terms of symmetric 
and antisymmetric influence coefficient matrices. 


DERIVATION OF THE ELECTRICAL CIRCUITS 


A direct-analog computer is essentially a large 
“tinker toy’ containing electrical elements which can 
be made to possess behavioral correspondence to levers, 
masses, springs, dashpots, and forcing functions. 
Electrical measuring instruments are used to record 
the response of the electrical (analog) system. Since 
the problem being analyzed is a static one, the analog 
model can be constructed from levers, springs, and 
forcing functions only. The analogous quantities 
given in Table 1 are used. 

Seale factors are normally used in the relationships 


TABLE 1 


Mechanical Quantity Electrical Quantity 


Displacement, U Voltage 
Rotation or Slope, @ Voltage 
Force, F Current 
Moment, Current 
Spring Constant, K (Resistance) ~! 


Leverage, and other similar effects Transformer 


between electrical quantities physically present in a 
computer and their analogs. The selection of scale 
factors is an important practical consideration that 
will not be discussed here because it is not essential. 
The reader is referred to reference 1 for details on scale 
factors. 

The method that will be used for deriving the elec- 
trical analogy starts with the observation that the 
lumped mechanical prototype consists of a (large) 
number of simple structural elements which can be 
grouped into two or three basic types. The basic 
types of elements considered here are axial load carrying 
members, horizontal shear panels, and vertical shear 
webs. The derivation of the electrical analogy for the 
complete structure is accomplished in two steps. In 
the first step, electrical analogies are derived for the 
basic types of simple structural elements. In the 
second step, the electrical analogies for the elements 
are connected together so that the equilibrium relation- 
ships and displacement compatibility relationships 
existing at the intersections of the structural elements 
are satisfied. 


Derivation of Electrical Circuit for the In-Plane Problem 


Consider the lumped model of a plane sheet shown 
in Fig. 5. In this model each shear panel is connected 
to adjacent axial load members at one point on each 
side of the panel. Hence the pair of axial load members 
that cross at the intersection of panels A, B, C, and D 
are unloaded except at their ends and may, conse- 
guently, be removed from the structure and analyzed. 
These members and their electrical analogy are shown 
in Fig. 6. The equations satisfied by the members, as 
derived from Eqs. (9) and (10) by making finite dif- 
ference approximations, are 


— t Ay ax Gt dy ~Y 
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Fic. 7 Electrical analogy for a shear panel. 
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Fic. 8. Circuit for in-plane problem. 


Use = [(Ax/Ay) (F, + 

u(F, + F,’)] (21) 
Uya — = (1/Et) [(Ay/Ax) (F, + F,’) 

+ (22) 
where 


F, = a,(tAy), F, = a,(tAx) 


F,’ = aTE(tdy)/(1 — 
F,’ = aTE(tAx)/(1 — 


In the analogous electrical circuit, Fig. 6(c) and (d), 
currents are represented by force symbols and voltages 
are represented by displacement symbols. Voltages 
at nodes are the electric potentials with respect to a 
common reference (ground). Eqs. (21) and (22) can be 
interpreted as the expressions for the voltage drops in the 
U, and U, circuits, respectively. The current flowing 
in R, is F, + F,’ and the current flowing in R, is FP, + 
F,’. Because of the presence of the unity turns ratio 
transformer, the current flowing in R, is (F, + F,’) — 
(F, + F,’). Summing the voltage drops across R, 
and R,, one obtains 


Use RAF, + + R,I(F; + F,’) 
(F, + F,’)] (23) 


If R, is made equal to [(Av/Ay) — u]|/Ft and R, is 
made equal to u/Et, Eq. (21) is obtained. The satis- 
faction of Eq. (22) by the circuit can be demonstrated 
in similar fashion. 

Shear panel A in Fig. 5 and its electrical analogy are 
shown in Fig. 7. The equations satisfied by the forces 
and displacements at the edges of the panel, as derived 
from Eqs. (3) and (6) by making finite difference 
approximations, are 


S,Ay = S,Ax (24) 


[(Use — Use)/Ay] + — Uya)/ Ax] = 
(1/GtAy)S, (25) 
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= Try(tAx), Sy = t2,(tAy) 


These equations are satisfied by the analogous elec- 
trical circuit shown in Fig. 7(c) and (d) where, as in 
the case of the axial load members, force symbols are 
used for currents and displacement symbols are used 
for voltage. The turns ratio of the ideal transformer 
coupling the U, and U, circuits is equal to Aw/Ay. 
Eq. (24) is, consequently, a statement of the relation- 
ship between the currents in the two windings of this 
transformer. The voltage across resistor R, is 


Ep, = Uys — Uy + (Ax/Ay) (Use — Usa) = 
S/R, (26) 


where 


If R, is made equal to Ax/GtAy, Eq. (26) becomes iden- 
tical to Eq. (25). 

With the electrical circuits for the axial load mem- 
bers and the shear panel derived above, the stage is 
set for the second major step, which is to interconnect 
a number of these circuits to form an analogy for the 
complete structure. It will be noted in Fig. 5 that, 
at point (a), two shear panels and two axial load mem- 
bers are joined together. The displacement com- 
patibility condition is siinply that U’,, is the common 
x displacement for these four members. The equilib- 
rium condition at point (a) is that the sum of the forces 
exerted on point (a) by the four adjacent members is 
zero. This condition is derivable from the displace- 
ment compatibility condition and the statement that 
the forces acting on point (a) do no work. In the 
electrical analogy, the displacement compatibility con- 
dition is satisfied by joining the terminals of the 
four circuits together so that their voltages are equal. 
Satisfaction of the no-work condition and, hence, of 
the equilibrium condition is assured since a solid con- 
nection consumes no power (the analogy of strain 
energy). In general, it is necessary to demonstrate 
either the displacement compatibility condition or the 
equilibrium condition, but not both, in deriving the 
electrical analogy of a built-up structure. 

In a similar manner, the U, circuits of two shear 


< 
—AXx — 
Fxg Fygh + Vgdx = O 
h Ura Uwe 
“| Ax h 
(a) VERTICAL SHEAR WEB (b) EQUATIONS 
} 
V, 
Va Ox 
Uws Uwe 
(d) ELECTRICAL ANALOGY 
| FOR 1/2 OF SHEAR WEB 
| ASSUMING ANTISYM- 
| METRY WITH RESPECT 
(c) ELECTRICAL ANALOGY TO NEUTRAL PLANE 


Fic. 9. Electrical analogy for a vertical shear web. 
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00 NOTES: 
1.FOR IN-PLANE PROBLEM, OPEN ALL SWITCHES (T), AND 
CONNECT ALL SYMMETRY SWITCHES TO S TERMINALS; 
THE CIRCUIT THEN REPRESENTS U, AND Uy. THE Uy 
CIRCUIT IS NOT USED 


2.FOR THE LATERAL LOAD ANALYSIS, CLOSE ALL SWITCHES 
(T). CONNECT SYMMETRY SWITCHES ACCORDING TO TYPE 
OF LOAD. (SYMMETRY IS DEFINED WITH RESPECT TO 
THE VERTICAL PLANE PASSING THROUGH SPAR 00.) 


3.—O— SYMBOL REPRESENTS METERING POINT FOR 
VOLTAGES NECCESSARY FOR ITERATION PROCEDURE 


4. @— SYMBOL REPRESENTS CURRENT GENERATORS 
TO BE SET FOR IN-PLANE PROBLEM 


5. ®— SYMBOL REPRESENTS CURRENT GENERATORS 
TO BE SET FOR LATERAL PROBLEM 
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Fic. 10. Circuit diagram. 
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Fic. 11. In-plane problem computer data. 


panels and of two axial load members are joined to- 
gether at point (b) in Fig. 5. 

The complete electrical analogy for the in-plane 
thermal stress problem is shown in Fig. 8. 

Note that although there is no direct electrical con- 
nection between the U’, and U, circuits, these circuits 
are magnetically coupled by transformers. 

Boundary conditions are treated in a natural manner. 
Along a clamped edge, voltages are constrained to be 
zero. A free edge is characterized by an absence of 
electrical constraints. At the plane of symmetry along 
, Spar (00), voltages corresponding to x displacements 
* are constrained to be zero and voltages corresponding 
to y displacements are unconstrained. 

The circuits for axial force members can be simplified 
at points along the free edges because the force normal 
to the edge is zero. For example, along spar (20f) 
in Fig. 8, F, is zero, so that Eq. (22) becomes 


Uya — Up = (1/Et) {(Ay/Ax) F, + 
Uya — Up = (1/Et) (Ay/Ax) X 


[F, + (1 u) 


Hence, an axial force member along spar (20f) is 
represented by a simple resistor with a constant current 
generator at each end. 

It may be seen from Fig. 6(c) and (d) that two con- 
stant current generators occur at the junction of axial 
force members, but since one generator inserts current 
into the node and the other withdraws current, the net 
effect can be represented by a single generator insert- 
ing a current proportional to the difference in F,’ (or 
F,’) at adjacent rib-spar intersections. Furthermore, 
since the temperature and the physical properties of 
the structure are uniform in the y direction, generators 
for the U, circuit are required at the outboard end 
only. 

A simplification of the circuit also occurs at the 
clamped edge since it is assumed that the strain parallel 
to the clamped edge is zero. Hence, along rib (00) in 
Fig. 8, it is found that, by setting U,, — U,, equal to 
zero in Eq. (21), 


Up = [(1 — u°)/Et] (Ay/Ax) (F, + F,') (28) 


Uva 


An axial force member normal to the clamped edge 
is, consequently, represented by a simple resistor whose 
value is (1 — yp?)Ay/EtAx. Ay is, in this expression, 
equal to one half of the rib spacing. 


Derivation of Electrical Circuits for the Lateral Load 

Problem 

An electrical analogy for the bending of flat plates 
is described in reference 9; it is shown in reference 8 that 
the electrical analogy for the bending of thin multicell 
shells with orthogonal ribs and spars has an identical 
form. Hence the circuits described in references 8 
and 9 are directly applicable to the lateral load problem. 

The electrical analogy for the lateral load problem 
can also be obtained by an extension of the analogy used 
in the in-plane problem. In the lateral load problem 
the cover sheets are still in a condition of plane stress, 
but they are fastened to an internal system of vertical 
shear webs. In the lumped model it is assumed that 
the vertical shear webs are connected to the cover 
sheets at points midway between ribs and_ spars 
[points (a), (b), (c), and (d) in Fig. 5], and that they 
are connected to each other at rib-spar intersections. 
The electrical analogy for a vertical shear web in which 
shear strain is neglected is shown in Fig. 9. Since 
the structure is symmetrical with respect to a plane 
midway between the cover sheets and the loading is 
antisyimimetrical with respect to the same plane, U,,’ 
is equal to — U,, in Fig. 9(a). This condition implies 
that the voltage at the center point of the vertical 
transiormer winding in Fig. 9(c) is zero, and then, 
consequently, an analogy that represents only the 
lower one-half of the shear web can be formed, as shown 
in Fig. 9(d). 

The complete analogy for the lateral load structure 
is constructed by adding electrical elements represent- 
ing vertical shear webs to the circuit previously derived 
for the in-plane problem. The complete circuit shown 
in Fig. 10 is used for both the lateral and the in-plane 
load problems, with the transformers that represent 
vertical shear webs being removed for the in-plane 
problem. The constant current generators appearing at 
rib-spar intersections in the U,, circuit simulate one- 


1000 2000 3000 4000 3010s 4010¢ 1020 2020¢ 3020; 4020; 


1000 ].045 .072 .100 -03! 064 .096 128 .02 .054 .090 123 
2000/.072 .206 .322 442 .064 .185 -309 434 .164 .290 418 
3000] .100 .322 .602 90 .094 306 .576 872 279 832 
4000} .13) 442 890 i400 -428 668 1.340 388 -280 
1010¢ | .O3! 064 .094 125.035 .064 .094 124 .059 .092 123 
2010¢] .064 185 -306 -428 .064 .308 426 .299 -420 
30105] .096 309 576 .094 .308 .578 864 .082 .290 
40104] .128 434 872 1.340 .124 426 .864 .105 399 832 -300 
1020¢| .02! 050 .078 10! .028 6.056 .082 105 .050 073 .093 114 
20204] .054 .279 -388 .059 -290 210 318 42! 


30207] .090 .290 542 812 .092 .299 558 832 .093 .318 602 872 


40204 123 418 832 1.280 .123 420 .848 1.300 «114 42 


Fic. 12. “Cold” structure symmetric deflection influence 
coefficients (all data presented in units of the dimensionless 
parameter )). 
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half of the lateral loads given by Eqs. (19) and (20). 
The increments in slopes appearing in these equations 
are obtained from measurements of voltage increments 
in the cover sheet circuits. The slopes are, of course, 
equal to the displacements in the lower cover sheet 
divided by the semidepth of the wing. 


Identification of Elements in the Complete Circuit 

Diagram 

The electrical circuit diagram which represents both 
the in-plane and the lateral problems is presented in 
Fig. 10. 

The functions of the elements in Fig. 10 are as 
follows: 


PandS 1to20: Vertical web (lateral problem 
only) 

PandS 21 to26: Poisson ratio (both problems) 

PandS 27to34: Horizontal shear panel (both 
problems) 

R1 to R6: R, for regions with Poisson ratio 
(both problems) 

R7 to R12: R, for regions with Poisson ratio 
(both problems) 

R13 to RIS: R, for regions with Poisson ratio 
(both problems) 


R19 to R24: R, for regions without Poisson 
ratio (both problems) 

R, for regions without Poisson 
ratio (both problems) 
R27 to R34: R, Horizontal shear panel 
flexibility (both problems) 
Metering dummies for purposes 
of measuring slope differences 
needed to calculate iteration 
currents in lateral load prob- 
lem 

Switches to disconnect web 
transformers in order to con- 
vert circuit to in-plane prob- 
lem 

Boundary condition symmetry 
switches 

Current generators used to in- 
troduce lateral loads 

Current generators used to in- 
troduce effects of thermal 
strain for the in-plane problem 


R25 to R26: 


M1 to M17: 


to Too: 


Si to Si! 
Gl to Gy: 


G17 to G23: 


THE COMPUTER SOLUTION 
The In-Plane Problem 


The self-equilibrating internal loads calculated by 
making appropriate measurements on the circuit of 
Fig. 10 are presented in dimensionless form in Fig. 11. 
It can be noted that the regions in the vicinity of the 
leading and trailing edges are in compression while the 
central region of the wing is in tension as expected. 
These data are then used in conjunction with Eqs. 
(19) and (20) and the lateral problem circuit to calculate 
the lateral influence coefficients. 


The Lateral Problem 


By proper switch connections, the circuit of Fig. 10 
can now be used to calculate the lateral deflection in- 
fluence coefficients for both the “hot” and ‘‘cold”’ 
structures. Data are presented in the form of deflec- 
tion influence coefficient matrices for symmetric and 
antisymmetric unit loads. The results are presented 
in terms of a dimensionless deflection parameter 4, 
where 


1000 2000 3000 4000 2010¢ 30104 4010, O20, 20204 30204 4020% 


1010¢ 0180 0169 0189 0192 0200 .0302 0345 .0362 
2010" | O C 0189 O46 0510 0544 0306 ogi2 000 
30104 | O 0189 0510 1000 0352 o9i9 1450 750 
4010¢] 0 0 0192 0544 000 590 0374 1020 1780 2520 
10204] 0 0 0200 0352 0374 05:9 0655 O76 0754 
20204] 0 0 .0302 o9i9 020 0655 580 880 2050 
30204] ie) 0345 450 780 O716 880 3040 3520 
40204] O 0362 000 750 2520 .0754 2050 3520 5090 


Fic. 13. “Cold” structure antisymmetric deflection influence 
coefficients (all data presented in units of the dimensionless 
parameter )). 


000 2000 3000 4000 2010% 4010, O20¢ 2020¢ 3020¢ 4020, 


1000 0 
| 0007 .00!5 .0020 .0024 .0016 .003i .0041 .0046 
2010¢ | .0037 .0054 .0066 .0028 .0075 .0107 0125 
30104 | 0020 .0054 .0084 .0109 .0035 .010 0160 .0200 
40104] © 0024 .0066 .0/09 50 .0038 .0N6 0260 
10204] © O O01 0028 .0035 .0038 .0036 .0058 .0070 .0076 
20204] 003! .0075 .O10 .0058 .0150 .02 0229 
30204] 0 004! .0107 .0160 .0198 .0070 .0201 .0314 .0376 
40204] 0 fe 0 0046 .0125 .0200 .0260 .0076 .0229 .0376 .048) 
92 


Fic. 14. Antisymmetric difference deflection influence co- 
efficient matrix, 8 = 3.192 (all data presented in units of the 
dimensionless paramete. 


1000 2000 3000 4000 1010¢ 2010¢ 40104 O20¢ 20204 30204 4020¢ 


1010¢ | O ° 0032 .0044 .0052 .0034 .0068 .0089 .0102 
20104 | O 0032 .00860 0142 006: .0162 .0234 0279 
3010, | O .0118 0186 0238 .0076 .022! 0356 0445 
40104} ie} 0052 .0142 .0238 .0325 .0084 .0256 .0432 .0574 


10204] O 0034 .006' .0076 .0084 .0078 .0128 o1s2 o168 


20204] ° 0068 .0162 022 0256 0126 .0334 .0442 .0508 
30204] O 0089 .0234 .0356 .0432 .0442 .0689 .0838 
40204] O 0102 0279 .0445 .0574 .0508 .0838 .1070 
pe 384 


Fic. 15. Antisymmetric difference deflection influence co- 


efficient matrix, 8 = 6.384 (all data presented in units of the 
dimensionless parameter )). 
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A = U,th?E/PL (29) 


of the structure. For symmetric boundary conditions, 
the total load applied to the complete structure is 2L. 
For antisymmetric boundary conditions, the total load 
is an up load, L, at the point of application (forward 
part of structure) and a down load, ZL, at the corre- 
sponding point on the aft part of the structure. 

Data are presented in Figs. 12, 13, 14, and 15 for 
three values of the dimensionless temperature param- 
eter 8, where 


B = (l/h)? aT, (30) 


In Figs. 14 and 15, for which 8 is not equal to zero, 
only the antisymmetric difference deflection influence 
coefficient matrices are presented. The difference ma- 
trix is defined as the difference between the ‘‘hot’’ in- 
fluence coefficients and the “cold” influence coeffi- 
cients. Symmetric difference influence coefficients were 
measured but are not presented since the applied tem- 
perature distribution does not significantly effect the 
symmetric results. The maximum effects in sym- 
metric cases are approximately one tenth of the maxi- 
mum effects in antisymmetric cases. 

It may be worthy of mention that the results of the 
analysis are consistent with conclusions based on ele- 
mentary beam theory that (1) self-equilibrating stresses 


and ZL is defined as the total load applied to one half 
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have no (or little) influence on the effective bending 
stiffness, and (2) self-equilibrating stresses of the type 
produced by aerodynamic heating produce a decrease 
in the effective torsional rigidity of the structure. 
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Investigation of Normal Force Distributions 
and Wake Vortex Characteristics of Bodies 
of Revolution at Supersonic Speeds' 


J. F. MELLO* 


McDonnell Aircraft Corporation 


SUMMARY 


The supersonic aerodynamic characteristics of inclined bodies 
of revolution at high angles of attack have been investigated in 
order to provide a more basic understanding of the body vortex 
wake flow and its relation to the problem of body-wing inter- 
ference. The results of wind-tunnel tests, whereby the normal 
force, pitching moment, normal force distributions, and the local 
flow properties in the vicinity of the body were determined, are 
discussed and analyzed. 

Comparisons of experimental normal force coefficient and 
center of pressure data with values calculated in accordance with 
theories which include methods for estimating the effects of vis- 
cosity show that the accuracy of these estimates is strongly de- 
pendent on the body fineness ratio and the angle of attack. 
Further comparisons of the distributions of theoretical and ex- 
perimentally derived cross-flow drag coefficients clearly show that, 
in general, the disagreement between experiment and existing 
theories is due to the inadequate prediction of the magnitude and 
distribution of the forces resulting from flow separation 

The circulation strengths of the concentrated vortices and the 
circulation strengths of the vortex feeding shects in the body 
vortex wake are determined by closed-contour velocity-perimeter 
integrations for paths enclosing the vortex or the feeding sheet. 
The values of vortex strength calculated in this manner are in 
close agreement with the values predicted by vortex strength 
formulas written for a simple theoretical model for which it 1s 
assumed that the cross-flow in any plane along the cylindrical 
portion of the body is represented by the steady incompressible 
potential flow about a cylinder, two symmetrical vortices of 
equal strength, and the attendant image vortices. However, 
in computing these strengths it is necessary to use the vortex 
locations and the viscous normal force distributions determined 
from experiment. 

The experimentally determined values of vortex strength are, 
in turn, used to calculate—by means of the aforementioned in- 
compressible cross-flow potential—the local flow inclination 
angles which are in good agreement with the measured values, 
except in the vortex core, in the vicinity of the feeding sheet, 
and in regions for which transonic cross-flow velocities are ex- 
pected. A consideration of these various regions with simple 
methods which account for the observed phenomena leads to 
substantial improvement in the agreement between theory and 
experiment. 

It is indicated that the complete vortex wake flow niay be ade- 
quately predicted for a body of revolution (for conditions repre- 
sented by the theoretical flow model), provided that the dis- 
tribution of the viscous normal force and the vortex locations are 
accurately known. 
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SYMBOLS 


body radius of the cylindrical section 

local cross-flow drag coefficient based on 
diameter 

cross-flow drag coefficient of an infinite circu- 
lar cylinder, per unit length, for steady- 
state flow, based on diameter 

local normal force coefficient per unit length 


— (2a/S) f C,, sin 6d0 
total normal force coefficient, 
l 
C, dz = N/q S 


pressure coefficient, p — po/dao 

body diameter, cylindrical section 

imaginary number, 

overall body length 

local Mach Number 

free-stream Mach Number 

cross-flow Mach Number 

local normal force per unit length, nm = dN /dz 

resultant normal force 

local static pressure 

free-stream static pressure 

local total pressure, corrected for shock losses 
and losses due to local flow inclination 

free-stream total pressure 

free-stream dynamic pressure, pp 

complex notation of a point in the cross-flow 
plane, s = x + = re”® 

reference area, 7a? 

time 

free-stream velocity 

total local velocity, Ui? = (Up cos ag + w)? + 
u? + y? 

cross-flow velocity components in x and y 
directions. respectively 

radial and tangential velocities in the cross- 
flow plane, respectively 

total local cross-flow velocity, V ut += 
V 0,2 + 062 

free-stream cross-flow velocity, Uo sin ag 

axial perturbation velocity 

model Cartesian coordinates, origin at ver- 
tex, z positive downstream along body axis 


model cylindrical coordinates, 6 = 0° in 
horizontal plane of model on right-hand 
side, looking upstream, while @ = 90° 


denotes the top location 

angle of attack of body 

tangential component of the local flow in- 
clination angle, tan a = vg/U) cos ag 

radial component of the local flow incli- 
nation angle, tan 8 = v,/U, cos ag 

total local flow inclination angle, tan « = 
cos ag 

nondimensionalized circulation strength 
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EXPERIMENT 
THEORY 

REF. 2 

- - REF 2+ KELLY 4 “he =1,20 
=---5---- - REF 2 + KELLY4 =.35 


° Mo=2.00 

Cn | 


O 8 16 24 
O p DEG. 
Fic. l(a). The normal force and center of pressure char- 
acteristics of a 30° cone-cylinder configuration (//d = 8.03). 
o = position of vortex center in the cross-flow 
plane, ~ + in (right-hand side) 
—t = position of symmetrically located vortex in 


the cross-flow plane, —& + in 
—a’/o, a?/é = image vortex location, + in; and + 
ini, respectively 


& on = Cartesian coordinates of right-hand real 
vortex 

gi, on = Cartesian coordinates of right-hand image 
vortex 

Po = free-stream air density 

W = compiex potential, + 

e = potential function in the cross-flow plane 

Vv = stream function in the cross-flow plane 

¥ = ratio of specific heats (y = 1.40 for air) 

In = natural logarithm 

RI = real part of complex function 

Subscripts 


p = potential contribution 

viscous contribution 

image 

conditions at axial station z from nose 


v 


(1) INTRODUCTION 


eos INTERACTION of the vortex wake, which has 
been observed on the leeward side of bodies of 
revolution at moderate angles of attack in a supersonic 
stream, with any attached lifting surfaces is an im- 
portant aspect of the total problem of wing-body inter- 
ference. It has been shown that this vortex wake at 
high angles of attack can have a large effect upon the 


forces and the distribution of forces developed on wings 
or tails immersed in this flow. Since methods! are 
available for the calculation of lifting surface loads and 
load distributions for nonuniform supersonic flows, the 
immediate problem is that of predicting the properties 
of the flow in the vicinity of the lifting surface. 

This paper presents the results of an experimental 
investigation of the aerodynamic characteristics of 
an inclined body of revolution and assesses the ade- 
quacy of available methods for estimating these char- 
acteristics. A 30° cone-cylinder configuration has been 
used to obtain normal force and pitching-moment data, 
normal force distributions, and local flow properties 
for these purposes. Comparisons of experimental and 
theoretical results are presented for a wide variety of 
conditions for which experimental Gata were obtained 


(2) ForcE AND MOMENT CHARACTERISTICS 


The normal force and pitching-moment character- 
istics of inclined bodies of revolution at supersonic 
speeds have been treated by several authors. Theories 
are available for the prediction of these characteristics 
in both the low angle of attack and high angle of attack 
regimes. The high angle of attack theories are semi- 
empirical in nature and consist of the addition of a 
viscous contribution to a potential solution which is 
usually taken to be linear with the angle of attack. 

In the present studies, the high angle of attack char- 
acteristics were of more immediate concern. There- 
fore, the potential solutions for the configurations dis- 
cussed herein were calculated using the hybrid theory 
of Van Dyke,? and comparisons with other theories 
were not made. The viscous contributions to the 


12 4 | 


Fic. 1(b). The normal force and center of pressure characteris- 
tics of a 30° cone-cylinder configuration (//d = 13.14). 
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body normal force magnitudes and distributions in all 
of the cases discussed were calculated by the method 
suggested by Allen,* and, for comparison purposes, by 
Kelly's method‘ as well. 

The variation normal force coefficient and center of 
pressure with angle of attack at Mach Number 2 for 
two of several 30° cone-cylinder bodies of varying fine- 
ness ratio which were tested are presented in Fig. 1. 
Also shown are the predicted values using the hybrid 
theory of Van Dyke to which has been added the 
viscous contribution obtained by the methods of Allen 
and Kelly. The results shown for Kelly’s method with 
c’a, = 0.35 and ca, = 1.20 are for the case of a turbu- 
lent boundary layer and for the case of a laminar 
boundary layer, respectively. 

These data illustrate how the methods of Allen or 
Kelly closely estimate the measured coefficients in cer- 
tain cases but, nevertheless, indicate that neither 
method is satisfactory in all cases. In particular, it is 
noted that the accuracy of both methods depends on 
the fineness ratio of the body and on the angle of 
attack. 

Because of the generally poorer agreement (except 
at very low angles of attack) shown for Kelly's method 
with Ca, = 0.35, subsequent calculations in this report 
are made with C'a, = 1.20 only. Furthermore, suffi- 
cient evidence is available to indicate that the form 
drag characteristics of a laminar-flow type separation 
are preserved (superficially, at least) for supersonic 
axial flow conditions, even though the boundary layer 
itself is turbulent. “The reason for this is the low criti- 
cal Mach Number ( ~ 0.4) of a cylinder and the presence 
of lambda or shoulder shocks at the side of the body 
(@ = 0°) at low angles of attack. 


(3) NORMAL ForcE DISTRIBUTIONS 


Static-pressure data were obtained for the 30° cone- 


cylinder configuration at M/, = 2 for angles of attack 
of 0°, 4°, 8°, 12°, 16°, 20°, and 23°. The data are 
reduced to coefficient form, and local normal force 
coefficients per degree of angle of attack are calculated 
at each station for all conditions tested. These data 
are presented in Fig. 2 in comparison with the Van 
Dyke solution. The inadequacy of linear potential- 
flow theories to predict the normal force distributions 
at all but very low angles of attack is shown very clearly 
in this figure. 

As previously mentioned, there are two approximate 
methods (Allen’s and Kelly’s) which yield satisfactory 
results in some instances. However, the accuracy of 
these methods, to a great extent, is dependent upon 
the fineness ratio and the angle of attack range. In 
order to understand these anomalies, the potential 
normal-force-per-degree distribution is first obtained 
from the data shown in Fig. 2. Since the local normal 


force coefficient slope, shown in Fig. 2 for ag = 4°, 
begins to rise from zero approximately 8.5 body diam- 
eters aft of the nose, it may be assumed that this is a 
viscous lift and that no viscous lift is obtained ahead 
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Fic.2. The axial distribution of local normal force coefficient 
per degree of angle of attack for a 30° cone-cylinder configura- 
tion (My = 2.00). 
of this point at this angle of attack. Therefore, the 
total distribution shown for ag = 4° ahead of this 
point is representative of the potential distribution. 

The viscous normal force coefficients are then ob- 
tained as the difference between the experimentally de- 
rived potential distribution and the measured total 
local normal force coefficients. These viscous com- 
ponents are, in turn, used to calculate the axial dis- 
tribution of cross-flow drag coefficient, Ca,: These 


coefficients are presented in Fig. 3 for ag = 8°, 16°, 
and 23°. Also shown are the cross-flow drag coefli- 
cients calculated by the methods of Kelly and of Allen. 
The difference between the value shown for the infinite 
cylinder® and that obtained by Allen’s method is due 
to the correction for finite body length. The correction 
applied is that estimated for the body of fineness ratio 
13.14. It may be seen that, in general, Kelly’s method 
is in good agreement with the data for the first four 
diameters only and thereafter diverges considerably 
from the measured results. On the other hand, it is 
clearly evident that Allen’s method at best provides a 
cross-flow drag coefficient which is close to the average 
for bodies of high fineness ratio at high angles of attack 
but would, in general, fail to predict accurately the 
distribution of the viscous normal force. These re- 
sults confirm the conclusions reached from analysis of 
the resultant force data that the accuracy of both meth- 
ods is dependent upon the body fineness ratio and the 
angle of attack. 

In the absence of a rigorous theory for the prediction 
of viscous forces and moments, it appears that more 
accurate predictions than provided by the two methods 
now available may have to depend on the use of corre- 
lated experimental results. The similarity noted in 
the data presented in Fig. 3 at higher angles of attack 
suggests that the data might be correlated by dividing 
the experimental results by the steady-state two- 
dimensional cross-flow drag coefficient for a circular 
cylinder, 
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The axial distribution of cross-flow drag coefficient on 
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a 30° cone-cylinder configuration (ag = 8°) 
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Fic. 3(c). The axial distribution of cross-flow drag coefficient on 
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These data, in the form Ca,/ ca, are presented in Fig. 
4 for all conditions for which data were obtained. It 
is evident that good correlation is achieved for the 
higher angle of attack range. 

These data also show that the viscous lift is first 
developed near the base of the body and moves forward 
and increases in magnitude with increasing the angle 
of attack. It is evident, therefore, that, in order to 
predict accurately the viscous lift and the center of 
pressure of viscous lift at all angles of attack, methods 
are required which are considerably more refined than 
any predicting a viscous contribution over the entire 
extent of the body, starting from zero angle of attack. 
In particular, it is indicated that the use of any existing 
method which is found to be successful at high angles 
of attack would tend to overestimate the lift-to-drag 
ratio for body alone configurations near the angle of 
attack for maximum lift-to-drag ratio. It is felt that 
additional data and an increase in the knowledge per- 
taining to viscous flow separation may provide theo- 
retical or semiempirical methods which would account 
for all of the observed phenomena. It may be noted 
that the rise in cross-flow drag coefficient at the lower 
angles of attack appears to be associated with the sta- 
tion at which the body vortex system first becomes fully 
developed. The mechanism for the genesis of the lift 
is discussed in a later paragraph on vortex origins. 

The correlated data, as represented by the solid 
curve in Fig. 4, were used to calculate the viscous nor- 
mal force and the center of pressure as a function of 
angle of attack for the configurations tested, which 
had the same nose shape as the static-pressure model. 
These results are presented in Fig. 5. The normal force 
coefficient data presented in this figure are consistently 
slightly higher than the values predicted by the corre- 
lation curve. Further analysis revealed that at least 
half of the difference shown was directly attributable 
to the difference between the force data and the values 
obtained by integration of the surface pressure; the 
contribution to the cross-flow lift arising from the vis- 
cous shear is neglected when relying solely on the 
integration of the pressures. The remaining differ- 
ence is within the accuracy to be expected. In any 
event, it is seen from the comparisons shown that the 
correlation curve used is by far preferable to the two 
existing methods (see Fig. 3) for the purpose of pre- 
dicting the viscous phenomena for the range of the 
body fineness ratios tested. 

In order to obtain some idea of the range of Mach 
Numbers, nose shape variations, and nose fineness ra- 
tios for which such correlation curves could be used 
for the estimation of the viscous characteristics, sev- 
eral other configurations were investigated. The re- 
sults obtained for three typical cases are presented in 
Figs. 6 and 7. The selected cases are: (a) the test 


configuration at My = 1.50; (b) a configuration of the 
same nose fineness ratio as the test configuration but 
having a tangent-ogive nose shape, and (c) a configura- 
tion with a tangent-ogive nose of double the fineness 
ratio of the test configuration. 
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Mo=2.00 
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Fic. 4. The axial distribution of cross-flow drag coefficient ratio 
(Mo = 2.00). 
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Fic. 5. The normal force and center of pressure characteristics 
of 30° cone-cylinder configurations of various fineness ratio. 
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Fic. 6. The effect of Mach Number on the normal force and 
center of pressure characteristics of 30° cone-cylinder configura- 
tion (//d = 13.14). 
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Fic. 7. The normal force and center of pressure characteristics 
of two ogive-cylinder configurations (My) = 2.00). 
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(6) SIMPLIFIED THEORETICAL 
MODEL 
> 
Fic. 8. Assumed body wake flow and simplified theoretical 
model. 


Briefly, these comparisons indicate that the correla- 
tion curve presented herein is affected only slightly by 
the free-stream Mach Number and the nose fineness 
ratio, at least for the range of variables for which data 
were obtained. The data presented by Perkins and 
Jorgensen,® wherein a similar analysis was made, 
showed an effect of nose fineness ratio on the viscous 
lift distribution. This effect was observed, however, 
for a body having a nose fineness ratio more than three 
times that of the body tested in the present investi- 
gation. 


(4) Bopy FLow CHARACTERISTICS 


The wing-body configuration tests conducted as 
another part of these studies have demonstrated (as 
have other similar investigations) that appreciable 
interference effects result from the interaction of the 
vortices shed from a body of revolution at high angles 
of attack with any attached lifting surfaces. The 
nature of this interaction is well understood and 
methods! for predicting the forces acting on the wings 
or tails are available, provided the strength and loca- 
tion of the vortex are known. The major objective of 
this account of the investigation is to provide more 
basic information on this latter phase problem which 
has received less attention. 

The body vortices and associated wake system affect 
the lift of a lifting surface immersed in this flow in two 
ways: (a) alteration of the flow inclination across the 
span, and, therefore, alteration of the magnitude and 
the spanwise distribution of the lift, and (b) changes 
in the local flow characteristics—i.e., Mach Number, 
static pressure, density, etc.—with associated changes 
in the magnitude and distribution of the lift. 
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Without an analytical procedure for determining the 
strength of the body vortices and the flow variations 


they produce, it is impossible to arrive at a reasonable. 


estimate of the interference effects generated in a sys- 
tem containing such vortices. 


Theoretical Considerations 


The three-dimensional steady flow associated with a 
body of revolution at moderate angles of attack has been 
visually studied by means of schlieren and vapor screen 
photographs by many experimenters, so that sufficient 
evidence is available to allow approximation of this 
flow by the model shown in Fig. S8(a). It is also known 
from visual studies and total pressure surveys that the 
boundary layer first begins to thicken on the lee side 
of the body near the base of the body. This phe- 
nomenon progresses as the angle of attack increases, 
until the boundary layer separates and two columnar 
regions of concentrated vorticity are formed. The 
axial location at which these vortices are first formed is 
also known to move forward as the angle of attack is 
increased. 

It is reasonable to assume from these observations 
that the flow of major interest here is completely 
divorced from any potential flow effects on the nose 
of the body. With this assumption and the omission 
of any effects of the feeding vortex sheet, the more 
tractable simplified model shown in Fig. 8 (b) may be 
used, and the resultant flow is that of simple line vor- 
tices symmetrically situated on the lee side of a circular 
cylinder in a uniform stream. At this point, it is 
convenient to take advantage of the well-worked 
out properties of incompressible cross-flow (slender- 
body) theory. This recourse to an expedient simpli- 
fication is sanctioned because a combination of non- 
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Fic. 9. Schematic illustration of cross-flow potential. 
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linear effects and compressibility effects would be p_=20° 
difficult to treat. Actually, since the vortex system is 30, | 
of very small span, in the real case the delay of inter- | | 
action between the vortices and between the vortices ‘. 
and the body, due to the supersonic axial flow, should 2 
not be of too great importance. This way of dealing | | 
with the cross flow does not necessarily limit the 9 
analysis to two-dimensional considerations because, HES i 
as will be shown later, the effects of vortex motion and ~10 
the rate of change of vortex strength on the axial t 
velocity component may still be accounted for in a ~20) dvdecrex! Loc ATION T } 
“hybrid” fashion—that is, the cross-flow velocities may | jm Tal 7 | po 
be added vectorially to the axial component of the — 10 20 30 40. “10 20 30 40 
free stream, and, in addition, any axial perturbation «fe 
velocities due to vortex growth and motion, as pre- Fic. 10. Experimental total pressure ratio contours and local 
dicted by potential theory, may be superimposed. 
y Pp y be sup Pp 
locity. These velocities are related to the potential 


Basic Formula 


The complex potential for the cylinder and two vor- 
tices in a uniform stream, together with two image 
vortices placed within the body in order to satisfy the 
boundary condition of tangential flow at the body, is 
given by Lagerstrom and Graham’ as 


W= = —iVo[s — (a?/s)] 
(iT /2m)In { [(s — (s + a?/o)]+ 
[(s + @) (s — a?/a)]} (1) 


The first term is simply the complex potential for the 
flow past a circular cylinder at the origin, and the 
second term gives the potential for four vortices; two 
of strength T located at o, and —a?/o, and two of 
strength —I located at —é and a*/é. A schematic 
illustration of such a flow is shown in Fig. 9 and serves 
the nomenclature used therein. The motion of the 
vortices in this flow may be determined simply by trac- 
ing the motion of one real vortex, since the other real 
vortex is symmetrically located, while the image vor- 
tices must be located so that the aforementioned 
boundary condition is always satisfied. For the vortex 
located at o, the velocities are given by the expressions 


dé/dt = — Vo[(2a*én)/(E? + + (T'/2m) X 
+ 9? — — + 9? — J+ 
+ 9? — a*)? + (2a) 
dn/dt = + — + — 
((1/28) + + 0? — a*)] — 

+ 9? + a*)]/[(? + — + (2b) 
The path of the vortex in three dimensions may be 
expressed by using the following relation for the axial 
velocity : 


dz/dt = Uo (3) 

and therefore 
di/dz (dé/dt)/(dz/dt) = (1/U ) (dé/dt) (4a) 
dn/dz = (dn/dt)/(dz/dt) = (1/Uo) (dn/dt) (4b) 


The cross-flow velocities are u and v, which are the 
velocity components in the x and y directions in the 
cross-flow plane, and w is the axial perturbation ve- 


equation in the following way: 
u = = RI[OW/dx] (5a) 
v = = RI[OW/cy| (5b) 


The aforementioned hybrid combination of the cross- 
flow velocities with the axial velocity component to 
obtain the resultant local velocity, U;, may now be 
expressed as 


U,? = (Uo cos ag + w)? + u? + v? (6) 


Using the adiabatic relationships for the flow of a 
perfect gas, the calculation of the other flow properties 
may be easily accomplished. In particular, a basic 
property, the local Mach Number, is given by 


= + [(y — 1)/2]Mo X 
[1 — (Us/Us)?]} (7) 
It should be noted that the axial perturbation ve- 
locity, w, arises from the fact that both the strength of 
the vortices and their positions are functions of the 
axial dimension z. The equation expressing this rela- 
tionship is 
w= RI[OW/dz] = + 
(OW/0a) (da/dz)| (8) 


The terms (dI'/dz) and (do/dz) may be evaluated by 
using the relationships for the variation of vortex 
strength and the vortex trajectory. The latter is dis- 
cussed above and a discussion of the former follows. 


Vortex Strength 


Consideration of the forces acting on the cylinder 
and the changes in momentum associated with the 
axial development of the body vortices for the simpli- 
fied flow model of Fig. 8(b) results in a simple rela- 
tionship between the strength of the vortices and the 
viscous normal force acting on the cylinder. 

The simplest way to arrive at this relationship is to 
consider the three-dimensional flow as an unsteady 
two-dimensional flow—i.e., ¢ = 2z/U>. For the as- 
sumed model, the impulse applied to the cylinder in 
time ¢ is given by the following, wherein ¢ = 0 denotes 
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Fic. 11. Circulation strengths of concentrated vortex and feed- 
ing sheet versus angle of attack. 


any point lying ahead of the point of origin of the vor- 
tices and viscous normal force; 


f (AN,/dz) (dz/Us) = (Nw/U) (9) 
0 0 


This impulse is equal to the change in momentum of 
the flow from time ¢ = 0 to time ¢. It is evident from 
the choice of the station for ¢ = O that the change in 
momentum of the fluid of the theoretical model is that 
change associated with the momentum of the vortices 
and the image vortices. Therefore, since it can be 
shown’ that the momentum of two vortices is equal 
to the product of the vortex strength, density, and the 
span of the vortices, the change of momentum of the 
fluid is given by 


Change in momentum = 2pI'(é — &;), (10) 


Thus, the finally sought result is achieved by equating 
the impulse to the change in momentum, or 


N,,/ Uo = 2oo0 — (11) 
so that with some rearrangement and substitution 
= — &)2] (12) 


It should be noted that (Cy,), is the viscous normal force 
coefficient for all the viscous normal force ahead of 
station z. 


Body Wake Surveys 


Two experimental investigations pertaining spe- 
cifically to the details of the body flow have been con- 
ducted as a part of this general study. One part of the 
part of the flow-field survey consisted of an exploratory 
study in which static and total-pressure measurements 
were taken circumferentially at various radial dis- 
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tances around a cone-cylinder body at AJ = 1.50 and 
2.00 for two body stations and at angles of attack from 
0° to 23°. The total-pressure contours thus obtained 
serve to indicate an approximate vortex location; 
this information, together with estimates of the vis- 
cous lift distribution, allows calculation of the vortex 
strength by means of the relationship previously de- 
rived. The values of vortex strength thus calculated 
are used to compute the static and total pressures on 
the basis of the adiabatic relationships described in the 
preceding section. Normal shock relationships are also 
applied in obtaining the pitot pressures. The close 
agreement of the measured pitot pressures over the 
greater portion of the survey plane with the predicted 
values calculated in the above manner provided the 
impetus to conduct a second investigation. 


The second set of wind-tunnel tests in the series was 
a more complete investigation in which calibrated cone 
probes and pitot-pressure probes were used for the sur- 
vey instrumentation. With this calibrated instru- 
mentation it was possible to define completely the true 
local flow conditions, including flow inclination, within 
the range of the calibration. The calibration itself 
was rather.extensive and the ranges through which the 
major parameters were varied are given in Table 1. 


It should be noted that, since the flow-inclination 
instrumentation consisted of conical probes provided 
with four orifices (orifices located 90° apart), and since 
the probes were calibrated for the complete range of 
roll-orientations, Oprobe, it was possible to obtain both 
components of the local flow inclination—1.e., down- 
wash and sidewash. 


Vortex Strengths 


The major purpose of this vector-field investigation 
was to confirm the hypothesis, by detailed measure- 
ments of the flow field, that the strengths of the con- 
centrated body vortices could be estimated from the 
viscous normal force distributions and vortex locations. 
It was also intended to provide more basic knowledge 
of the body flow, the nature of the vortex feeding-sheet, 
and the vortex starting phenomenon. 


The general nature of the body-alone flow is shown 
by the total-pressure-ratio contour plot and flow- 
inclination survey plot presented in Fig. 10 for ag = 
20°. The total-pressure plots are most useful for lo- 
cating the vortex core and the feeding vortex sheet in 
the survey plane. The flow inclination results are 
presented in vector form where the tail of the vector is 
located at the point at which data were obtained. The 
magnitude and direction of the vector correspond to 
the magnitude and direction of the local flow inclina- 


TABLE 1 
Range of Calibration of Instrumentation 
Minimum Maximum 
M 1.50 2.50 
€ 0° 45° 
Oprobe 0° 45° 
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tion. For the sake of clarity, all the data* are not 
shown in these figures. 

These plots are typical of the data obtained for all 
cases for which a region of concentrated vorticity was 
in evidence. The total-pressure plots are used to de- 
termine the path of integration used in the velocity- 
perimeter integration performed to obtain the vortex 
strength, [. The contour of integration is chosen to 
consist of rays and arcs along which data were obtained. 
Therefore, the vortex strength is calculated using the 
following expression: 


+ v,dr +f verid@ (13) 
02 


In practice, the length of the integration path be- 
tween the limits of integration was equally divided by 
the points at which data were obtained and a finite- 
step summation process was substituted for the indi- 
cated integrations. Several closed-contour integra- 
tions were made for each case to ensure the determi- 
nation of the true vortex strength and feeding-sheet 
strength. For some cases, it was not possible to com- 
pute the vortex strength without resorting to an extrap- 
olation procedure. That is to say, when the vortex 
lies close to the body and use of the body surface for 
the inboard circumferential contour is required, the 
measured variation of the tangential velocity com- 
ponent, v, is extrapolated to the body surface for use 
in the integration procedure. 

The calculated values of the circulation strength of 
the body vortices and of the vortex feeding sheets are 
presented as functions of angle of attack in Fig. 11. 
It may be noted that, in accordance with expectations, 
the circulation strengths of the vortex and feeding sheet 
increase with both increasing angle of attack and in- 
creasing body station. 

The experimental viscous normal force distributions 
discussed in the preceding section and the vortex loca- 
tions determined from the total pressure ratio con- 
tours are used to compute the circulation strength 
of the concentrated vortices for the various body 
stations by Eq. (12). These predictions of concen- 
trated vortex strength are presented in Fig. 12 in 
comparison with the values obtained from the closed- 
contour integrations of the cross-flow velocities deter- 
mined by the detailed flow measurements. It is seen 
in this figure that the theory predicts values which, 
in general, are in close agreement with the measured 
results. The cases that are not in satisfactory agree- 
ment with the theoretical predictions are those cases 
representing conditions near the vortex starting point 
for which both normal force coefficient, Cy,» and the 


vortex span difference, ( — &,), are very small. It 
has been found that for these conditions the ratio of 
*In general, for the sake of brevity, only typical selections 
from the data obtained are presented. A more complete docu- 
mentation is presented in reference 8. 


these two small numbers can vary considerably with 
changes in either or both of the above-mentioned 
quantities, even though the changes are within the 
accuracy of the data. 

The experimentally determined values of concen- 
trated vortex strength are generally smaller than the 
theoretical predictions, as attested by the comparison 
presented in Fig. 12. It is believed that this trend may 
be due primarily to the fact that the strength of the 
feeding sheet has been ignored while all of the viscous 
normal force has been used in Eq. (12). The use of 
a more general equation based on a finite number of 
discrete vortices in order to represent the feeding sheet 
may improve the agreement shown. It was felt, 
however, that there was not sufficient data available 
to warrant such a calculation, in view of the assump- 
tions that would have to be made concerning the rela- 
tive strengths and locations of the feeding-sheet vor- 
tices. Another reason for the generally higher theo- 
retical prediction may be the supersonic nature of the 
axial flow. That is, the viscous normal force at any 
axial station may be related more accurately to the 
vortex strength existing further downstream in accord- 
ance with propagation processes at work in supersonic 
flow. In any event, any such effects not considered 
in the elementary analysis given here are indicated to 
be of higher order. 


Local Flow Inclinations 


Calculations were made of the local cross-flow veloci- 
ties using the experimentally determined vortex 
strengths and locations for several cases for which 
concentrated vortices are present in the body flow. 
The cross-flow velocities were calculated by Eqs. (5a) 
and (5b), which reduce to the following specific forms 
of two concentrated vortices symmetrically positioned 
above the body: 
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Fic. 12. Comparison of theoretical concentrated vortex 
strengths with concentrated vortex strengths obtained from ex- 
perimental data. 
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Fic. 13. Comparison of experimental and theoretical flow in- 
clinations (ag = 8°) 
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clination distributions (ag = 8°). 


40 


3.0 


20 


Fic. 14. 


Fic. 16. 


OF THE AERO/SPACE SCIENCES—MARCH, 


10 20 30 40 


Mo=200 

| \exPeriment *°, \\ THEORY 
\\ 

ead 


10 20 30 49 


Comparison of experimental and theoretical flow in- 
clinations (ag = 20°). 
Mo=2.00 ap=20° 3/£=6.57 
ANGLE THEORY EXPERIMENT 
a ° 
+ ,o. 
6=0° a 
2 


N) 
° 


& OR @ DEG 


Comparison of theoretical and experimental flow in- 
clination distributions (ag = 20°) 


Rig 
€, is 
W, 
cor 
loci 
any 
ver} 
the 


the 

the: 
fron 
rela' 


It 
prov 
and 
(usu 
and 
sons 
for t 
trate 
spec’ 
wher 
whic 
lated 


corre 
flow 
rays 
as 5 
clarit 


} | wh 
20 4 ] 

q 

was 
cal 

RY -20 0 20 

++ MH 3.0 

Ala | |_| 

1.0 

-20 20 “20 20 


INVESTIGATION OF NORMAL FORCE DISTRIBUTIONS 165 


u/ Vo = —[2a*xy/(x? + y?)?] — X 
— — ((y — 1) + 
— 2)/Rs*] — — (14a) 


= 1 + — y?)]/(x? + + 
{ — — — &)/R*] + 
+ &:)/Rs?] — +  (14b) 


where 


(x — &)*? + (y ni)” 
= (x + &)? + 
R? = (x + §)*? + (y — »)* 


and (£, 7) and (&;, ;) are the coordinates of the right- 
hand vortex and image vortex, respectively. 

These velocities are used to obtain the resultant cross- 
flow velocity, V, and then the local flow inclination, e, 
and the direction of the flow, 6. The relationships 
used in these calculations are: 


(V/Vo)? = (u/Vo)? + (v/ Vo)? (15) 
= tan~! [(V/Vo) tan ag] (16) 
6 = [v/u] = tan— [(v/Vo)/(u/Vo)] (17) 


Rigorously speaking, the local flow inclination angle, 
e, is also influenced by the axial perturbation velocity, 
w, given by Eq. (8). However, it is felt that the added 
computational effort involved in calculating this ve- 
locity is too great to be practicable in most cases. In 
any event, sample calculations have shown that w is 
very small relative to lL) cos ag for the greater part of 
the flow field. 

Modification of the cross-flow velocity relationships 
was made in the region of the vortex core. The classi- 
cal assumption that the core rotates as a solid leads to 
the following relationship applied to the region within 
the radius of the core, r,, the extent of which is estimated 
from the experimental flow inclination data. These 
relationships are: 


u/Vo = —(T/2rVo) [(y — n)/r.?] (18a) 
v/Vo = (T/2rVo) [(x — &)/r,?] (18b) 


It should be noted that these relationships do not 
provide for a transition region between the viscous 
and potential flows. The omission of such a region 
(usually very small) greatly simplifies the calculations 
and has little effect on the overall results. Compari- 
sons between the theoretical and experimental results 
for two typical cases, one without and one with concen- 
trated vortices, are presented in Figs. 13 and 14, re- 
spectively. The results are presented in vector form 
where the tail of the vector is located at the point for 
which the flow inclination was measured and calcu- 
lated. The magnitude and direction of the vector 
correspond to the magnitude and direction of the local 
flow inclination. Although data were obtained along 
rays with sector angles chosen in increments of as small 
as 5° (A@ = 5°), the comparisons, for the sake of 
clarity, are made for circumferential increments of 


15°. The theoretical values are presented for the loca- 
tions corresponding to those for which flow inclination 
measurements were made. 

Additional comparisons are presented in Figs. 15 and 
16 for several selected circumferential rake locations 
(@ = O°, 45°, 60°, 75°). The conditions selected 
pertain to the survey station at z/d = 6.57 for angles 
of attack of ag = 8° and 20°. These latter compari- 
sons are for the components of the local flow inclina- 
tion, @ and 8, where a and @ are the local angles of 
inclination which correspond to the tangential and 
radial velocity components, v. and v,, respectively (see 
Fig. 9). 

In general, the agreement between the theoretical 
results and measured data is as good as may be ex- 
pected in view of the fact that the theoretical methods 
used ignore mauy facets of the problem. In particular, 
it is noted that the two major regions for which the 
agreement between theory and experiment deteriorates 
with increasing angle of attack are near 6 = 0° and 
45°. For @ = 0° this discrepancy may be attributed 
to increasing effects of compressibility with angle of 
attack; this trend is exhibited most severely at the side 
of the body. For @ = 45° the increase in the circula- 
tion strength of the vortex feeding sheet with increas- 
ing angle of attack is primarily the cause for the noted 
variation in the flow inclination. The omission of the 
effects of the feeding sheet at the higher angles of attack 
is, therefore, the primary reason for the lack of agree- 
ment in this region. These effects will be discussed 
further in the following sections. 

Agreement between experiment and theory was also 
noted to deteriorate with increasing axial distance from 
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Fic. 17. The effects of compressibility on the flow inclination 
distribution at various angles of attack. 
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Vic. 18. Typical effects of the vortex feeding sheet on the flow 
inclination distributions (ag = 20°). 


the nose. In addition to the aforementioned phe- 
nomena, it is believed that the increasing effects of 
viscous flow and the elongation and spread of the 
concentrated vortex region have a direct bearing on 
the agreement to be expected. 


Compressibility Effects 


In the calculation of the cross-flow velocity distribu- 
tions, compressibility effects can to a certain extent be 
accounted for in a simple manner, at least for the re- 
gion to the windward of the cross-flow sonic line—i.e., 
the locus of points at which M, sin e = 1.00. For this 
region it may be assumed that the effects of the vortices 
are not felt, so that the use of the compressible poten- 
tial, as given by Kaplan® for the cylinder-alone in a 
uniform stream, is therefore the best representation of 
the flow. This procedure has been applied for 6 = 0° 
to the previously presented data and other cases, and 
the results are presented in Fig. 17. Also shown for 
purposes of comparison are the incompressible theo- 
retical results for the intermediate station which ap- 
proximate the incompressible results for the other two 
body stations. It is seen that, in general, the afore- 
mentioned procedure results in some improvement in 
the agreement of the theorctical results with the experi- 
mental data. 


Vortex Feeding Sheet Effects 


The omission of the vortex feeding sheet in the cross- 
flow potential used to estimate the cross-flow velocity 
distribution naturally results in appreciable disagree- 
ment between theory and experiment in the vicinity 
of the sheet. It would be possible to account for the 
effects of the vortex sheet if the precise orientation and 
distribution of circulation strength along the sheet were 
known. Unfortunately, the probe spacing of the 
instrumentation used was not sufficiently small to 
allow a practical assessment of these physical char- 
acteristics of the sheet. However, since it is possible 
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Fic. 19. Ratio of the circulation strength of the concentrated 


vortex to the total circulation strength versus viscous normal 
force coefficient. 


to determine the total circulation strength of the sheet 
and since the approximate locus of the sheet is indi- 
cated by the total pressure contour plots, calculation 
of the effects of the vortex sheet was attempted for a 
typical case. 

The effects of the vortex feeding sheet were estimated 
by including in the mathematical model a number 
(four) of small discrete vortices and the respective 
image vortices for each sheet. The location of each 
vortex was chosen so as to provide, as nearly as possible, 
equal vorticity for each segment of the sheet. The 
total circulation strength of the sheet was equally 
apportioned among the vortices used to simulate the 
sheet. The results of the calculation are presented in 
Fig. 18. Itis evident that for the case chosen the agree- 
ment between theory and experiment is substantially 
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improved in the vicinity of the vortex feeding sheet 
by accounting for the presence of the sheet in the man- 
ner discussed above. 

The vortex feeding sheet effects were not investi- 
gated further due to the approximate nature of the 
method used and due to the fact that the computational 
work required is nearly proportional to the number of 
vortices included in the cross-flow potential. In addi- 
tion, theoretical methods are not available for the pre- 
diction of the circulation strength of the vortex sheet. 
The circulation strength of the sheet, however, may be 
estimated empirically from experimental results, such 
as are available from this investigation. For example, 
the ratio of vortex strength to total (vortex plus feed- 
ing sheet) circulation strength is given in Fig. 19. 
Such a plot allows the determination of the strength 
of the sheet as a function of viscous lift coefficient ahead 
of the vortex axial station under consideration. It is 
evident that the data are well correlated in this man- 
ner. Additional experimental work may provide a 
broader correlation of this type. 


Vortex Paths 


The experimentally determined vortex locations for 
the configuration tested are presented in Fig. 20. The 
calculation of the vortex paths starting from a known 
vortex location was attempted by the use of Eqs. (2) 
and (4), together with the knowledge of the normal 
force distributions on the body. In general, contrary 
to the experience of Jorgensen and Perkins, !° the vortex 
trajectories could not be adequately predicted. Ex- 
amination of the possible reasons for this anomaly has 
indicated that the vortex trajectory is very sensitive 
to the accurate determination of the vortex location 
selected for starting the trajectory, although not very 
sensitive to vortex strength. This latter result con- 
curs with the findings of Jorgensen and Perkins.’ It 
is, therefore, reasonable to assume that the major rea- 
son why it was not possible to predict adequately the 
vortex trajectories in the present investigation is the 
lack of ability to measure the vortex locations within 
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cumferential variation of the static pressure coefficient. 


tolerable limits with the probe spacing used, but it is 
not attributable to inadequacy of the theory. 

Although a trajectory could not be adequately pre- 
dicted, the vortex locations presented are sufficiently 
accurate for calculation of local flow properties. For 
this reason, any correlation of the results obtained is 
desirable. The vortex spanwise coordinates, £/a, 
shown in Fig. 20, are all noted to approach a value of 
0.750, except near the vortex starting location. It is 
indicated, therefore, that £/a = 0.750 is a satisfactory 
spanwise coordinate for most computational purposes. 
An attempt was made to correlate the vortex coordinate 
n/a as a function of the angle of attack and the viscous 
lift acting on the body ahead of the axial station in 


question. It is found that the expression 


n/a = (22/d) tan (kag) (19) 


collects the data obtained in the present investigation, 
where is a function of the viscous normal force coeffi- 
cient, Cy,,.. This correlation is presented in Fig. 21. 


Vortex Origins 


The prediction of the vortex origins, like that of 
many other viscous phenomena, is not possible by 
theoretical methods at the present time. Analysis of 
the static pressure distributions on the body for all 
conditions tested indicates that the vortex formation is 
associated with the station at which a second minimum 
first appeared in the pressure coefficient curve (C, ver- 
sus 6). 

This phenomena is illustrated by the pressure coeffi- 
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Fic. 23. Region of indicated vortex origins versus angle of attack. 


cient data presented in Fig. 22. The data presented 
are typical of that for all conditions tested. It is seen 
in the data for these adjacent stations that the curve 
for the more forward station does not have a second 
minimum while the curve for the more aft station has 
a second minimum at 6 = 50°. This second minimum 
in the pressure coefficient data is to be expected if a 
vortex bubble is formed on the body at that axial 
location, as illustrated by the sketch in Fig. 22. It is 
evident, therefore, that a vortex bubble first was 
formed at some station intermediate to those for which 
data have been plotted in this instance. An area of 
incipient vortex formation, defined in this manner, is 
therefore indicated for the conditions tested. This 
region is presented as a function of angle of attack in 
Fig. 23. Included in this figure is the band presented 
in reference 10 as the region of vortex origins for an 
ogive-cylinder body. The region of vortex origin for 
both bodies is seen to be quite similar. 

The relationship of the vortex origins to the axial 
development of viscous lift is indicated by the symbols 
denoting a value of constant viscous normal force 
coefficient (Cy,, = 0.175). It is seen that the station 
at which this value of viscous normal force coefficient 
is achieved is in fair correspondence with the axial 
location of the vortex origins. 


(5). CONCLUSIONS 


(1) The axial distribution of cross-flow drag coeffi- 
cients for bodies of revolution at high angles of attack 
in supersonic flow in general will differ considerably 
from the distributions predicted by the methods of 
Allen’ and Kelly. 

(2) The use of correlated experimental results for 
the axial distribution of cross-flow drag coefficients in 
conjunction with the hybrid theory of Van Dyke? pro- 
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vides a satisfactory method for estimating the normal 
force and pitching moment characteristics of similarly 
shaped bodies of revolution. Furthermore, such cor- 
related results can be extended, within reasonable 
limits, to provide satisfactory estimates of the viscous 
force and moment characteristics for configurations 
and Mach Numbers other than those for which data 
are available. 

(3) It has been confirmed by detailed flow meas- 
urements that the circulation strengths of the concen- 
trated vortices in the body vortex wake flow can be 
estimated, provided the viscous normal force dis- 
tributions and vortex locations are known. 

(4) The flow inclination distributions in the flow 
field around a body of revolution, in general, can be 
closely predicted by the theoretical potential cross- 
flow model consisting of a cylinder and two symmetric 
vortices (plus image vortices) in a uniform stream. A 
closer agreement with experiment can be achieved in 
the vicinity of the vortex core, in the vicinity of the 
vortex feeding sheet, and in the regions affected by 
cross-flow compressibility by accounting for these phe- 
nomena with certain simple refinements to the basic 
theoretical model. 

(5) Several of the physical properties of the body- 
alone wake vortex flow are found to correlate in a 
manner dependent primarily on the distribution of vis- 
cous lift on the body. The present methods of corre- 
lation should be applied to a wider range of body 
configurations and Mach Numbers to determine their 
limitations. 
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On the Flow of a Sonic Stream Past an 
Airfoil Surface’ 


COLIN S. SINNOTT* 
National Physical Laboratory, Teddington, England 


SUMMARY 


This study of the flow about an airfoil in a near-sonic stream 
indicates the important factors determining the pressure dis- 
tribution on the airfoil. Analysis of the Mach wave pattern 
suggests that the supersonic domain of the flow can be derived 
from two simple-wave flows, one arising from the Mach waves 
reflected at the sonic line and the other from the changes in air- 
foil surface slope. The compressive effect of the reflected Mach 
waves is determined quantitatively as a function of airfoil lead- 
ing-edge geometry from an analysis of measured pressure dis- 
tributions for uncambered airfoils; and it is shown how this can 
be superimposed on the wave system from the curved surface to 
give an equivalent simple-wave flow over the airfoil. 

An application of this scheme to the calculation of the pressure 
distribution over an airfoil in a sonic stream gives results in good 
agreement with experiment. 


SYMBOLS 
c = airfoil chord 
H = stagnation pressure 
M = Mach Number 
p = static pressure 
x,y = coordinates of airfoil surface 
X =(x — xcr)/(xre — Xcr) 
a = airfoil incidence 


bre = flow deflection at trailing edge 
6 = surface slope measured from airfoil crest 
Prandtl-Meyer angle corresponding to local Mach 


o = 
Number on an airfoil surface 
2 = effective compression angle 
= (x/c); ie., a function of airfoil 
geometry near the leading edge 
Subscripts 
0 = free-stream conditions 


TE = conditions at trailing edge 
CR = conditions at crest 


(1) INTRODUCTION 


i is FLOW of a sonic stream past an airfoil surface 
is of considerable interest, for it is known that the 
resultant pressure distribution on the surface is insen- 
sitive to small changes in free-stream Mach Number. 
Because of this, an estimate of the sonic pressure dis- 
tribution over an airfoil can be a useful guide to the air- 
foil characteristics throughout the transonic speed 
range. Some typical results illustrating the phenom- 
enon are given in Fig. 1, where the pressure dis- 
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tribution over an airfoil is shown for a number of free- 
stream Mach Numbers. It is seen that once a shock 
wave is established in the flow the pressure distribution 
head of it changes little with increase in Mach Number 
and the consequent rearward movement of the shock 
wave. When the shock wave has reached the trailing 
edge of the airfoil surface the pressure distribution over 
the whole surface is insensitive to further increase in 
free-stream Mach Number until some slightly super- 
sonic value is attained. Except for very thin airfoils 
the surface shock wave will reach the trailing edge at a 
high subsonic free-stream speed, so that the range of 
Mach Number over which the pressure distribution is 
approximately constant will include the sonic value; 
this has been called the “sonic range.’ The associated 
independence of local Mach Number and free-stream 
Mach Number is known as the ‘‘Mach Number freeze.”’ 
Tests on thin airfoils suggest that the surface shock 
does not reach the trailing edge until a slightly super- 
sonic free-stream speed is attained. In these cases, 
the development of the flow is not essentially different 
from that described above, for although the sonic-range 
will be entirely supersonic the Mach Number freeze will 
still be established at points upstream of the shock wave 
at lower speeds. 

A number of approximate methods of calculating 
sonic-range pressure distributions have been suggested, 
and purely theoretical results have been obtained for 
double-wedge and flat-plate profiles. The semiempiri- 
cal method of Lindsey and Dick! requires a knowledge 
of the sonic point position, the pressures from there to 
the trailing edge being calculated by simple-wave 
theory with a modified expansion rate. Gullstrand? 


olK 


Fic. 1. Typical variation of pressure distributions with free- 
stream Mach Number for a 10 per cent thick airfoil at zero 
incidence. 
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Expansion waves ------- » Compression waves 


Continuation of the compression waves after reflection 
at the aerofoil surface is omitted for clarity 


Fic. 2. Mach wave pattern of the sonic flow past an airfoil 
surface. 


and Bévierre* employ exact simple-wave theory from a 
chosen position on the airfoil surface and treat the up- 
stream flow by other methods. There are important 
criticisms of these approaches; first, it is not yet pos- 
sible to determine the sonic point on an airfoil with use- 
ful accuracy, and second, modification of the super- 
sonic simple-wave expansion rate does not represent 
the physical mechanism whereby the leading-edge sub- 


Expansion waves 


Compression waves 
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Fic. 3. Illustration of how the compression effect over the 
rear part of an airfoil arises from an expansion ahead of the 
crest. 
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sonic region affects the following supersonic flow. Also, 
use of exact simple-wave theory is physically incorrect. 
It is felt that for these reasons none of the above meth- 
ods could be applied with confidence. The theoretical 
solutions for a double-wedge‘ and a flat plate’ appear 
satisfactory but cannot be directly extended to conven- 
tional airfoil sections. 


In this paper, an alternative semiempirical method is 
suggested in which results obtained from an analysis of 
measured sonic-range pressure distributions are used in 
conjunction with simple-wave theory modified to rep- 
resent the actual Mach wave pattern of the flow, the 
modification being partly derived from the theoretical 
solution for a double-wedge. By this means, the major 
part of the supersonic flow over an airfoil surface can be 
estimated. No method for calculating the pressure dis- 
tribution near the leading edge is given, but in many 
cases the calculated distribution can be extrapolated 
forward to the stagnation point with sufficient accu- 
racy. An example of the application of the scheme is 
given and the results agree well with measured values. 


In Section (5) there is some discussion of the physical 
mechanism whereby the sonic-range flow field is estab- 
lished, attention being given to the relative importance 
of the leading and trailing edge in influencing the pres- 
sure distribution on an airfoil. 


(2) DESCRIPTION OF FLOW PATTERN 


The flow pattern associated with the near-sonic flow 
about an airfoil surface is illustrated in Fig. 2. Condi- 
tions upstream of the leading edge are determined by 
the subsonic or supersonic nature of the free stream. 
In the subsonic case, there is a smooth deceleration to a 
stagnation point near the leading edge; alternatively, 
if the upstream flow is supersonic, there must first be a 
discontinuous deceleration to subsonic flow via a shock 
wave. However, in either case the flow immediately 
ahead of the airfoil is nonuniform normal to the free- 
streain direction. The flow on the airfoil surface, ac- 
celerating from the stagnation point, reaches local sonic 
speed at some point. Similarly, sonic speed is reached 
on streamlines adjacent to the surface, the locus of the 
sonic points forming the sonic line. Because of the 
nonuniformity of the flow the sonic line is curved, as 
are the expansion waves which originate from the curved 
airfoil surface bounding the supersonic flow behind the 
sonic line. Some of the expansion waves which meet 
the sonic line and are reflected as compression waves 
reach the airfoil surface, and so give an apparent reduc- 
tion in the Prandtl-Meyer simple-wave expansion rate 
over the rear of the airfoil. Each of the compression 
waves reaching the surface is again reflected regularly 
as a compression; hence, the magnitude of the compres- 
sive effect is double that of the expansive effect which 
gave rise to the original expansion wave. It is probable 
that some multiple reflections occur in the neighborhood 
of the sonic point. The supersonic flow over the airfoil 
terminates in a shock wave at the trailing edge; the 
significance of this shock and the associated flow de- 
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flection in relation to the flow downstream of the trail- 
ing edge has been discussed in reference 6 and is recon- 
sidered in the present context in Section (5) of this 
paper. 

In a sonic stream, the Mach wave originating from 
the airfoil surface and meeting the sonic line at in- 
finity defines the ‘‘limiting Mach line’’; clearly dis- 
turbances behind this line cannot affect the upstream 
flow. Consideration of the Mach wave pattern suggests 
that the limiting Mach line must leave the airfoil sur- 
face ahead of the point where the surface tangent is 
parallel to the free stream. This point has been called 
the crest of an airfoil, and as the position of the limit- 
ing Mach line cannot be determined easily, the fol- 
lowing discussion and analysis will be confined to the 
flow behind the crest. This flow is entirely supersonic 
in character and can be treated as arising from the 
superposition of two simple-wave flows, representing 
the two families of characteristics in the classical treat- 
ment of supersonic flow. One is due to changes in air- 
foil surface direction and the other to the compres- 
sion waves coming from the sonic line. This secondary 
wave system is illustrated in Fig. 3, where it is shown 
that the origin of the compression waves which meet 
the surface between its crest and trailing edge is a 
wedge of expansion waves arising from the curved air- 
foil surface ahead of the crest. In Fig. 4 the contribu- 
tions of the two simple-wave flows to the pressure 
variations of the surface of an airfoil are illustrated. 
Fig. 4(a) shows the pressure distribution for a double- 
wedge section; here the supersonic expansion takes 
place at a single point (the shoulder of the wedge) and 
there is no further change in flow direction until the 
trailing edge is reached. The pressure changes between 
the crest and trailing edge can arise only from compres- 
sion waves emanating from the sonic line. Fig. 4(b) 
shows how the simple supersonic flow over a curved air- 
foil surface is modified in the same way, but in this case 
part of the pressure changes are due to the expansion 
round the curved surface of the airfoil between its 
crest and trailing edge.* 

The factors determining the magnitude of the com- 
pression effect must be the shape and position of the 
sonic line, which is itself dependent on the neighboring 
surface shape, the airfoil shape in the region where the 
expansion compression waves originate, and the change 
in surface direction between the crest and trailing edge, 
for this will affect the path of the reflected Mach waves. 
In the next Section it is shown empirically how the 
total compression between crest and trailing edge can 
be estimated, and that the theoretical distribution of 
the compression for a double-wedge‘ can be used for 
curved airfoils. From these results it is shown in Sec- 
tion (4) how the compression effect can be applied as a 
correction to airfoil surface slope distribution and an 
equivalent simple-wave flow obtained. 


* In this case, of course, the flow from the sonic point onward is 
affected by the interaction of the sonic line and Mach waves, 
but, as has been stated, to simplify the problem only flow down- 
stream of the crest is considered here. 
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(3) DERIVATION OF SEMIEMPIRICAL METHOD FOR 
ESTIMATING SONIC-RANGE PRESSURE DISTRIBUTIONS 


In the description of the sonic-range flow pattern 
about an airfoil surface it is noted that the flow over the 
surface terminates in an oblique shock wave at the 
trailing edge. This feature has been investigated by 
Holder.* He found that the local Mach Number ahead 
of the trailing-edge shock appeared to be directly re- 
lated to the flow deflection at the trailing edge (deter- 
mined on the assumption of no downwash). This 
result is illustrated in Fig. 5, where the Prandtl-Meyer 
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Fic. 6. Empirical relation for the flow at the crest of 
an airfoil surface in a sonic stream (for definition of symbols 


see Fig. 5). 
Fig. 7 (a) 
¥ 
Wre Fz | 
vo <|A oY 4 
x | 
0-2 0-4 0-6 08 10 
6% RAE 104 a=-03 0 NACA x 
1-7 9 
3:7 6 22 
4% RAE 104 a= 0 + 10%RAE 102 a=-0:13 A 
23 < 
10% Oouble - wedge (theory) © 
Fio.7(b) 
os 
0-2 04 06 0-8 1-0 
Fic. 7. Distribution of compression effect. 


angle corresponding to the trailing-edge Mach Num- 
ber is plotted against the assumed flow deflection for a 
number of airfoils at incidence. The mean curve is 
shown by Holder to be similar to that given by the 
oblique-shock relation for the assumed flow deflections 
and a downstream Mach Number of about 1.08, the 
remaining pressure recovery to free-stream conditions 
taking place along the wake. The relation between 
Holder’s result and those of the present paper will be 
considered later. 

Consideration of the Mach wave pattern has sug- 
gested that the effect of the waves reflected back to the 
airfoil ahead of its trailing edge can be represented as a 
correction to the surface slope distribution. It will be 
convenient to consider the compression effect of this 
wave system in terms of the equivalent Prandtl-Meyer 
flow deflection angle, as this can most easily be related 
to the changes in surface slope. It is noted above that 
the Prandtl-Meyer angle (wz) just upstream of the 
trailing-edge shock wave is related to the trailing-edge 
flow deflection (67z)—which is also the change in sur- 
face direction between the crest and trailing edge. 
Hence, to determine the total compression angle from 
crest to trailing edge, given by [irz — (wrx — “«cr)|, 
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it will be sufficient to determine wee. This should be a 
function of the upstream geometry of the surface. 
Measured values of wer for a number of airfoils were 
plotted against a mean of the surface curvature ahead 
of each airfoil crest. There appeared to be a con- 
sistent relation for each airfoil at different incidences, 
and the introduction of airfoil leading-edge radius as a 
further parameter unified the results for all the airfoils. 
The most consistent correlation is shown in Fig. 6; 
here wer is plotted against the parameter x, which is 
defined as at x = 0.02c)/ [A(x/c) for @ = 2° to@ = 
4°] where @ is the surface slope measured forward from 
the crest. This change in x/c for a prescribed change 
in surface slope is a convenient measure of surface 
curvature. The value of y/c at x = 0.02c is preferred 
to the nominal leading-edge radius as it takes account 
of the variation in leading-edge angle of biconvex and 
similar profiles. 

From Figs. 5 and 6 the total compression from crest 
to trailing edge can be estimated; there remains the 
problem of the determination of the distribution of 
this compression. It has been emphasized that the 
pressure distribution over this part of an airfoil will be 
partly determined by the changes in surface direction 
and partly by the Mach waves reflected at the sonic 
line back onto the surface. To dissociate these two 
eifects, measured pressure distributions have been an- 
alyzed to give the increment in equivalent Prandtl- 
Meyer angle over we, at positions along the airfoil 
chord. The difference at a given point between this 
quantity (w — wer) and the surface slope is then the 
compression angle measured from the crest and will be 
denoted @; clearly — (wre — wer)| = To 
facilitate comparison of the compression distributions 
estimated in this way they are plotted in Fig. 7(a) in 
the form @/@rg against (x — vcer)/(Xre — Xcr). Re- 
sults taken from the theoretical solution for a double- 
wedge‘ are also shown on that Figure, and it is seen 
that these are in good agreement with the mean of the 
measured values. The considerable scatter of the 
experimental values bears no consistent relation to the 
parameter x, and may be due to inaccuracies in the 
analysis of the experimental results. It is therefore 
suggested that the double-wedge distribution be used 
in the calculation of pressure distributions, and, for 
convenience, this curve is repeated in Fig. 7(b). 

The method for the prediction of sonic-range pres- 
sure distributions based on the foregoing analysis is 
discussed in its application to a particular airfoil in 
Section (4). 

The following points concerning the experimental 
results used in the analysis should be mentioned. 
The completeness of the Mach Nuinber freeze within the 
sonic range is uncertain, and, to avoid discrepancies 
from this source, the results used were obtained at a 
nominal free-stream Mach Number of unity, with the 
exception of those for the 10 per cent RAE 102 airfoil. 
This was tested in a solid-walled tunnel which has been 
found to give a complete Mach Number freeze in the 
sonic range, so that results obtained at a high subsonic 
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Mach Number with the surface shock at the trailing 
edge should not differ from those at 4, = 1. The in- 
fluence of wind-tunnel wall configuration on sonic-range 
pressure distributions is shown to be small in references 
7 and 8 for the wall shapes considered here. Because 
of this, it appears consistent to compare results for the 
6 per cent RAE 104, 4 per cent RAE 104, and NACA 
0010 airfoils, which were tested in slotted walls of 
1 30th open area, with those for the 4 per cent Bicon- 
vex and NPL 491 airfoils, which were measured with 
1 11th open area walls. These differences in test condi- 
tions do not appear to be important in the analysis. 
However, certain NACA measurements at transonic 
speeds do not conform with the relations established 
in Figs. 5 and 6. It is thought that this may partly be 
due to the low aspect ratio of the ‘‘two-dimensional”’ 
models used, and to other differences in test conditions. 


Insofar as the theory of this paper depends on a 
relation between airfoil geometry and local pressure at a 
point, it is comparable with the well known Karman 
transonic similarity law.® This states that, for affinely 
related airfoils in a sonic stream, the local pressure at a 
point is proportional to r”*, where 7 is the thickness/ 
chord ratio. The dependence of local pressure on 7 
implied by the curve of Fig. 6 cannot be easily investi- 
gated because x is not a simple function of r. How- 
ever, it is proposed to study the relation to the trans- 
onic similarity law by testing airfoils of the same basic 
section (RAE 104), but with values of 7 of 0.04, 0.06, 
and 0.10. The results for the 4, 6, and 10 per cent 
RAE 104 section given in Fig. 6 were obtained at dii- 
ferent incidences. 


(4) APPLICATION OF METHOD 


In this Section the calculation of a sonic-range pres- 
sure distribution for a particular airfoil (NPL 491) is 
described. The airfoil was designed to have a small 
trailing-edge angle, so that shock-induced boundary- 
layer separation would be absent for a range of inci- 
dence at all speeds. It is wedge shaped over the rear 
60 per cent of its chord. This is of interest since con- 
sideration of the changes in airfoil surface direction 
only would result in the prediction of constant pres- 
sure over the wedge-tail, whereas the measured dis- 
tributions show a rise in pressure towards the trailing 
edge. An analysis of the tests on this airfoil is given in 
reference 10. 


The curve of airfoil surface slope @ against x ¢ is 
shown in Fig. 8. From this the trailing-edge flow de- 
flection and the value of x can be obtained for a given 
incidence; values calculated for a = 1° are given with 
Table 1, where the calculation of the sonic-range pres- 
sure distribution is set out. Column 1 gives the chord- 
wise stations at which the pressure is to be calculated, 
starting at the crest corresponding to a = 1°, Column 2 
gives these positions in terms of XY = (x — xer)) (Xre — 
Xcr), and, in Column 3, the corresponding values of the 
surface slope measured from the crest are given. The 
values of x and ér¢ are then used to obtain wer and wre 


| 
e° 
| 
| 
| 
| 
2 ‘= 
| 
| 
| | 
4 
02 0-4 0-8 0 


z/e 


Fic. 8. Surface slope distribution for the NPL 491 airfoil. 


from Figs. 5 and 6, respectively, and hence @rg. This 
is the total compression from crest to trailing edge. 
The values of @ at the positions of Column 2 are calcu- 
lated with the use of Fig. 7(b); Column 4 lists these 
values. The differences between the values of @cr — 0 
and @ are given in Column 5; these represent the effec- 
tive increments in flow deflection to be added to wer 
to give Column 6, which lists the equivalent simple- 
wave flow deflection angles when the compression ef- 
fect has been taken into account. Column 7 gives the 
corresponding pressure distribution. 

In Fig. 9(a), the predicted and measured pressure 
distributions for an incidence of one degree are shown, 
together with that obtained on the assumption of 
Prandtl-Meyer flow from the crest. When the allow- 
ance for the reflected Mach waves is included, good 
agreement with experiment is achieved; this is fur- 


TABLE 1 


Calculation of the Sonic-Range Pressure Distribution for NPL 
491 Airfoil at a = 1° 


For this airfoil (y/c)y=0.02- = 0.0096, and from Fig. 8 it is 
found that 


(xcr/Cla=1° = 0.151 


and (x/c)gp=3° — (x/cig=5° = 0.036 
hence, x = 0.267; bre = 2.51° 
Figs. 5 and 6 give wre = 3.57° and wer = 3.90°. The total 


equivalent compression angle from crest to trailing edge is then 
given by 


@re = bre — (wre — wer) = 2.51° — (3.57° — 3.90°) = 2.84° 


1 2 3 4 5 6 7 
Col. 3 wCR + 

x/c X @e° @° Col. 4 Col. 5 p/ He 
0.151 0 0 0 0 3.90 0.405 
0.2 0.06 1.00 0.43 0.57 4.47 0.393 
0.3 0.18 2.07 0.97 1.10 5.00 0.382 
0.4 0.29 2.40 1.38 1.02 4.92 0.384 
0.5 0.41 2.50 1.70 0.80 4.70 0.388 
0.6 0.53 2.51 2.03 0.48 4.38 0.395 
0.7 0.65 2.51 2.31 0.20 4.10 0.401 
0.8 0.76 2.51 2.51 0 3.90 0.405 
0.9 0.88 2.51 2.75 —0.24 3.66 0.410 
1.0 2.5 2.84 —0 3.57 0.412 
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Fig 9 (a) 
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Fic. 9. Predicted and measured pressure distributions for the 
NPL 491 airfoil. 


ther illustrated in Fig. 9(b) where results for a = 0° and 
a = 3° are given. 

Although for this thin airfoil the surface shock wave is 
observed to be slightly forward of the trailing edge 
when .\/) = 1, whereas in the estimation of the pressure 
distribution it was assumed to be at the trailing edge, 
the results are comparable because of the persistance 
of the Mach Number freeze ahead of a shock wave when 
My is below the sonic range (see Fig. 1). The results 
for the NPL 491 airfoil shown on Fig. 5 were obtained 
on this premise by extrapolating the measured pres- 
sure distribution to the trailing edge. Values of wre 
so obtained differ little from those measured with the 
shock wave at the trailing edge when M, > 1. 


(5) CONCLUDING REMARKS 


An analysis of the Mach-wave pattern of the sonic 
flow about an airfoil surface has led to an understand- 
ing of the factors governing the flow, and it has been 
possible to represent these in a simple form so that a 
quantitative analysis of experimental results could be 
made. The empirical relations established by this 
analysis can be written as follows: 


wre = const. + m bre (1) 
wer = const. + n x (2) 

which together give the compression angle, 
are = const. + mx + (1 — m) re (3) 


Of these three relations only two appear to be inde- 


pendent, and it is of interest to consider which are thie 
fundamental relations. 

Earlier it is noted that the trailing-edge relation (1) 
appears to satisfy approximately the downstream free- 
stream conditions, but it is difficult to see the physical 
mechanism by which these can govern the supersonic 
flow upstream of the trailing edge, particularly for 
free-stream Mach Numbers above unity when the 
supersonic flow extends to infinity in a lateral direc- 
tion.* It may be that the development of the whole 
flow to satisfy the trailing-edge conditions arises from 
the instability of alternative flows. The physical mech- 
anism whereby airfoil leading-edge geometry plays a 
large part in determining the overall flow pattern is 
easily understood, and relations (2) and (3) demon- 
strate how quantitative results can be obtained from a 
simple representation of this geometry. It is perhaps 
significant that, in the theoretical studies of sonic flow, 
it has not been necessary to prescribe trailing-edge flow 
conditions nor to ensure that the downstream free-steam 
conditions are satisfied. In these studies, the trans- 
onic flow up to the limiting Mach line is treated as an 
entity, and then the supersonic flow to the trailing edge 
is obtained by the method of characteristics. Even so, 
values of wry derived from the solutions for a double- 
wedge and flat plate at incidence are in fair agreement 
with the experimental results of Fig. 5. 

The above considerations suggest that the leading- 
edge geometry of an airfoil is fundamental in deter- 
mining the overall flow in a sonic stream and that the 
correlation of wrz with trailing-edge flow deflection 
is a consequence of this. However, there must clearly 
be some means whereby the downstream free-stream 
boundary conditions are satisfied, for the pressure rise 
through an oblique shock to give the required flow 
deflection at the trailing edge may not be sufficient to 
attain the appropriate downstream pressure. This 
difficulty can be overcome by supposing that the re- 
sidual pressure rise takes place under the influence of 
the Mach waves reflected at the sonic line to down- 
stream of the airfoil trailing edge. The interaction of 
such wave systems from the two surfaces of an asym- 
metric airfoil configuration may result in small changes 
in flow direction downstream of the trailing edge, which 
may be necessary to eliminate a downwash at the trail- 
ing edge. However, analysis of measured results indi- 
cates that such a downwash will be very small. Since, 
to satisfy the free-stream conditions downstream, it is 
necessary that the total effect of expansions and com- 
pressions behind the sonic line is zero, it appears that 
the free-stream conditions must to some extent deter- 
mine the sonic-line/limiting-Mach-line configuration. 
It seems reasonable to conclude that (1) the free-stream 
conditions are dominant in determining the sonic-line 
configuration far above (and below) the airfoil, so that 
the total reflected expansion—i.e., the limiting Mach 


* At the time reference 6 was written, data were available 
only for free-stream Mach Numbers below unity, and the physi- 
cal nature of the flow at low supersonic speeds was not consid- 
ered. 
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line position—is such as to ensure that the correct 
downstream velocity is attained, but (2) that the lead- 
ing edge of the airfoil dominates the flow near the 
surface. It may be that the influence of the free-stream 
boundary conditions extends sufficiently towards the 
surface to affect the compression distribution over the 
surface and accounts for the small dependence of the 
compressive effect on 67, shown by Eq. (3). The as- 
sociated correlation between trailing-edge pressure and 
flow deflection presumably arises from the general 
similarity* of the airfoil profiles considered in reference 
6, and the fact that for these the bulk of the reflected 
Mach waves return to the plane of the airfoil ahead 
of the trailing edge, so that the pressure rise down- 
stream of the trailing edge is always small. The 
approximate correlation between wrg and 67, would 
then follow from the presence of an oblique shock at the 
trailing edge related to trailing-edge angle and to only 
slightly different downstream velocities. In a real 
flow, the change in displacement thickness of the wake 
downstream of the trailing edge will also have some ef- 
fect on the downstream pressure variations. De- 
flections of the streamlines due to reduction in the 
displacement thickness of the wake must in some way 
(at present unknown) be related to the ‘“‘total effect of 
expansions and compressions.” 

These tentative conclusions do not invalidate the 
analysis of the preceding Sections of this paper. With- 
out the mean curve of Fig. 5 it would be difficult to 
estimate empirically the influence of 67g on @rg. In- 
deed, because of the tentative nature of the above 
discussion, it is thought better to base the method 
for estimating sonic-range pressure distributions directly 
on the correlation of experimental results shown in 
Figs. 5 and 6, rather than to use the deduced correla- 
tion between @rgz and x and bry. Of course, this cor- 
relation is implied by use of Figs. 5 and 6, but the value 
of the method as such is not dependent on the interpre- 
tation of the flow mechanism leading to the results. 
However, it must be borne in mind that the scheme is 
based on a somewhat restricted selection of experi- 
mental results, and that certain NACA transonic tests 


* In the sense that the majority of the sections considered had 
blunt leading edges and conventional maximum thickness posi- 
tions. 


do not accord with these. Further experimental work 
is being undertaken to study this difference and also 
to adapt the method to consideration of cambered 
airfoils. 

Knowledge of the sonic-range pressure distribution 
can be a useful guide to the onset of shock-induced 
boundary-layer separation effects on an airfoil. Future 
work will be directed towards determining the pressure 
distribution on an airfoil with the surface shock up- 
stream of the trailing edge. f 
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On the Base Pressure Resulting From the 
Interaction of a Supersonic External 
Stream With a Sonic or Subsonic Jet 


SUMMARY 


It is shown that the two-dimensional base pressure problems 
relating to base bleed into the wake of blunt-trailing-edge air- 
foils, or the interaction between an external supersonic or sonic 
slipstream with a sonic or subsonic jet stream of a jet engine, can 
be calculated by theoretical considerations. Constant-pressure, 
isoenergetic, turbulent mixing between the streams and the stag- 
nant fluid in the wake is assumed. The theoretical calculations 
are in good agreement with the experimental results. 


SYMBOLS 
C = U/tmar = {1 + [2/(R — 1)y2}} Crocco Number 
k = Cp/Cp, ratio of specific heat 
Gs = mass rate per unit width flowing between jet-boundary 
streamline and streamline s 
H = step height 
M = Mach Number 
P = absolute pressure 
T = absolute temperature 
¥ u = velocity component in x direction 
se x,y = coordinates in the intrinsic coordinate system 


X, Y = coordinates in the reference coordinate system 
= o(y/x), dimensionless coordinate 


n 
0 = streamline angle 
p = density 
o = similarity parameter for homogeneous coordinate y/x 
= dimensionless velocity 
Subscripts 
1, 2,3, 4 refers to conditions at cross sections indicated in 
Fig. 3 
a refers to conditions in free stream adjacent to the 
dissipative regions 
b refers to conditions at the base or in the wake 
A d refers to conditions along the discriminating stream- 
line 
ae i refers to conditions along the jet boundary stream- 
line 
iS 0 refers to stagnation value 
; s refers to streamline 
Superscript 


’‘ refers to the jet stream or lower stream 


INTRODUCTION 


gb TWO-DIMENSIONAL FLOW MODEL, as shown in 
Fig. 1, simulates the problem of base bleed into 
the wake of blunt-trailing-edge airfoils, or the inter- 
action between the external slipstream and the jet 
. stream of a jet engine or rocket flying at sonic or super- 
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sonic speed. Experimental investigations dealing with 
a preliminary phase of such a problem, where the stag- 
nation temperatures in both streams are nearly equal, 
were done notably by Fuller and Reid! and Wimbrow.’ 
It is intended to show to what extent the base pressures 
for problems of this type can be calculated by theoreti- 
cal considerations. Such treatment is based on con- 
cepts emerging from a theory of base pressures, as de- 
veloped by Korst ef a/., at the University of Illinois** 
(see references 3, 4, and 5). 


Originally dealing with the problem of a two-dimen- 
sional supersonic flow past a back step, Korst devised a 
simple flow model which encompasses essentially as com- 
ponents: the approaching supersonic flow expanding 
around the trailing edge, constant-pressure turbulent 
jet-mixing along the wake,f recompression at the end 
of the wake, and conservation of mass within the wake. 
Extensions of the theory included the considerations of 
base bleed (see references 3, 4, 5, 7, and 8), and the inter- 
action of two supersonic streams at the base (see refer- 
ences 9 and 10). 


Interaction between the dissipative flow region and 
the adjacent free stream is such that the latter defines 
the pressure field in the sense of boundary-layer con- 
cepts. The free stream follows a Prandtl-Meyer rela- 
tionship for the expansion around the trailing edge, and 
goes through an oblique-shock recompression at the 
end of the wake. The boundary of the corresponding 
inviscid jet (defined as the hypothetical frictionless jet 
which has the same approach Mach Number, expands 


** This portion of the work was partially supported by the 
Office of Scientific Research under the Contract AF 18(600)-392. 

7 The case with laminar mixing was considered in a similar way 
by D. R. Chapman.* 
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Fic. 1. The two-dimensional flow model. 
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through the same pressure ratio, and is influenced by 
the same geometry of physical boundaries as the actual 
viscous jet) serves as a reference system of coordinates 
and establishes the shape of the wake (inviscid flow 
model). The dissipative flow regions are superimposed 
on the jet boundaries of the inviscid flow model. A re- 
striction of the flow model to a constant-pressure dead- 
air region is often justified, and is here assumed to be 
applicable to the present problem at the outset. Thus 
the theory of two-dimensional compressible jet-mixing 
at constant pressure*!' may now be employed. 

If the thickness of the boundary layer approaching 
the trailing edge may be considered to be small com- 
pared to the significant wake dimensions, a considerable 
simplification of the problem is achieved by application 
of the “‘restricted theory’’* which allowed a quantita- 
tive solution of the base pressure problem. Again, the 
applicability of the ‘‘restricted theory” is here assumed 
for the problem under consideration, so that the solu- 
tions may be interpreted as being reached asymptoti- 
cally for high Reynolds Numbers. 

Recompression at the end of the wake utilizes the 
concept of the discriminating streamline within the 
shear flow region, which attains just sufficient energy to 
penetrate into the high-pressure region behind the ob- 
lique shock at the end of the wake by complete isen- 
tropic (although irreversible-diabatic) conversion of its 
kinetic energy. The requirement of conservation of 
mass within the wake is here to be satisfied by balancing 
the contributions of both streams to the wake within 
the concept of an “open wake.’’* 

It shall be noted that an extension to nonidentical 
stagnation temperatures in both streams principally re- 
quires the application of nonisoenergetic mixing con- 
cepts* and use of an energy balance for the wake. Such 
a treatment, which is considerably more tedious in 
numerical calculation procedures, would also allow the 
prediction of base temperatures. 


THEORETICAL FLOW MODELS 


The present analysis will be applied to a fixed base 
geometry as shown in Fig. 1, so that for a given slip- 
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stream approach Mach Number 1/, 2 1, the primary 
variable will be the pressure ratio Py’ P;. As men- 
tioned before, the stagnation temperatures in both 
streams are assumed to be the same; also, both streams 
have the same composition and identical (and constant) 
specific heat ratios (k = 1.4). 

It is now instructive to inspect a typical experimental 
result (Fig. 2) to illustrate the appearance of three dis- 
tinctly different regions—namely, region I for Py’/ P, < 1, 
region II for 1 < Po’/P,; < ~2.35, and region III for 
P,'/P,; 2 ~2.35. While Fig. 2 represents experimental 
data obtained at the small supersonic facility in the 
University of Illinois, the same general trend was re- 
ported by Fuller and Reid! and Wimbrow.” 

The different regions as shown in Fig. 2 can be asso- 
ciated with different idealized flow patterns as illus- 
trated by Figs. 3 (region III), 4 (region I), and 5 (re- 
gion II). It is immediately recognized that the prob- 
lem associated with region III (Fig. 3) is essentially a 
special case (7,’ = 1, and the presence of the lower 
wall) of the one investigated by Korst and Tripp.° 
Note that in this specific case, the slipline angle 4;, 
downstream of the oblique shocks is positive, as shown. 

Region I, on the other hand, may be studied on the 
basis of the base-bleed concept, which was also de- 
veloped before for negligible’ and finite’ bleed velocities. 
The former will be directly applicable for values of 
P.’ P, near unity, and the influence of finite bleed 
velocities is disregarded here. As shown in Fig. 4, the 
entire mass of the jet flow is added to the wake and must 
be pumped out by the external slipstream. 

Region II, associated with the flow model as shown 
in Fig. 5, is treated with the assumption that the recom- 


Inviscid flow model for region III. 
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Fic. 5. The flow model for region IT. 

pression of the jet stream reaches the static pressure 
level impressed by the external slipstream turning into 
a direction parallel to the lower wall. Only a portion of 
the jet stream which is slowed down sufficiently by the 
mixing process along the wake will be added to the wake 
while the external slipstream has to pump mass out at 
the same rate. 


THEORETICAL CALCULATION OF THE JET INTERFERENCE 
EFFECTS 


In order to calculate the jet interference effects, it is 
necessary to treat the flow components within the dif- 
ferent flow models on a quantitative basis. 

(1) The streams approaching the base are considered 
to be uniform and parallel to the wall (symmetry line 
of the jet flow). 

The external flow will then be specified by 1; 
(21), pr, To, k, while the jet flow is given by MM,’ < 
1, po’, To, k’ = k = 1.4. 

(2) Expansion at the base. The external stream will 
undergo a Prandtl-Meyer expansion in each of the flow 
regions I, II, and III, while the jet flow will follow 
such a relation within region III only. In regions I and 
II, essentially parallel outflow from the jet nozzle may 
be assumed (although this is unimportant for cases 
within region I). 

(3) The flow of the streams along the wake with mix- 
ing. As the wake is assumed to be of constant static 
pressure, the construction of the boundaries of the 
“corresponding inviscid jets’’ can be accomplished (for 
a given wake pressure) by potential flow concepts 
(method of characteristics) in the following cases: 


External slipstream: regions I, II, and III 
Jet stream: region III* 


* Boundaries of the ‘‘corresponding inviscid jet’’ for the under- 

expanded sonic stream were calculated by the method of charac- 
teristics for two-dimensional flow on the Illiac (Electronic Digi- 
tal Computer, Engineering Research Laboratory, University of 
Illinois). 
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In absence of external interference, this boundary of 
the external slipstream will be a straight line with the 
angle of 

Along the jet boundaries, isoenergetic turbulent mix- 
ing at constant pressure is assumed to take place and 
the theoretical results obtained previously on this com- 
ponent are directly applied.* '' The flow solution in 
the mixing regions is based on a highly simplified equa- 
tion of motion, but is interpreted to hold within an 
“intrinsic system of coordinates”’ (x, y) which is subse- 
quently localized with respect to the ‘‘reference system 
of coordinates” (XY, by application of the momentum 
integral (see Fig. 3). 

Within the ‘‘restricted theory” the velocity distribu- 
tion in such a case can be conveniently represented by 


¢ = (1/2) (1 + erf n) (tf 


The mass rate per unit width, 


Vs 
G, = f pudy 
vj 


flowing between the ‘‘jet boundary streamline,” y,; or 
n) = nj(Cx_) (which separates the fluid of the approach- 
ing stream from the fluid in the dead-air region) and 
the streamline s (specified by the coordinates y, or 7,) 
can be conveniently expressed in dimensionless form by 


Cute) = 
n) 
Ti(ns, Ti(n), (2) 


where 


Coq?) = f [¢dn/(1 — Cx ?¢*) | (3) 
has been obtained by numerical] calculation on the 
Illiac for 21 values of the parameter C., (see reference 
12), and is given in graphical form in Fig. 6 as J;[¢(n), 
C:,7|.. o is the empirical similarity parameter in tur- 
bulent mixing supposedly dependent upon Mach 
Number. Here, the relation suggested by Korst and 


7 Equations are written for the external stream but are equally 
applicable to the jet stream when used with prime as superscript. 
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Fic. 6. Auxiliary function J; (¢, Coq) for »mpressible isoener- 
getic turbulent jet-mixing at constant pressure (for fully developed 
profiles ). 
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A 
Theoretico! Colculotions 
— 0 =i2 2.758 
/P, = Po /P,= 0.268 —o— Experimento! Results Reported 
Obtained by the by Fuller ond Reid 
Theory of Bose Pressure (Ref.3) (Fig.13 in Ref. 1) 

' 2 3 4 5 6 7 P,/P, 
Fic. 7. Comparison of experimental results and theoretical cal- 
culations. M, = 2.3, H/H' = 2. 

Tripp? ¢ = 12 + 2.758M was used in the calculation. 


For the treatment of the jet flow in region I, mixing 
concepts are not needed, since the jet velocity is as- 
sumed to be small. 

(4) Recompression at the end of the wake. With 
the static pressure ~; produced by the oblique shock in 
the external slipstream, the escape criterion’ defines 
a discriminating streamline d with the dimensionless 
coordinate ng, so that 


Poa,/ Ps = 1 (4) 
where 
Poa,/ Ps = + [(k 1) (5) 


This relation, which will allow the evaluation of AJ, 
and, consequently, Cyg and ¢3«,* is needed for the ex- 
ternal flow (regions I, II, and III) and for the jet 
stream (regions II and III only). 

Thus, for a selected base pressure p,) = pz: = ps, the 
streamline angles of the inviscid flow model will be 
determined and the static pressure behind the recom- 
pression zone will be fixed. While this is evident for 
regions I and II, where the shock deflection angle must 
be equal to 6.,, a trial-and-error procedure is required 
for region III so that a 64, is found for which Py = P,’. 

Since the base pressure p, was so far assumed to be 
known, the principle of conservation of mass within 
the wake must now be employed to find the correct 
solution. 

(5) Conservation of mass in the wake. 
that for steady state 


Ga + Gia’ 


It is evident 


Il 


(6) 
or 

Gi’ Ga = —] (7) 
which can be expressed, with the help of the relations 
given above, by 


Gy’ Coq" sin O24 


Ge Can (1 — (H/H’) — Ty) 


=-1 (8) 


* Note that the discriminating streamline d will coincide with 
the jet boundary streamline j for closed-wake problems only. 


for region I, and 


Ga Coq (Lig — Ly) 
for regions II and III, where the term x’/x corresponds 
to the ratio of the rectified length along the mixing 
paths of the streams in the inviscid flow model (ap- 
proximated by cos 42, in case of region IT). 

Problem solutions for given height ratio /7/H’ and 
values of 1, and P,’/P, are obtained by a trial-and- 
error procedure where values of p,/); are assumed and 
the relation for mass conservation has to be satisfied 
as a check. 


THEORETICAL AND EXPERIMENTAL RESULTS 


It was first attempted to produce a theoretical solu- 
tion for a flow case investigated experimentally by 
Fuller and Reid (test on afterbody No. | in reference 1). 
The geometry of this afterbody was similar to that of 
Fig. 1 with /7/H/’ = 2, and the Mach Number of the ex- 
ternal stream approaching the base was reported as 
being 14, = 2.3. Comparison between the theoretical 
curve calculated by the method given in the previous 
sections showed good agreement within region I but 
only qualitative agreement for regions II and III (see 
Fig. 7). A careful study of the test setup used by Fuller 
and Reid revealed that the shock originating at the lead- 
ing edge of their model, as well as the expansion fan 
originating at the shoulder of the leading wedge, were 
reflected from the tunnel wall and interfered with the 
wake, so that their experimental results do not seem 
to be representative for undisturbed conditions of the 
external flow. 

In order to produce test conditions which were free 
of such a deficiency, thus allowing a better evaluation 
of the theory, an experimental program was conducted 
at the University of Illinois. The configuration of the 
flow channel used in this investigation is shown in 
Fig. 8. Air for both streams was taken from the same 
supply line (7) = 7)’) and the variation of P,’/P; was 
achieved by appropriate individual throttling of the 
valves leading to the two streams. The external flow 
(upper stream in Fig. 8) was simulated by a sharp- 
corner nozzle designed for parallel outflow which pro- 
duced an exit Mach Number of 1/, = 1.92 of reasonable 
uniformity, as checked by pressure readings in the exit 
cross section. The jet stream (lower stream in Fig. 8) 


was produced by a converging nozzle. H/H’ was 
1.94 as mentioned before. 
PR RE TAP: 
0.6 LB/SEC = Potm 
AIR SUPPLY 
4 iif 


Fic. 8. Experimental configuration (1/2-in. width). 
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Agreement between experimental data and theoretical 
calculations, as evidenced by Fig. 2, is quite satisfac- 
tory, except for the region around Py’/P,; = 1, where 
the test configuration of the internal flow will fail to 
simulate correctly the condition for undisturbed ex- 
ternal flow.* 

It is noteworthy to mention that the theory does not 
introduce any empirical quantity other than o. That 
the uncertainty in dealing with the latter is of little 
practical consequence is demonstrated by plotting into 
Fig. 2 another theoretical curve obtained for ¢ = 12, 
a value well established for J —> 0. 

From this method of calculation, it can be predicted 
that the base pressure ratio would have a relative maxi- 
mum for a pressure ratio of Py’/P, equal unity, if one 
neglects the irreversibility of the shock recompression 
The actual irreversibility of the shock of moderate 
strength will cause the relative maximum to occur for 
P'/P; less than but close to one. This result is at 
variance with the observance of Fuller and Reid! and 
Wimbrow? who locate the relative maximum of the base 
pressure ratio for values of Py’ P; below unity. 

It is also interesting to point out that the geometrical 
height ratio // ///’ has little effect on P,,/P, in region IT. 
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An Exact Solution of the Navier-Stokes 
Equations 


Herbert E. Hunter 

Graduate Student, California Institute of Technology, 
Pasadena, Calif. 

October 13, 1958 


ie THIS NOTE the exact solution of the Navier-Stokes equations 
presented by W. Chester in reference 1 is generalized to in- 
clude all values of y. The notation of reference 1 will be fol- 
lowed, except as noted. 
In rectangular coordinates the solution presented in reference 
1 may be written as 
uw, = U(1 — y?/a?) 
(1) 
fp: = pe — (x — UNS 
It can be shown that Eqs. (1) are a nontrivial solution only if 
y = 1.5. It seems likely that this solution is part of a more 
general solution which has other terms involving (y — 1.5). 
We shall assume that the more general solution is of the form: 
= [u(y, z, 0), Under this assumption, the 
continuity equation is satisfied identically and the remaining 
Navier-Stokes equations reduce to 


Op/Ox = + | — p(du/dt) 
= dp/ds = 


(Op/dt) + u(Op/dx) = — 1) + 
u(Ou/Oy)? + 


0, O], p = const 


Integration of the momentum equations gives 

p= p(x, 
= + — p (Qu/dt)} x + G (t) (3) 
After examining the details of the solution presented in reference 
1, we are led to trying a solution of the form: 


(4) 
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and to the assumption that 
Ci = (2nU*t/a?) + b) (5) 
where “; = solution from reference 1 (Eqs. 1); u’ = u"(s, t, y); 
and C,’ = (7,2). 
Combining the last of Eqs. (2) with Eqs. (3-5) and taking into 
account the fact that u’ is independent of y, we get the following 
two equations: 


— p(du’/dt) = —(4 U/a?) (y — 1.5)p 


= (4 U?/a?) (y — 1.5) w + (6) 
(y — 1) [(4 nU/a*)u’ + 
The last of these and Eq. (5) imply 
(Ou’/dz) [(02u’/ds?) + (2u/a’)] = O 
Solving this we get 
u’ = — (Us?/a*) | + (7) 


where 7 is either 0 or 1. 
Substituting the above equation into the first of Eqs. (6) and 
solving it, we get 
u’ = nU [Coz — (2?/a?)] + (QuU/a%p) (2y — 3 — (8) 
where C2 is a constant. 
Solving the last of Eqs. (6) for C,’ we get 
Ci’ = ([4(y — 1)(2y — 3 — n)p/pat) t? + 
[4(y — 1.5)/a?] + (y — 1) C2 n?} t) (9) 
Therefore the more general solution for u is 
u = — [(y2/a®) + (ns?/a®) |] + ness + 
(2Qu/a%p) (2y — 3 — n)t} (10) 
Substituting this into Eq. (3) we get 
p = po — [4(y — X 
{x — [(2y — 3 — n)pU/pa?|t? — [1 + (coan/2)2]Ut} (11) 
The case = O is clearly a two-dimensional flow. The 
case nm = 1 is best interpreted by means of cylindrical coordinates, 
in terms of which it becomes 
u = U[l — (r?/a?) + cor sin + (y — 
= po — — 1)#U/a?] X 
{x — — 2)uU/pa?|t? — [1 + (cxa/2)2] Ut} 
The solution for » = 0 represents plane flow between two in- 
sulated parallel flat plates which have a time-dependent distance 
between them. Eqs. (12) represent flow through an insulated 
expanding or contracting circular pipe. The direction of motion 
of the walls in either case is dependent on y. It is easy to show 
that the above solutions reduce to those presented in reference 1 
for y = 1.5 and y = 2, respectively. That is, the distance 
between the plates (or the radius of the pipe) remains fixed. 


(12) 
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A Note on the Skin-Friction Coefficient for 
Compressible Couette Flow 


Robert M. Inman 
Assistant in Engineering, Department of Engineering, University of 
California at Los Angeles, Calif. 


October 10, 1958 


SYMBOLS 


= skin-friction coefficient 
specific heat at constant pressure 
distance between plates 


Cy 
Cp 
h 


M = Mach Number 
Pr = Prandtl Number 


q = heat-transfer rate 

R = gas constant 

Re = Reynolds Number 

T = fluid temperature 

u = axial component of velocity 

U = velocity of moving plate 

7 = ratio of specific heats 

uw = viscosity of the fluid 

£ = dimensionless velocity, u/U 

p = density of the fluid 

7 = shear stress 

w = exponent for viscosity power law 
Subscripts 

w = conditions along stationary plate 

« = conditions along moving plate 

ry = recovery conditions 


¥ HAS BEEN SHOWN! that, for compressible Couette flow, the 
consideration of viscosity and heat-transfer effects results in 
the following equations for the temperature and velocity dis- 
tributions within the fluid, respectively: 


CAT — Tw) + (Pru?/2) = —(qw/tw)Pru (1) 


1 
U) f, {1 (Prqe/ twU)(7 1)M,*(1 = 
(Pr/2\(y¥ — 1)M,%1 — (2) 
where the variable £ = u/U has been introduced. w arises from 
the assumption that ~ 
Eq. (2) has been solved for the simple case of w = 1 and qy = 0 
in reference 1, and the solution is given as 


= 1+ (Pr/3) (7 — 1)M.? (3) 

The skin-friction coefficient, Cy; = rtu/(p./2)U?, was shown 
to be 

(C;Re/2) = 1+ (Pr/3)(y — 1)M,? (4) 


where Re = 

We want to point out that a more general expression for the 
skin-friction coefficient can be obtained when the heat-transfer 
rate, gw, is not assumed as zero. The assumption that w = | 
will again be made. 

For arbitrary heat transfer, Eq. (2) is integrated easily, to give 
(ruh/poU) = 14+ (Prqu/2Utw) (y — 1)M.2 + 

(Pr/3)(y — 1)M.? (5) 
Now putting the upper-wall conditions into Eq. (1) gives 
— To) = Pr(qu/tw)U + (PrU?/2) (6) 
or 
Pr(Qu/twV) = — To.) — (PrU*/2)|/U? (7) 
With Eq. (7), Eq. (5) becomes 
(rwh/toU) = 1+ (Tw — + 
(Pr/12)(y — 1)M,2 (8) 
where the definition 7,2 = U?/(yRT,.,) has been employed. 

The skin-friction coefficient becomes 
(CyRe/2) = 14+ (Te — To)/2T + (Pr/12)(y¥ —1)M.? (9) 
Eq. (9) permits determination of Cy, for all ratios of (7./T,, ). 
In particular, if the lower wall is insulated (g, = 0), so that it 
attains recovery temperature 7,, given by 

(T;/T) = 1+ (Pr/2)(¥ 1I)M (10) 
then Eq. (9) reduces to Eq. (4) upon the use of Eq. (10). 

Eq. (9) points out the strong dependence of the skin-friction 

coefficient on the heat-transfer rate as well as the Mach Number. 
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The Stagnation-Point Boundary Layer on a 
Rotating Hypersonic Body* 


Sinclaire M. Scala and Joseph B. Workman 


Missile and Space Vehicle Department, General Electric Company, 
Philadelphia, Pa. 


October 2, 1958 


SYMBOLS 


B; = Xi/m 

f = similarity stream function 

fn = u/ue 

hy = enthalpy of species f, including chemical 

pu/ Pwhw 

Lij = Lewis Number 

L* = thermal Lewis Number 

m; = molecular weight of ith species 

m = 2 Xim; mean molecular weight 

M; = chemical-source term, moles per unit volume, per unit time 

= rpm 

7 = temperature 

Pi = Prandtl Number 

x, y, @ = body-oriented coordinate system 

u,v, w = velocity components 

Vo = flight speed 

r = radius from axis to surface of body 

RB = nose radius 

& = angular velocity of body, radians per sec 

X; = mole fraction of ith species 

u = viscosity 

p = density 

” = similarity variable 

u = (w/wr) dimensionless velocity 
Subscripts 

crit. = critica! 

e = outer edge of boundary layer 

u = wall 

” = denotes differentiation with respect to 7 


gosta NOTES have reported on catalytic effects' and surface 
combustion? at hypersonic speeds. To date, however, the 
effect of the rotation of a body about its axis has been neglected. 
The purpose of this note is to present some results of an investi- 
gation into the effect of spin on the hypersonic laminar boundary 
layer. 

The equations pertinent to the study of the laminar boundary 
layer on an axially-symmetric rotating body traveling at hyper- 
sonic speed include the conservation of mass, momentum, energy, 
and chemical species. For steady flow, these may be written: 
Continuity 

(0/dx) (pru) + (0/d¥) (prv) = 0 (1) 


Momentum x 


p|u(Ou/dx) + v(Ou/dy)] = 
—[(dp/dx) — p(w?/r) (dr/dx)] + (0/dy) [u(Ou/dy)] (2) 


Momentum 

p(u} [O(rw)] /dx} + vf [O(rw)/dy]}) = (0/dy) (uf [O(rw)] /dy¥}) 
(3) 

Energy 


p|u(Oh/dx) + v(Oh/dy)] = u(dp/dx) + 
ul(Ou/dy)? + (Ow/dy)?] + (0/dy) [K(OT/dOy)] — 


(0/dy){ he(umi/Pr) — 
(L;7/m;T) (OT /dy)|} (4) 
Chemical Species i 
p(u(0B;/dx) + v(0B;/dy)] + 


(L;7/m;T) (0T/dy)|} = Mi (5) 


* This analysis is based upon work performed under the auspices of the 
U. S. Air Force, B.M.D. Contract No. AF 04(645)-24. 
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and represent a generalization of those given, for example, by 
Chu and Tifford.* 

One observes that the effects of body rotation do not appear 
explicitly in either the global continuity equation or in the con- 
servation of species. Further, at a stagnation point the viscous- 
dissipation term u[(0u/dy)? + (Ow/dy)?] vanishes. Thus, upon 
introducing a similarity transformation and restricting the dis- 
cussion to the stagnation point, the effects of body rotation appear 
only in the momentum equations: 


(lfan)n + + 1/21 — + = 0 
(6) 


+ ft, — = 0 (7) 


for the x and ¢ directions, respectively. 

In order to obtain an indication of the combination of physical 
and thermodynamic conditions under which the effects of body 
rotation can be neglected, one may take as an approximate cri- 
terion, the condition 


(1/2) 0.01 (8) 


It is not difficult to show that a first-order relation between nose 
radius Rg, flight speed V,,, and stagnation-point velocity gradient 
du,/dx, is given approximately by 


(Rp du,/dx)"*?| |5V, X 1074 (9) 


to within plus or minus 25 per cent, at flight speeds above Mach 
10, regardless of altitude where all quantities in Eq. (9) are in ft. 
per sec. units. 

Thus, since # < 1.0, Eq. (8) becomes 


Nerit. S 3.4 K 10-5 (r.p.m. sec.?/ft.) (10) 


It can be immediately concluded that the effect of spin may be 
completely neglected in the analysis of hypersonic stagnation- 
point heat transfer for small vehicles, since, for example, a body 
having a 1-ft. nose radius when traveling at 25,000 ft. per sec., 
must rotate at speeds above 14,000 r.p.m. before the structure 
of the boundary layer is affected. 
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On Simititude in Rarefied-Gas 
Aerothermodynamics 


F. Marcel Devienne 
Director, Mediterranean Laboratory of Thermodynamic 
Research, Nice, France 


October 8, 1958 


i RAREFIED-GAS AEROTHERMODYNAMICS, similitude is indi- 
cated by several numbers without dimensions, one of which is 
the Knudsen Number, A, equal to the ratio of the mean free path 
of the gas molecules to a characteristic dimension of the body. 

When we operate in the free-molecular regime—i.e., for K 
values which are theoretically superior to 10*— the accommoda- 
tion coefficient intervenes in the formulas. Thus, the energy 
received by the unit of surface of a flat plate moving perpendic- 
ularly to its plane is given by the formula: 


* The researches carried out by Stalder and his collaborators and by our- 
selves have shown that the free-molecular flow regime takes place for Knud- 
sen Number values inferior to 10. In the case of a cylinder, Stalder finds a 
Knudsen Number equal to 2; for the flat piate we have found a Knudsen 


Number equal to 3.'» 
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W = (apv,8/2)[(s/2)®(s) + + 


where a is the accommodation coefficient of the molecules in 
relation to the surface of the plate, p—the specific mass of gas at 
the pressure considered, v,—the most probable speed of the mole- 
cules, and s—the ratio of the speed of displacement U to v,.4(s) 
represents the error function with upper limit s.2~* 

The accommodation coefficient depends on numerous factors 
and, in particular, on the pressure. The influence of the pressure 
moreover, is, very complex; in particular, in the case of pres- 
sures less than a few microns of mercury it depends on the out- 
gassing, which, in turn, depends on the various physical factors 
to which the surface is submitted and, in particular, on the pres- 
sure. We have also brought into evidence the influence of the 
displacement speed of the body in a rarefied gas.? 

It is therefore impossible to apply the results obtained in the 
free-molecular regime at different pressures in an absolutely 
certain manner. Thus, if we operate on a wire having a very 
small diameter under a pressure of the order of 50 microns, it is 
impossible to apply these results in the case of a cylinder having a 
diameter 500 times larger in the case where the pressure is only 
1/10 of a micron. The accommodation coefficient of gas mole- 
cules can vary between these two pressures in considerable pro- 
portions—for example, 10 to 1. 

This shows the indispensable precaution which should be taken 
with a view to generalizing the results of rarefied-gas aerother- 
modynamics. 
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Shock-Tube Testing Time 


G. F. Anderson 
Brown University, Providence, R.I. 
October 10, 1958 


HE FLOW Of a real gas initiated by a strong shock in a constant- 
gyi shock tube is complicated by the marked attenu- 
ation of the shock wave. 

In a theoretical investigation of attenuation effects—it was 
hoped that the investigation would provide some information on 
the state of the gas at various locations in the test slug—the 
conservation of mass equation led to an explanation of the differ- 
ence between the ideal theoretical test time and the experi- 
mentally obtained test time. 

The result of this approach shows that for a given diameter 
shock tube, there is a length of tube at which the testing time 
reaches a maximum, and beyond this length the testing time de- 
creases. The difference between the actual test time and the 
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theoretical test time is explained by the loss of test gas to the 
boundary layer behind the interface rather than to mixing of the 
driver gas and test gas at the interface. This concept is impor- 
tant when one uses the procedure of reference 1 to determine 
the state of the gas at various points in the test slug. 

The time rate of change of the mass of gas contained in the 
test slug, Fig. 1, is given by 


bi 
(d/dt)m, = ar*p(t, — U2) — 2rp(r — y) (we — u)dy 
0 


If it is assumed that the average density fz in the test slug does 
not vary strongly with time, it may be stated that the test time 
will reach a maximum when 


2rp(r — y) — u)dy = — = Us 
0 


which reduces to 


pilts/ potty = (2/r) (8;* — for r 
= [1 — 
( 


-f [1 — (p/p2)|dy 
0 


= 1/(M, — 1) 
And, therefore, the criterion for 7,,,, is given by 


1/2(8;* — 8°) = M, —1 


where 


For strong shocks, 


Assuming a turbulent boundary layer, let 


u/us = 
= (T./T2) (1 + agl/7 — bg?/7) 


Then 
(8:* — = 7(T2/T.) (Is — 


2"dz 
0 1 + etc =— bs? 


a =(T,/Te) — 1 ((T2/Tw) — > 1 
b = [(T,/T2) — 1) (T2/Tw) = (0)m, >1 
Therefore, the criterion for 7,,,, becomes 
— = M, 1 
For 6; not small compared to 7, the above relationship becomes 
(7/8:)/{14(T2/Tw) — Iz) — (6:/r) (is — |} = 1 
The test time at any tube station may be calculated if the as- 
sumption is made that 1, = constant. This is not a realistic 


assumption but will show the order of magnitude of the boundary- 
layer effect on testing time. 


where 


Feorr/ Ttheo — 1)/x,] x 


Xs 
f 7(T2/Tw) (Te — Iz) (26;/r)dx,, 6; >r 
0 


or for 6; not small compared to r, 


Teorr/ Ttheo = 1)/x,] x 


Xs 
0 


which may be reduced to 
teorr/Ttheo = 1 — (Kix,4/> — Kox,8/5) 


where K, and A> are functions of the shock Mach Number, the 
tube radius, and the initial tube pressure (P,) before the shock. 
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FIG 2 
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These parameters are plotted in Fig. 2 for two tube sizes and 
P, = 10cm. Hg. 

An example has been calculated for 17, = 6.0 and P; = 10 cm. 
for the two tube sizes and the results are plotted in Fig. 3. The 
50 ft. station, measured from the diaphragm, is indicated on the 
result for the 4-in. tube. It is seen that little increase in test time 
can be expected beyond this station. This agrees well with ex- 
periments with a 4-in. tube. Furthermore, the corrected test 
time is approximately one half of the ideal theoretical test time 
which agrees with the experiment. The same holds true for the 
15-ft. station in a 1'/2 in, diameter shock tube. 
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Correlation of Base Pressure and Wake 
Structure of Sharp- and Blunt-Nose 
Cones With Reynolds Number Based on 
Boundary-Layer Momentum Thickness 


R. Lehnert and V. L. Schermerhorn 
Naval Ordnance Laboratory, White Oak, Silver Spring, Md. 
October 10, 1958 


T HAS BEEN ESTABLISHED in the past that there is a certain 
relationship between base pressure and boundary-layer be- 
havior.'~* The base-pressure and wake-flow conditions were 
found to be dependent upon the local flow characteristics at the 
surface of supersonic vehicles directly upstream of the base or of 
the region of wake-flow separation. 

In order to use existing data on cones and other shapes to pre- 
dict wake angle and base pressure on blunt bodies, an attempt was 
made recently at the Naval Ordnance Laboratory to establish 
a unique relationship between given local flow conditions at the 
downstream end of sharp and blunt cones at supersonic speeds 
and the corresponding wake-flow conditions with zero heat trans- 
fer. The models chosen for this investigation were a 10° semi- 
vertex angle, sharp-nose cone and a 10° semivertex angle, spheri- 
cal-nose cone with a 0.5 bluntness ratio. The bluntness ratio is 
defined as the ratio of the nose radius to the base radius. Both 
cones had a base diameter of 6 in. They were sting mounted 
(1.5-in. sting diameter) in the test section of the Naval Ordnance 
Laboratory’s 40 & 40 cm. Aeroballistics Tunnel No. 2. The 
sharp cone with surface finishes of 15 and 125 microinches (both 
designated as ‘‘smooth’’) was tested at Mach Numbers of 2.74 
and 2.76, respectively, while the blunt cone was investigated at a 
Mach Number of 4.74 with different surface conditions—i.e., 
“smooth” (125 microinch finish) and “‘rough.’’ The rough sur- 
face was achieved by application of evenly distributed carbo- 
rundum powder (No. 24 grain) over the entire wetted surface 
area. 

Surface pressure-distribution and base-pressure measurements 
were carried out and total pressure-wake surveys and boundary- 
layer profiles near the base were taken at zero angle of attack. 
The variation of the boundary-layer profiles immediately up- 
stream of the cone bases was accomplished by varying the fol- 
lowing parameters: pressure gradient (sharp cone, blunt cone). 
surface conditions (smooth, rough), and Reynolds Number 
(wind-tunnel supply pressure and temperature). The local 
Mach Number on the surface at the downstream end of the cones 
was the same for both models due to the proper choice of the 
ambient Mach Number. 

Simultaneously with the experimental investigation, a calcula- 
tion of the boundary-layer development along the cone surfaces 
was performed by a stepwise solution of the boundary-layer 
momentum equation, using the procedures from references.‘ 
Prevailing pressure gradients and surface roughness, as well as 
the systematic variation of supply conditions, were taken into 
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account for each individual case tested in the wind tunnel. The 
calculated Reynolds Numbers based on boundary-layer momen- 
tum thickness at the end of the cones were in good agreement 
with the ones obtained by evaluation of the momentum thickness 
from boundary-layer surveys. 

The measured base pressures are presented as p,/p. = f( Reg) in 
Fig. 1 where p, = base pressure, p. = static pressure on the sur- 
face at the downstream end of the cone, and Reg = Reynolds 
number based on boundary-layer momentum thickness near the 
end of the cone. Fig. 2 shows the corresponding wake angles as a 
function of Reg obtained from wake surveys, schlieren photograph 
evaluations, and Prandtl-Meyer expansion computations. 

From these results, and in view of the model geometry and the 
parameters considered, it can be concluded that the base pressure 
ratio p,/p, and the corresponding wake angle are unique functions 
of Reg, independent of the history of the boundary-layer whenever 
the local Mach Number on the cone surface near the base is the 
same. A more detailed description of this investigation is given 
in reference 9. 
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The Role of Unit Reynolds Number in 
Boundary-Layer Transition 


Jack D. Whitfield and J. Leith Potter 
Gas Dynamics Facility, ARO, Inc., Tullahoma, Tenn. 
October 13, 1958 


XPERIMENTS HAVE DEMONSTRATED that variation of the 
free-stream unit Reynolds Number, U/v, often produces 


changes in the usual forms of Reynolds Number used to describe 
the mean location of boundary-layer transition. Recognition of 
the possible influence of the unit Reynolds Number is important 
in the examination of data which are sensitive to the location of 
transition. There is doubt as to the status of the L’/v effect 
ie., whether it is associated with boundary-layer stability or 
entirely related to such factors as surface roughness or leading- 
edge bluntness. This note shows how a possible basis of the 
unit Reynolds Number effect is revealed by examination of the 
stability of laminar boundary layers. Subsonic flow is assumed, 
since it is desired to compare the results with experimental 
data; and data from supersonic flow experiments may include 
leading-edge shock-wave phenomena which confuse the search 
for pure effects.! 

A measure of instability, which, we assume, is a measure of 
closeness to transition, may be formed by the ratio of virtual 
shear stress to laminar shear stress, 7*. Thus, for a semi- 
infinite flat plate at zero angle of attack in a uniform sub- 
sonic stream, we define 

r* = 7,/71 = pu,v;/p(du/dy) 
or, from Liepmann (reference 2), 
= —(2/C,,)Kb(u’/ U)o2A? (1) 


where C,, = laminar friction coefficient, AK = correlation coeff- 
cient = u,7,/u'v’, b = v'/u’, (u'/U)o? = energy of frequency most 
amplified, and A = amplification ratio. Now, if the boundary 
layer is equally stable or unstable at two different conditions, 
(1) and (2), we may write 


m1*/72* = 1 (2) 


Following Tetervin’ and using an expression for amplification 
ratio given by Schlichting,’ it can be shown® that Eqs. (1) and 
(2) result in the equation 
[( Res*)1/(Res*)o]? — 1 + 
({In[(u’/U)o, 1/(u'/ U)o, / [S(Res*)s?]) 
{ /2S(Res*)2?} = 0 (3) 


where Res* is the boundary-layer Reynolds Number based on 
displacement thickness and S is an amplification parameter. 
Eq. (3) expresses a relation between initial energy levels at the 
beginning of amplification of the frequency of maximum ampli- 
fication and the boundary-layer Reynolds Number. One solu- 
tion to Eq. (3) is obvious when (u’/U)o,1 = (u’/U)o, 2 then 
(Res*), = (Res*)e. A more general case exists when the initial 
energy levels are not equal, say (u’/U)o,1 > (u'/U)o, ». ~Then it 
follows from Eq. (3) that 


(Re s*):2/( Res*)o? + [1/2 S(Res*)s?] In [( Res*):/Res*)o] <1 (4) 
A necessary condition for Eq. (4) is 
(Res*): < (Res*)s 


and it follows, for the conditions initially specified, that the 
relationship between the transition Reynolds Numbers, based 
on wetted length to transition, is 


(Re,) < (Re, )2 (5) 


There is experimental evidence that the critical frequency 
increases with U/yvy and that u’/U commonly decreases with 
frequency in the usual range of frequencies where maximum 
amplification occurs. 7 Therefore, we conclude that it is 
possible and even most likely that (u’/U)o,1 > (u’/U)o, 2 when 
(U/v), < (U/v)o. Thus, if the previously undefined conditions 
(1) and (2) refer to (U/v), and (U/v)s2, Eq. (5) results. However, 
it must be kept in mind that there are other factors influencing 
boundary-layer transition, and any one of these may dominate 
the process. For example, with increasing values of U/v bound- 
ary-layer thickness becomes small and roughness is more im- 
portant. If roughness dominates, then (Re;): > (Rez) when 
(U/v), < (U/v)s. 

The effect of U/v appears to depend on the initial energy 
level of the critical frequencies of oscillation which will receive 
maximum amplification. It seems academic whether we regard 
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U'/v as a fundamental factor or one only indirectly effective 
through its connection with frequency and energy levels of 
disturbances in laminar boundary layers. The significant point 
is the apparent probability that Re; generally will not remain 
constant with varying U’/v even when the influences of all other 
factors are invariant. 

The authors have received a report® including systematic, 
subsonic data showing variations of Re, with U’/yv in support of 
the writers’ earlier work® which is the basis of this note. 
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The Heating of Slabs With Arbitrary Heat 
Inputs 


Theodore R. Goodman 

Senior Aerodynamics Engineer, Allied Research Associates, Inc., 
Boston, Mass. 

October 14, 1958 


_ THERMAL HEAT FLUX at the surface of airplane and missile 
structures may become very large at high speeds, necessitat- 
ing the inclusion of thermal stresses in the structural analysis. 
The surface temperatures may, in fact, reach such high values 
that it is necessary to account for surface radiation. Further- 
more, even if only aerodynamic heating is accounted for, by as- 
suming the heat flux to be proportional to the difference between 
the wall temperature and an ambient temperature, nevertheless, 
the heat-transfer coefficient and ambient temperature will both, 
as a rule, be time dependent. The transient temperature dis- 
tribution in the skin or near the surface can be determined by 
solving the linear heat conduction equation. However, if surface 
radiation is important, the boundary conditions are nonlinear, 
and if the heat-transfer coefficient and ambient temperature are 
time dependent, no general solution is known. The entire prob- 
lem becomes even more cumbersome if the thermal properties of 
the skin, namely heat capacity, and thermal conductivity (and 
even density) must be considered temperature dependent, for 
then the differential equation itself becomes nonlinear. 

The standard method for solving the practical problems of 
heat conduction for which no analytical solution is possible is to 
integrate the heat conduction equation numerically, using a finite 
difference procedure. Recently, however, the author developed 
an approximate analytical technique for solving unsteady heat 
conduction problems.! This technique, known as the “heat 
balance integral,’ reduces the problem to an ordinary differential 
equation which, for many cases of interest, can be solved in closed 
analytical form. In reference 1 the emphasis was placed on prob- 
lems for which there was a change of phase. For aeronautical 
application, the heat balance integral method can be applied to 
heat-transfer problems in the skin when the heat flux at the sur- 
face is an arbitrary function of temperature and time. Thus, 


all the boundary conditions discussed in the first paragraph can 
be handled within the framework of a uniform theory. It will be 
assumed that the thermal properties are not temperature de- 
pendent, and that the solid slab is semi-infinite. Both of these 
restrictions can be removed, and some problems relaxing one or 
both of these restrictions have been set up and solved in reference 
2. The present note is, in fact, a condensation of reference 2. 

It is worthy of note that when the method of the heat balance 
integral is applied to problems whose exact solution is known, 
the error which is introduced by the essentially approximate 
nature of the method is, at the most, a few per cent. 


THE HEAT BALANCE INTEGRAL TECHNIQUE 


The heat conduction equation for the temperature, u, with 

constant thermal properties is 

x(O02u/Ox?2) = Ou/Odt (1) 
where « is the diffusivity; x and ¢ are the space coordinate and 
time, respectively. 

Assume a semi-infinite slab extending over positive x. Initially 
the temperature is zero. The surface condition is 

(Ou /dx)(0, t) = —f[u(O, t), (2) 

The method of the heat balance integral consists in assuming 
a polynomial temperature profile between the surface, x = 0, and 
a thermal layer, 6 (t), consistent with the boundary conditions. 
At x = 6 the boundary conditions are: 

u[d(t), t] = 0 (3) 
(du /dx)[6(t), t] = 0 (4) 

By differentiating Eq. (3) with respect to time and applying 

Eqs. (1) and (4), the following conditions can be derived: 
(07u/Ox?)(6, = O (5) 

Eqs. (2)-(5) constitute four conditions which the temperature 
profile must satisfy. Assuming, therefore, a cubic profile, it must 
take the form, 

u = (2/6*)(6 — x)8 (6) 
where 6 = 32/[f(z, (7) 
and z is surface temperature. 

Eq. (7) expresses a relationship between the surface tempera- 
ture and the thermal layer, 6. There is, therefore, still one un- 
known function of time in Eq. (6), and this will be obtained from 
the heat balance integral. 

Integrating Eq. (1) from x = 0 tox = 6 there is obtained 


x[(Ou/dx)(5, t) — (Ou/dx)(0, t)] = d9/dt (8) 

where udx (9) 
0 


The left-hand side of Eq. (8) is immediately determined from 
Eqs. (2) and (4). The quantity, @, may be obtained in terms of 6 
or z by substituting Eq. (6) into Eq. (9) and integrating. There 
is then finally obtained the following differential equation for 2: 


(4/3)«f(s, t) = (d/dt){s?/[f(z, t))} (10) 


The proper initial condition is z = 0 whent = 0. If the func- 
tion, f(z, ¢), is dependent on both z and ¢, this equation must be 
integrated numerically. Nevertheless, this constitutes a tre- 
mendous saving in effort when compared to the finite difference 
techniques which would ordinarily have to be performed on the 
original heat conduction equation. There are two cases for which 
Eq. (10) can be integrated directly, and these are when f depends 
solely on z, and when f depends solely on ¢. 

Case (a)—If f is independent of f, the solution is readily ob- 


tained, 
(4/3)xt = — (1) 


Case (b)—If f is independent of z, the problem is linear and can 
be solved exactly by more conventional methods. Nevertheless, 
the simplicity of the solution given here warrants attention. The 
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result is 
t 1/2 


APPLICATIONS 


If f is constant Eq. (12) reduces to 
= V4/3fV (13) 


The exact solution appears in Carslaw and Jaeger,’ and the re- 
sult is 


at (14) 


The numerical factor V 4/3 = 1.15, whereas V4/x = 1.13. 
Hence, the two solutions agree to within 2 per cent. 

As an example of Eq. (11), suppose f(z) = h(V — z). The 
result is 


(4/3)h?xt = (1/2){[1 — (s/V)]-? — 1} + (2/V)] (15) 


The exact solution is given by Carslaw and Jaeger (reference 3, 
p. 53) and a plot of this solution is so given (p. 54). To the scale 
with which this graph is plotted, it is impossible to detect any 
error in Eq. (15). 

Eq. (11) may be applied to a nonlinear problem. Suppose 


= (H/k)(z + Vo)4 (16) 
This represents a slab of absolute temperature Vo radiating 
into a medium of absolute temperature zero. The result is 
H2/k?) Vo® = (3/2n8)[(1/6)n? — (3/7)n + (1/4)] + (1/56) 
(17) 
where n=1+2/Vo 
This problem has been solved by Jaeger‘ by integrating the heat 
conduction equation numerically. On the graph which Jaeger 
presents it is impossible to detect any difference between his solu- 
tion and Eq. (17). 
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The Effects of Constant-Pressure Sound Fields 
on Free-Convection Heat Transfer From a 
Horizontal Cylinder 


Jack P. Holman and Tipton P. Mott-Smith 

Fluid Dynamics Research Branch, WADC, Wright-Patterson Air 
Force Base, Ohio 

October 15, 1958 


es QUANTITATIVE INFLUENCE of strong sound fields on free- 
convective heat transfer was investigated by installing a 
3/4 in. electrically heated cylinder in a reverberant chamber 
providing a constant sound-pressure field. The test results are 
shown in Fig. 1. Wall temperatures ranged from 208°F. to 
270°F. Ah represents the increase in the free-convection heat- 
transfer coefficient over the value obtained without the presence 
of the sound field. An increase in the heat-transfer coefficient 
did not occur until the sound-pressure level was greater than 
approximately 135 db. An effect of frequency was observed at 
sound-pressure levels slightly greater than this value; however, 
at higher sound levels the effect of frequency was very small. 
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Fic. 1. Increase of heat-transfer coefficient with sound level 


In an effort to explain the increase in the heat-transfer co- 
effictent when the sound level is raised above a certain critical 
value, a survey of the literature on the effects of sound on the 
flow field about a cylinder was made. References 1—3 show that 
the presence of an external sound field may alter the flow field 
about the cylinder due to a phenomenon called acoustical stream- 
ing. The mathematical relationships that describe the detailed 
behavior of acoustical streaming are based on the following 
assumptions: a plane sound wave is propagated through the 
fluid, the wave length of the sound wave is large compared to the 
radius of the cylinder, and the amplitude of the sound wave is 
small compared to the diameter of the cylinder. The physical 
significance of the various acoustical-streaming parameters 
may be discussed in terms of the typical acoustical-streaming 
flow field shown in Fig. 2. 


CYLINDER 


Fic. 2. Typical acoustical-streaming flow field. 
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Two types of acoustical streaming are usually observed: 
an inner vortex system and an outer vortex system, as shown 
in Fig. 2. The inner system is usually described in terms of 
its thickness corresponding to the distance from the surface of 
the cylinder to the streamline y = 0. This thickness is called 
the DC boundary-layer thickness, dpe. Raney, et al.,? have 
established the existence of a universal thickness profile defined 
by: 


Spc/a = fla/b4c) 


where 64¢ is called the AC boundary-layer thickness and is 


defined by: 
b4c = (v/a) 


where v is the kinematic viscosity of the fluid, w is the frequency 
of the sound wave, and a is the radius of the cylinder. It can be 
shown that a/64¢ is proportional to the ratio of the wave forces 
to the viscous forces. The essential features of the universal 
thickness profile are that 6p)¢ increases with an increase in v 
and decreases with an increase inwora. The detailed acoustical- 
streaming profiles, as determined theoretically, are in excellent 
agreement with experiment as long as the amplitude of the sound 
wave is less than a certain critical value.* It has been shown? 
that if the amplitude is raised above the critical value, the thick- 
ness of the DC layer decreases very rapidly. The critical value 
increases with an increase in the radius of the cylinder. 

The rate at which heat is transported from a cylindrical surface, 
due to free convection, in the presence of a sound field, is de- 
pendent on the superimposed free-convection and sound-field 
currents. In accordance with the above discussion, one would 
not expect any significant change in the heat-transfer rate with an 
increase in the sound intensity until a certain critical value of the 
amplitude of the sound field had been reached. With an increase 
in the sound level above the critical value, the thickness of the 
DC layer would be decreased with a corresponding increase in the 
acoustical-streaming velocity in the boundary layer. Hence, it 
would be expected that the heat-transfer rate would be increased. 
The effect of frequency on the DC layer thickness would be 
observed in the neighborhood of the critical sound-pressure level, 
but at higher sound levels the frequency effect would be very 
small in comparison to the amplitude effect so that the heat- 
transfer rate for sound levels much greater than the critical value 
would be relatively independent of frequency. In addition to 
the direct effects of the increased boundary-laver velocities, the 
acoustical streaming may serve as a trigger mechanism for 
transition to turbulent flow. Consequently, the increase in the 
heat-transfer rate may be dependent on the Grashof Number also 
and, therefore, on the surface temperature difference. This 
dependence will be investigated in future experiments since it 
may comprise an important part of the effects of sound fields on 
heat transfer. 

The sound field in the reverberant chamber used in the above 
heat-transfer tests does not correspond to the plane field which 
has been treated theoretically in references 1-3 since it is a 
scalar field, as opposed to the vector field of a conventional sound 
wave. However, it is expected that the presence of large ampli- 
tude constant-pressure sound fields would have the same general 
effect on the flow field in the vicinity of the cylinder as in the 
case of a plane wave. Therefore, from a qualitative standpoint, 
the analysis of the effects of constant-pressure sound fields on the 
heat-trausfer rate is essentially the same. 
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A New Solution to the Laminar Boundary- 
Layer Equations 


H. A. Hassan 
College of Engineering, Baghdad, Iraq 
October 16, 1958 


N A RECENT WORK, GoOrtler! has presented a new power series 
for solving plane and steady laminar boundary-layer equa- 
tions for a certain class of outer pressure distributions. The ob- 
ject of this note is to present a solution of the same equations 
for a class of outer pressure distributions which has not been 
covered by Gértler. This is accomplished by using the trans- 
formation 
= exp [—(1/») f dx}, = y U,.(x)/v 
g(t, n) = (x, y)/v (1) 
and assuming a power-series solution | U’a(x) is the given outer 
velocity distribution]. The new power-series solution has all 
the desirable. qualities of Gértler’s series and, in addition, its 
zeroth-order term, which is governed by a nonlinear total differ- 
ential equation, can be given in closed form. 
The differential equation for steady two-dimensional laminar 
boundary layer can be written as 
Voy — Vor Voyy = Ueol(x) + v (2) 
with the comma denoting partial differentiation and the prime 
denoting total differentiation. 
The boundary conditions can be written as 


W(x, vy) = v(x, ¥) = ¥,2(x, y) = 0 for y = 0 
V(x, y) U(x) for y— « (3) 


Carrying out the indicated transformation in Eq. (1) and sub- 
stituting in Eq. (2), one obtains 
+ BCE) (1 — Gn?) + — GE Gn] =O (4) 
with 
B(é) = —éEd[In /dé (5) 


where l has the dimensions of velocity and can be chosen as 
the unit velocity or any other representative velocity. 
The solution of Eq. (4) will be assumed as 


0 
in the region of convergence of 
B(é) = ZB, é", Bo > 0 (7) 
0 
Substituting Eqs. (6) and (7) into Eq. (4), the following equa- 
tions for the coefficients ¢, will be obtained: 
+ — = 0 (8) 
and 
+ (n — 2Bo)go'en’ — G0” Gn = Ri-i, n>1 (9) 


with 
n-1 


Ri-1 Bn ( — 1) + (fo — kt 
A 


n-ln-—m 
k=1 m=1k=1 
kee (10) 
k=1 


The boundary conditions for the functions ¢, are 


¢o(0) = (0) = 0, go (2) = 1; 
on(0O) = on'(0) = = O, (11) 


Eqs. (8) and (11) give 
go’ = 3 tanh? { + tanh! 02/3} —2 (12) 
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The functions ¢g, (n > 1) satisfy the linear equations (9) and can 
be expressed as linear combinations of universal functions which 
can be tabulated once and for all for a given value of 8. This 
is accomplished by letting! 

= Bf, = Bi? fir + Be fr, 

= + BiBofis + etc. (13) 
The governing differential equations for the functions f can be 
obtained as follows: letting 


(f) = + (n 280) go’ f’ go” f (14) 
the following differential equations for the functions result [by 
using Eqs. (9), (13) and (14)]: 

1 = Li(fi) => go’? 

n = 2, fr) go? — 1, 

LAfu) = (Bo — 1)fi’ + fifi", ete. 

The boundary conditions for all the functions f are 

S(O) = = =0 (16) 

The class of outer velocity distributions corresponding to Eq. 

(7) can be easily obtained from Eqs. (1), (5), and (7). Eqs. (5) 

and (7) give 


= | — = 


(15) 


Ue~ © exp [ - auer/n U = const. (17) 


n 
It follows from Eqs. (1) that 
dt/dx = (18) 
therefore, Eqs. (17) and (18) give 


x= én (19) 
or? 
t= b, Bo (20) 


Combination of Eqs. (17) and (20) leads to the desired expression 
for Uwo(x)—i.e., 


Uol(x) = x71 co +0 (21) 
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On the Supersonic Flutter of Two-Dimensional 
Infinite Panels 


John W. Miles* and William P. Rodden** 


University of California at Los Angeles, 
and North American Aviation, Inc., Respectively 


October 16, 1958 


iene FLUTTER Of infinite, flat panels on two-dimensional 
periodic supports was studied originally by Hedgepeth, 
Budiansky, and Leonard.! Miles,? comparing their results with 
those for an infinite unsupported panel, questioned several of their 
basic assumptions. Leonard and Hedgepeth® then offered further 
justification of these assumptions on the basis of a “‘static air 
force approximation.’”’ The problem has since been reviewed by 


* Professor of Engineering, University of California at Los Angeles, and 
Consultant, Thompson Ramo Wooldridge Inc. 

** Research Specialist, Vibration and Flutter, North American Aviation, 
Inc., Los Angeles. and Lecturer in Engineering, University of California 
at Los Angeles. 
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Fic. 1. Region of single-degree-of-freedom instability. 


Rodden.‘ The primary purpose of this note is to emphasize 

that much of the apparent conflict among the discussions of 

references 1-3 may be resolved by distinguishing between (a) 

single-degree-of-freedom instability and (8) coupled-mode instability. 

To define our terms, let the panel motion w(r, ¢) be approxi- 
mated by a set of N primitive modes w,(r) according to 

N 
wr, t) = A,(t) wr(r) (1) 


let a, denote the Laplace (or Fourier) transform of a,, and 
let the equations of motion be reduced to the algebraic set 


N 
= O(m =1,..., N) (2) 
1 


by any of the usual methods—e.g., Galerkin’s; the character of 
the permissible motions then will be determined by the N eigen- 
values (usually complex frequencies) for which the determinant 
Of Cm, Vanishes—viz., |Cmn| = 0. If the motion determined by 
the (single-degree) approximation c,, = 0 is unstable, we define 
the instability as type a, while if the primitive motions corre- 
sponding to each of two or more modes are stable but the result- 
ing, coupled motion is unstable, we define the instability as type 
8B. To be sure, these definitions—although mutually exclusive 
are not comprehensive, but we may cover all possibilities in the 
designation (NV, P, R), where P denotes the number of primitive 
modes that are unstable and R the resulting number of unstable 
eigenvalues for the coupled set; types a and @ then are (1, 1, 1) 
and (N, O, R) with N > 2 and R > 1. We remark that the 
coupling of two or more unstable modes may yield only stable 
motions; e.g., (2, 2, 0); but then it seems likely that an improved 
approximation based on a greater number of modes would yield 
unstable motions. We also remark that if N is too small {- 
instability may be predicted where, in fact, it cannot exist. 

The instability considered in reference 2 is associated with a 
harmonic, traveling wave (waves of different wave number being 
uncoupled) and is, therefore, type a. The analysis of reference 
3, on the other hand, rules out a-instability a priori by its neglect 
of aerodynamic damping; accordingly, it can predict only 8- 
instability. 

It is established in reference 4+ that a necessary condition for 
a-instability of the fundamental mode sin(ax/c) is 


(R+12< MP 74241 (3) 
where k = we/xa, M = Mach Number, w» = angular frequency, 
c = support spacing, and a = sonic velocity; the k, \/ region 


described in Eq. (3) is shaded in Fig. 1. We emphasize, however, 
that Eq. (3) is not sufficient for a-instability, since the allowable 
values of k depend on the equations of motion; moreover, the 
linearized aerodynamic theory breaks down as k ~ M — 1 (but 
the equations of motion may bound k away from M — 1 in such 
a way as to prevent this breakdown—e.g., reference 2). Similar 


boundaries may be established for other modes,‘ and we conclude 
that a-instability is especially likely in (although not restricted 
to) 1<M<V2ifk< M-1. 

The analysis of reference 4 achieves an explicit stability bound- 
ary for a two-mode analysis of a periodically supported panel 


tion 
Supe 
6 | 
Mid 


has 
trav 
equi 
sma 
dow 
whe 
able 
is di 
the 
dist: 
equi 
equi 
O 


= Cu 
cS 
are 
> ins 
det 
en 
S 0), 
effi 
loo 
2.0 
che 
Th 
wh 
lus, 
1.3 
sta 
be 
fig 
inv 
tru 
stal 
to 1 
1 
Ana 
of tl 
2 
of t 
Lon; 
ers’ 
Cali 
Har 
Con 
Oct 
iq 
| 
B 


asize 
is of 
(a) 
ality. 
roxi- 


READERS 


by neglecting the (still) air beneath the panel (see Nelson and 
Cunningham’ with regard to this approximation). The results 
are in agreement with those of reference 1 but permit rather more 
insight by virtue of their greater simplicity; in particular, they 
demonstrate that the second stability boundaries of Fig. 9, refer- 
ence 1, correspond to transition to a stable motion of type (2, 2, 
(0), thereby affording a probable explanation for the destabilizing 
effect of structural damping (but, as noted above, these second 
loops are not likely to be significant approximations). 

An analysis of a two-bay continuous panel at M = 1.56 and 
2.00 also is presented in reference 4 and tends to support the 
choice of modes for the periodically supported panel in reference 
1, the spectra of the two configurations being somewhat similar. 
The mean stiffness requirement found at 7 = 1.56 and 2.00 was 


— v)] (E/q)} (t/c) = 0.5156 (4) 


where B = (1M? — 1)!/2, » = Poisson’s ratio, E = Young’s modu- 
lus, g = dynamic pressure, and ¢ = panel thickness. At WV = 
1.3 and V2, on the other hand, a@-instability was found and no 
stability boundary obtained. The detailed analysis of this con- 
figuration will be published following the completion of a similar 
investigation for a four-bay panel, for which the frequency spec- 
trum more closely resembles that of the periodically supported 
panel. 

We conclude that the problem of single degree of freedom in- 
stability warrants further investigation, especially with respect 
to the approximation of two-dimensional flow.® 
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A Note on Blast-Wave Theory 


Harvey Einbinder 
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October 21, 1958 


— THEORY, first developed by G. I. Taylor! to de- 
scribe the shock wave produced by an atomic explosion, 
has been applied to the flow about slender blunt-nosed bodies 
traveling at hypersonic velocities. As Hayes? has shown, the 
equations of this theory are equivalent to those of hypersonic 
small-disturbance theory. The theory is valid for distances 
downstream that are not too close to the stagnation region but 
where the strong-shock conditions hold. It is, therefore, suit- 
able for analyzing the inviscid flow that occurs outside the 
boundary layer along the body. The simplicity of the theory 
is due to the assumption of flow similarity, which requires that 
the variation in the transverse flow field be independent of the 
distance downstream. When this assumption is satisfied, the 
equations of motion can be reduced to ordinary differential 
equations, whose solution is much simpler than the more general 
equations. 

Of special interest is the case of constant energy release which 
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corresponds to a blunt-nosed body of constant drag, such as a 
hemispherically-capped cylinder, where the drag of the body 
provides the transverse velocity field. For this case, the equa- 
tions of motion were first integrated numerically by Lin.* Re- 
cently they have been integrated analytically by Kubota,‘ by 
following the procedure of Latter, who solved the equations 
for the spherical blast wave. 

The purpose of this note is to draw attention to some of the 
consequences that can be derived from Kubota’s analytic solu- 
tion. According to his solution for the axially symmetric case, 


@ = sz (1) 
F = | — 2)” ~ "J — 
where DAy + 1) (2 
= {ev — — o)g?} 1/8 


where = 4(y — t+! (3) 


where ¢, F, and y are dimensionless variables that represent the 
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transverse variation of the velocity, pressure, and density, 
respectively. y = C,/C,, s is an auxiliary variable, and z = 
y/R, where y is the transverse distance from the axis of the body 
and R is the distance to the shock front. When 0 < z < 1, 
then 1/y <s < 2/(y +1). The variable s is defined implicitly 
by the equation 


Vilys — 1)/on) (4) 
where ¢ = 4/[‘y + 1)*s(2 — ys)] and om = (y — 1) + 
(y + 1) 
The expressions for F and y can be simplified by expressing them 
in terms of s by eliminating ¢ and z with the aid of Eqs. (1) and 
(4). Thus, 


F = s{l(v — 1) (2 — ys)I/[21 (5) 
G = (F/s)?/7 [(ys — (6) 


where G = omy, so that at the shock front where z = 1,G = I’ 
and F = 2/(y + 1). These functions are shown in Figs. 1 and 
2 for y = 1.40 and y = 1.18, corresponding to normal and dis- 
sociated air, respectively. It can be seen from these figures that 
there is a slow decrease in the pressure, F, in going from the shock 
toward the body. In contrast, there is an exceedingly rapid de- 
crease in G, the density. The form of this decrease can be ex- 
hibited explicitly by means of an approximate formula for G 
valid when z < 0.8. 

Whenz = 1, Vi = 1, while whens = 0, Ve = IV y/(y +1) 
= 1. Therefore, Ve can be replaced by unity in Eq. (4). Using 
this equation to eliminate (ys — 1)/o, in Eq. (6), and replacing 
F/s by its value at z = 0, vields the desired expression for G. 

G = (7/2)/2 (g £08) (7) 


Now the exponent of z in this equation, 2/(y + 1), will be large. 
For y = 1.40, it equals 5, while for hypersonic flow where y is 
closer to unity, it will be even larger. Thus, when y = 1.18, G 
varies as the eleventh power of z. This extreme variation in the 
density between the shock front and the body does not seem 
plausible. Furthermore, as y approaches unity, the exponent 
tends toward infinity so that the solution breaks down. 

The rapid decrease in the density implies a corresponding in- 
crease in the temperature since by the equation state, 7 ~ P/p ~ 
F/G. Thus, according to blast-wave theory there is a very low- 
density, high-temperature layer in the inviscid flow just outside 
the boundary layer. 

The solution for the flow field which has been obtained is based 
on a constant y, which is reasonable as long as the flow does not 
depend critically on the value of y. However, according to the 
theory, the density and the temperature are very sensitive to the 
value of y in the hypersonic regime. Therefore, it is not justified 
to neglect its variation with temperature. 


Kubota* ® has attempted to confirm blast-wave theory by 
measuring the shape of the shock front and the pressure along 
the surface of a hemispherically-capped cylinder for M = 7.7, 
However, these quantities are relatively insensitive predictions 
of the theory. Consequently, such measurements do not pro- 
vide a very adequate test of the theory. 

Because of the complex interaction that occurs between the 
body and the shock front, blast-wave theory fails to provide an 
adequate description of the hypersonic flow about biunt-nosed 
slender bodies. The presence of such interaction is demon- 
strated by the characteristic lines which are reflected a number of 
times between the body and the shock front in passing from the 
stagnation region downstream. This is the primary cause of the 
failure of blast-wave theory. 
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Errata—‘‘Some Panel-Flutter Studies Using 
Piston Theory’’* 


David J. Johns 

Department of Aircraft Design, The College of Aeronautics, 
Cranfield, Bletchley, Bucks., England 

January 20, 1959 


I WOULD LIKE to point out an omission and an error in my 


recent paper. 
On page 681, the definition of p,(x) should be: p,(x) = 
Pi, pe pew? with Z replaced by Z, in Pr. 
On page 684, the values of 2gc*/./D in Fig. 5 should be factored 
by (800/3/9). 


* Johns, David John, Some Panel-Flutter Studies Using Piston Theory, 
Journal of the Aero Space Sciences, Vol. 25, No. 11, pp. 679-684, 715, 
November, 1958 
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